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PROLOGUE 


This report summarizes the major points of a doctoral dissertation 
submitted in partial fulfillment of the requirements for the Doctor of 
Philosophy degree, in the Horace H. Rackham School of Graduate Studies 
at The University of Michigan. Copies of the dissertation are available 
for study at the Library of Congress and the University Library. 
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CHAPTER I. INTRODUCTION 


1.1 Background . 

1.1.1 So lar Activity . The sun is a restless star. Its activity 
is manifested by multiple characteristics: it may be as emission of 

particles or of electromagnetic waves, it may vary slowly or abruptly and 
violently, it may be very weak or extremely strong. 

We can observe the efflux of matter from the Sun that ranges from 
the steady solar wind up to the explosively ejected packets of high- 
energy particles. Also from the Sun we can observe the widest electro- 
magnetic spectrum of any celestial object. It spans at least fourteen 
decades of wavelengths that range from about the 10 8 cm of hard X-rays 
up to 10 s cm of the Type III burst radio emission. 

These forms of activity exhibit different degrees of interdependence. 

The solar radio emission can be classified according to frequency 
range and spectrum, to temporal characteristics, to polarization, or to 
morphology. Commonly a mixture of these elements is used since they are 
not completely related. 

The broadest classification is based on intensity-time behavior: 

(l) The emission of the undisturbed Sun and (2) the bursts. The first is 
divided into two parts : A background component which remains at a constant 

level for months or years, and a slowly varying component which changes 
gradually as the active centers are carried by the Sun's rotation. The 
steady background is sometimes called the "quiet sun" and it corresponds 
to the thermal radiation of the Sun, it excludes the active centers or 
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other sources of transient perturbations; it changes with the solar cycle. 
The slowly varying component, or S-component, is also of thermal origin 
but it is associated with sunspots and chromospheric plages or active 
centers. Sometimes the designation of quiet sun seems to be used in a 
.more general context in place of undisturbed sun. This radiation has 
been observed in the range between millimeter and meter wavelengths and 
it originates approximately between the lower chromosphere and one solar 
radius above it. 

The second group in the broad classification of solar radio emission, 
the radio bursts, refers to the active sun. They are observed between 
millimeter and kilometer wavelengths and their characteristics depend on 
the region of the frequency spectrum where they are found; this in turn 
is related to the dominating mechanism of generation. In contrast with 
the quiet sun the bursts are primarily of non-thermal origin and they 
are usually associated with explosive events like flares. 

The original burst classification started with three groups: the 

types I, II and III that were followed later by the types IV and V. At 
first all these types were observed in the meter range but later studies 
have revealed that the Type IV is present in the centimeter range while 
the Type III have now been observed at wavelengths in the kilometric range. 
Other bursts of a more uncertain classification are also observed in the 
centimeter range. This gave rise to a more loose grouping in terms of 
wavelength range (microwave, meter wave, etc.). 

Review articles on solar bursts have been published, for example by 
Wild, Smerd and Weiss ( 196 '-;), Maxwell (1969) and Boischot (1067a). The 
whole subject of solar radio emission is covered at an elementary level by 
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A. g. Smith (1967) and exhaustively at an advanced level in books by 
Kundu (1965) and by Zheleznyakov ( 1970 ). 

The radio bursts occur in the mass of gas that surround the Sun. 
Several names for this gas are found in the literature for example solar 
atmosphere, solar corona, solar plasma and solar wind; we will use them 
indifferently. Sometimes the term corona or inner corona is used to 
designate the part within probably ten solar radii while the terms outer 
or super corona refer to the part beyond that approximate limit. 

The important characteristics of the solar bursts are: (a) the 

frequency range in which they are observed during their lifetime, (b) the 
instantaneous bandwidth, that is, the frequency range in which they are 
observed at any instant, (c) the spectral structure that shows whether 
the emission is concentrated in certain frequencies (harmonics) or in a 
continuous range, (d) the type and degree of polarization of the radiation, 
(e) the time of start at the different frequencies, (f) the time structure 
at the different frequencies, (g) the peak intensity, (h) the size, position 
and motion of the emitting region. 

Which characteristics can be measured depends on the observing 
instrument. Basically the receiver may be tuned to one or more fixed 
frequencies or its bandpass may be swept rapidly and continuously in 
frequency. In the first case, we obtain bursts profiles (plots of intensity 
versus time), in the second, we obtain dynamic radio spectrograms (plots 
of frequency versus time, modulated by intensity). The antenna may be 
a simple antenna or an array. A large array is required to measure size 
and position of the source. 

The sequence of Roman numerals used to designate burst types does 
not bear any generic relationship, however Types IV and V are usually 
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observed after Types II and Irrespectively. While I, II and IV have a 
well defined spectral structure, IV and V belong to the continuum class. 

Bursts of Type I are short lived (tenths of seconds), occur in rapid 
succession and are superimposed on a more or less steady background that 
may last for hours or days. The bursts themselves are confined to a 
narrow band of frequencies while the background emission is a broad 
continuum. The combination of the two phenomena has been called "noise 
storm." They are often 100$ circularly polarized and their emission 
mechanism is not yet definitely identified. They are not associated with 
flares but with active centers. Their existence has been reported at 
frequencies as low as 0.2 MHz (Fainberg and Stone, 1971b). 

A very useful concept in dealing with Type II and III bursts is 
that of frequency drift rate (FDR). It is the rate of change of burst 
arrival time at the different frequencies and it is expressed in MHz sec -1 . 
In the classical Types II and III bursts the FDR is negative, that is, 
the lower the frequency the later the time of arrival. In what follows 
this sign will be implicitly understood. 

This driftrate is a direct consequence of the emission mechanism of 
these kinds of bursts: a perturbation originating on the sun that travels 

outward through the solar atmosphere producing radio waves at frequencies 
related to the fundamental or the second harmonic of the local plasma 
frequency. Since the density of the solar plasma decreases with increasing 
distance the bursts exhibit a negative frequency drift rate. The magnitude 
of the FDR depends, among other things, on the velocity of the perturbation 
and on the gradient of electron density distribution along the path of the 
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The Type II bursts are characterized by a low FDR, the highest 
being of the order of 1 MHz sec at about 200 MHz. For this reason 
they have also been called "slow drift" bursts. Their detection has been 
claimed at frequencies as low as 0.2 MHz (Slysh, 1967c) but without 
convincing supporting facts. The exciting perturbation is believed to 
be plasma shock waves emanating from flares and with speeds of a few 
thousand km sec . The spectral characteristic is a well-defined frequency 
structure prominent at the fundamental and, often, at the second harmonic. 
The radiation is unpolarized or very weakly polarized. In the meter range 
they last for about 10 minutes. 

lyP e HI Bursts . These are the most common of all the solar 
phenomena. They constitute the only type that is observed at frequencies 
as low as 0.050 MHz covering a range of at least four decades. Haddock and 
Alvarez, (1969) . Their principal characteristic is a high FDR in the 
meter range, of about 30 MHz sec -1 at J00 MHz, decreasing with increasing 
wavelength. They are called "fast drift" bursts. 

Since the Type III bursts span a large frequency range it is not 
surprising that their characteristics vary drastically with frequency. 

For example, between a few hundred MHz and 0.050 MHz the duration of the 
bursts increases from tenths of a second to thousands of seconds, the 
intensity of one burst may increase by a factor of 10 5 or 10 s , the frequency 
drift rate decreases from approximately 10 2 to 10" 4 MHz sec _1 . 

The first report in the scientific literature on solar radio 
emission is that of Reber (19^). The discovery of solar radio emission 
was made during World War II by Hey in 19^-2 but for military security 
reasons it was not reported until later. Hey (19^-6). Also in the same year 
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and independently Southworth observed quiet sun radio emission that he 
did not report until 1945. Before this many reports were made during 
the sunspot maximum of the previous decade by ham operators, in Q.S.T., 
etc. The history of the pre-discovery period has been summarized by 
A. G. Smith (1967). During the rest of that decade many papers were 
published; for example by Appleton (19^5), Shklovskii (1946), Pawsey, 
Payne-Scott and McCready (1946), Denisse (1946), Khaikin and Chikachev 
(1947)> Allen (1947b), Hey, Parsons and Phillips (1948), Payne-Scott 
(1949), etc. These workers recognized the existence of sudden short 
enhancements that were called bursts or outbursts. Of special importance 
was the observation of Payne-Scott, Yabsley and Bolton (1947), confirming 
Bowen's verbal report that there were "successive delays between the onset 
of the burst on 200, 100 and Mc/s..." 

The decade of the 50 's was very prolific in solar work. This was 
carried out especially by the Australian scientists and it started with 
a series of fundamental papers by Wild and McCready (1950) and Wild (1950a, 
1950b). The first paper defines the Types I, II and III; the second is 
related mainly to Type II and the third to Type III. They constitute the 
first comprehensive study of solar bursts. These studies were performed 
between 70 and 150 MHz that is a very small fraction of the frequency range 
covered by these bursts; nevertheless many of the approximate formulas 
found by Wild to represent some of their characteristics have been shown 
by us to be valid in the kilometer wavelength range. The success of Wild's 
observations was due to a large extent to the introduction of the 
"panoramic receiver" or dynamic spectrograph. 

The discovery of the delays in the arrivals of a burst at different 
frequencies aroused interest in the mechanism of excitation. Payne-Scott, 



Yabsley and Bolton (19^7 ) wrote "the successive delays . . . suggest that 
the outburst was related to some physical agency passing from high-frequency 
to lower-frequency levels. If we assume, following ideas suggested by 
Martyn, that the radiation at any frequency originates near the level 
where the coronal density reduces the refractive index to zero, we can 
derive a rough estimate of these heights from electron-density data given 
by Baumbach." They derived velocities between 500 and 750 Km sec 1 and 
they noticed that they "are of the same order as those of prominence 
material . . . and auroral particles ..." Because of the magnitude of the 
speeds we can infer that they observed Type II bursts. 

It is interesting to notice that earlier Shklovskii (19^6), 
apparently without knowing about the observed frequency delays, had 
proposed that "these proper oscillations can be excited by the streams of 
charged particles moving through the corona with a velocity larger than 
the velocity of sound." This paper seems to have been unknown to Payne- 
Scott, Yabsley and Bolton and has remained almost unnoticed. 

Wild (1950b) examined these ideas and concluded that "the explanation 
of the frequency drift of Type III bursts in terms of the motion of exciting 
agencies is consistent with observational evidence. With the assumed 
electron density . distribution, [Baumbach-Allen ] velocities of between 
about 2 x 10 4 and 10 5 Km sec 1 would be required to account for the 
observed drift rates." 

The presence of double peaks on single-frequency observations 
was recognized as early as 19^9 by Payne-Scott, and it was interpreted as 
due to an echo of the original burst (from the characteristics of the 
data we have inferred that she observed Type III bursts). This was probably 
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the same phenomenon that Wild, Roberts and Murray (1954) interpreted later 
as harmonic pairs in Type III bursts. It was to them the justification 
of the application of the critical frequency hypothesis to this kind of 
bursts. This hypothesis has been called the Local Plasma Hypothesis and 
we have accepted it in this work. 

The second harmonic has had important implications in the theoretical 
studies of the emission mechanisms. One of the main results of this work 
is the unexpected discovery that the relative amount of second harmonic 
radiation observed at the Earth increases as the frequency decreases. We 
have called this the second harmonic phenomenon or second harmonic prevalence. 

Two other mechanisms have been proposed besides the local plasma 
hypothesis. The first, by Slysh (1967), postulates that the bursts are 
emitted at a frequency that is higher than the local plasma frequency and 
that it is sensitive to the local temperature. Slysh’ s paper contains 
some errors that were partially corrected by Papagiannis (1971). A version 
completely corrected by us still is inconsistent with our results. Our 
discovery of the second harmonic phenomenon certainly affects Slysh' s 
argument against the local plasma hypothesis. 

A second alternative to the local plasma hypothesis has been 
offered by Kuckes and Sudan (1969, 1971). They propose that the emission 
is related to coherent synchrotron deceleration of 100-keV electrons that 
excite radio waves at the electron cyclotron frequency and its harmonics. 

A possible weakness of this theory rests in the magnitude of the required 
magnetic field, that may exceed the observed values. The second harmonic 
phenomenon found by us, that places the origin of most hectometric and 
kilometric wavelength bursts further out than was previously thought, may 
also effect this theory. 
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According to the local plasma hypothesis fast particles excite 
plasma waves that in turn produce radio waves . An important and 
influential theoretical paper is that of Ginzburg and Zheleznyakov ( 1958a). 
According to these authors the plasma waves are produced by a Cherenkov- 
type process and are amplified by the two -stream instability resulting 
from the interaction of the fast charged particles with the ambient 
plasma. The mechanism by which the longitudinal waves couple into 
transverse waves with the observed characteristics is not completely 
understood. The same theoreticians proposed that the conversion of energy 
from plasma to electromagnetic waves is achieved by the scattering of 
'coherent" plasma waves on small-scale inhomogeneities of the background 
plasma. The slow thermal motion of the heavy ions produces quasi-neutral 
random fluctuations SN ' of the plasma density while the light and mobile 
electrons produce spatial electric charge fluctuations BN". "Since the 
variation of plasma density occurs rather slowly, scattering by BN' does 
not result in any significant change in the wave frequency ... so that 
the electromagnetic waves formed by this scattering have the same frequency 
as the incident plasma wave. This scattering process is called Rayleigh 
scattering. The radiation in this case is emitted incoherently and mainly 
in the forward direction" (Kundu 1965, p. 59). The space-charge fluctuations 
BN" of the electron density can be plasma waves present in the background 
plasma with frequencies close to the plasma frequency, waves with 
other frequencies are more strongly damped. The interaction of coherent 
plasma waves set up by the exciting particles with the space-charge 
fluctuations produces an electromagnetic wave at twice the plasma frequency. 
This process is called combination scattering. According to Ginzburg and 
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Zheleznyakov, the fundamental of Type III bursts is generated by Rayleigh 
scattering and the second harmonic by combination scattering. This 
theory does not require a magnetic field. But is was considered later by 
them in connection with noncoherent mechanisms (Ginzburg and Zheleznyiakov, 
1961). They conclude, among other things, that solar bursts of Types II 
and III may not be related to the "noncoherent" emission of plasma waves 
in an isotropic plasma. The conversion of plasma waves into electromagnetic 
waves in a strong magnetic field was further studied by Tsytovich and 
Kaplan (1969). 

Type III bursts are elliptically polarized. This phenomenon has 
been treated theoretically by Yip (1970c). 

During the last few years several theoretical contributions have 
been made, for example, by Yip (1970a, b and c), Melrose (1970), 

D. F. Smith (1969, 1970), Zheleznyakov and Zaitsev (1970a and b). All 
these authors adopt scattering as the basic process of energy conversion. 

A radical departure from these ideas was taken by Kuckes and Sudan (1971) 
who invoke coherent synchrotron deceleration of high-energy electrons. 

Reviews on the mechanisms of solar radio bursts have been given 
by Wild, Smerd and Weiss (1963), and by Takakura (1966, 1967). 

As a direct consequence of the local plasma hypothesis it is 
expected that the relaxation of the plasma oscillations set up by the 
passage of the exciter particles is governed by the physical condition 
of the local plasma. For this reason the study of the time of decay of 
Type III bursts is of special interest. 

The importance of Type III bursts lies in the fact that they can 
be used to study "in situ" some properties of the solar plasma at different 



- 11 - 


levels. Except for deep-space probes this may be one of the few methods 
which can do this. Other radio methods to study the solar plasma involve 
the observations of the scattering of radio waves from extra-solar radio 
sources, Hewish (1955)., Baseylan and Sinitsin (1971) or from space probes, 
Stelzried et al. (1970)* Also other techniques have been used, radar 
observations, James (1970), and the observation of dispersion of pulsar 
pulses, Hollweg (1968), Goldstein and Meisel (1969). The usefulness of 
these methods has been limited in determining the electron density 
distribution because they give only integrated densities along the line 
of the sight. Some of the techniques can give information about the size 
of density inhomogeneities. 

Another radio method of studying solar plasma, recently proposed 
by Ferencz and Tarcsai (1970) (Consists- of the observation of an anomalous 
frequency decrease of the 21-cm emission of Taurus A when it passes close 
to the Sun. 

The earliest radio method was suggested by Vitkevich in 1951- It 
consists of the observation of the occultation of a radio source by the 
solar corona. The method would allow the determination of electron density 
electron temperature and magnetic fields of the corona. The observation 
of changes in the apparent angular size of the radio source would permit 
the study of coronal inhomogeneities. Later that year, and independently 
of Vitkevich, Machin and Smith (1951) proposed a similar but less 
comprehensive method. 

The Type III bursts can be used as "depth probes" because the 
frequency of observation determines the electron density of the emission 
level, provided that the mode is known to be fundamental or second harmonic 
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The height of the level above the photosphere is undetermined; however, 
because of the negative density gradient, the lower the frequency the 
higher the level. 

The early observations of Type III bursts were performed in the 
meter-wavelength range. This corresponds to sources a few solar radii from the 
Sun. Interferometric observations made by Wild, Sheridan and Neylan 
(1959) showed that the exciters velocities ranged between 0.2 and 0.8 of 
the velocity of light in vacuum, and that the deduced electron densities 
were about one order of magnitude larger than the standard Baumbach-Allen 
model. Similar results by other workers led to the idea that the bursts 
are produced in dense streamers. 

High-frequency observations of the early workers revealed also 
that the time decay of the bursts was nearly exponential, Williams, (19U6) ; 
Payne-Scott^ (1949) . The last author was aware that this phenomenon could 
lead to determination of the temperature of the solar plasma, however she 
only mentioned it. As far as we know the first calculations of temperature 
using this property were carried out by Boischot, Lee and Warwick (i 960 ). 

The values derived by them were in agreement with those determined by 
optical means. 

Observations from Spacecraft . Because of the shielding 
effects of the Earth ionosphere on the waves at frequencies below 
approximately 5 MHz, the ground-based radio observations of Type III 
bursts yielded information from levels not higher than that reached by 
optical means. The advent of Earth artificial satellites permitted 
observations at longer wavelengths. The solar radio observations from 
space were initiated by Hartz (1964) who obtained information down to 
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1.5 MHz. Since then lower and lower frequencies have been reached. Table 
1.1 shows a complete list of solar radio astronomy instruments aboard 
artificial satellites or space probes and some of their characteristics. 

They are presented in order of decreasing lowest frequency of observation. 

Very low frequency observations from spacecraft have to overcome 
many problems. Among them we should mention the low directivity attainable 
with electrically short antennas and its impedance matching over a wide 
range of frequencies. The matching difficulties result from the change 
in the antenna electrical parameters with electrical length as well as 
from the modifications of these parameters by the surrounding plasma. 

The extension of the observations to successively lower frequencies 
has allowed the possibility to probe further and further out into the 
solar corona. Our discovery of the second harmonic phenomenon about doubles 
the distance associated with Type III bursts below approximately 1 MHz. 

Up to the advent of space observations, and particularly by OGO-V, there 
was a large distance gap that extended between 20 to 215 solar radii between 
the realm of the optical observations and that of the Earth neighborhood 
spacecraft observations. This gap was filled by the radio observations 
of Type III bursts from satellites. 

As the lowest frequency of observation was lowered some discrepancies 
with the measurements made of the solar atmosphere by other means began to 
stand out. Ground-based radio observations indicated that the electron 
densities along the exciters' path was about one order of magnitude larger 
than predicted by the current solar wind models. This was explained by 
assuming that the sources of the bursts occurred within dense streamers. 

The existence of these streamers at distances larger than about 100 solar 
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radii had to be invoked in order to interpret the very low frequency 
bursts in the usual way. Our interpretation based on the second harmonic 
phenomenon has removed this condition for distances larger than a few 
tens of solar radii at the same time that it maintains the constancy of 
the velocity of the exciters at values in agreement with those obtained 
from direct observations. 

Another important discrepancy refers to the temperature of the 
solar corona. The values deduced from the decay time of the bursts 
observed by Hartz in 196^, assuming collisional damping, were lower than 
those expected from solar wind models. This discrepancy increased as 
lower frequency data became available. The closing of the gap between 
the Sun and the Earth has permitted us to demonstrate quantitatively that 
this temperature discrepancy is the same that solar wind theoreticians have 
found in the last few years between theory and experiment near the Earth. 
The reason for this discrepancy is still unexplained. 

Low-frequency observations have also allowed Fainberg and' Stone 
(1971a) tomeasure inhomogeneities in the ambient plasma^and the solar 
wind speed. 

Observations in this range have brought the need to consider the 
shape of the exciters 1 path and its orientation with respect to the Sun- 
Earth direction. Ground-based measurements studied levels so close to the 
Sun that it was a good assumption to make that the path of the common 
Type III burst exciters is radial. In working at distances larger than 
100 solar radii the assumption may not be justified. Though we have no 
definite evidence against the radial trajectory we have chosen a spiral 
path determined by the shape imparted to the interplanetary magnetic field 
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by the Sun's rotation. Preliminary data from IMP -VI look very promising 
to elucidate this and other important aspects of Type III bursts, such 
as the size of the sour ces ; Haddock (1971). 

The very low frequency observations will allow the study of the 
associations between the Type III exciters and the high-energy particles that 

are being ejected from flares and passing by the Earth. In this way the 
controversial identity of the exciters may be revealed. 

Finally, since the Earth is affected so directly by the Sun's 
activity, the study of Type III bursts may be of practical importance in 
obtaining a better understanding of the solar-terrestrial relationships. 

1.2 Outline of the Research 

The data of the OGO-V satellite give plots of intensity versus time 
at the eight frequencies of observation and with a time resolution of 
9 seconds. Of the many bursts observed we selected 6.4 that 
drifted down to and below 0.55 MHz. From them we measured: 

(a) Time of occurrence, 

(b) Amplitude and 

(c) Time of decay. 

The basic assumptions in the analysis of these data are: 

(A) Local plasma hypothesis, 

(b) The exciter particles travel with a constant speed along an 
Archimedes spiral situated on the plane of the ecliptic, 

(C ) The radio waves travel from the source to the Earth as in vacuum, 

(D) The fundamental and second harmonic emission are generated 
simultaneously and in the same small region and 
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(E) Except for their role in guiding the exciter particles we do not 
consider the effect of magnetic fields. 

In Chapter 2 we describe the instrumentation of the OGO-V radio 
astronomy experiment. We discuss the receiver preflight and inflight 
characteristics. Attention is given to the study of the spacecraft-antenna 
system in order to perform approximate analysis of the burst intensities. 

We examine the data acquisition system and how, together with the receiver 
transfer characteristics, it affects the large signals that will be of 
special importance to us. At the end of this chapter we describe the 
satellite orbit and discuss the effect of orbit and spacecraft motion on 
the orientation of the antenna pattern. 

Chapter 3 deals with the presentation and characteristics of the 
data. These, plus scientific considerations, .lead us to establish a 
selection criterion. We emphasize throughout the whole Work that the 
validity of some of our results may be limited by this selection. 

Chapter ^ is devoted to the data analysis proper and it is divided 
into five parts. In the first part we discuss what can be measured and 
learned from the data supplied by our OGO-V experiment. The second part 
presents the analysis of the bursts' start time. The delays in the arrivals 
at the different frequencies are expressed as drift rates. We combine the 
results from OGO-V bursts with as many results as we could find in the 
literature and we obtain a simple law of FDR versus frequency valid at 
least between 550 and 0.075 MHz. By using this law, the basic assumptions 
and the electron densities measured near the Sun and near the Earth we 
derive an expression of electron density versus distance. This is used as 
an intermediate step to arrive at a simple semi -empirical formula in 



-18 


agreement with the electron densities observed near the Sun and near the 
Earth, consistent with the observation of drift rates of Type III bursts 
and that is not in conflict with theoretical solar wind models. Other 
applications of the formula are discussed. 

We describe a new technique devised to study the bursts ' time of 
arrival. This consists in predicting the arrival times at the Earth, for 
the fundamental and second harmonic emission. A comparison with the 
observations shows unexpectedly that at approximately 1 MHz tHe mode 
of the observed radiation changes from predominately fundamental to 
predominately second harmonic. This preponderance of second harmonic 
increases with decreasing frequency and at 0.05 MHz most of the emission 
appears to be harmonic. Extrapolation of these results to higher frequencies 
agrees with ground-based data. Next, we make direct use of this phenomenon 
to determine velocities of the exciters. By assuming a density model as 
thin as the solar wind we find that we can deduce values similar to those 
obtained from direct ground-based observations. We discuss the effects of 
the basic assumptions on the derived velocities. 

The third part of Chapter 4 contains the study of the burst amplitudes. 
The absolute quantitative aspect relies on the validity of the antenna 
model discussed in Chapter 2. The power spectrum is usually double* or 
single peaked. We suspect that the double peak may be related to the 
second harmonic emission. The power levelof the burst maximum is found to 
be several orders of magnitude larger than at frequencies observed from 
ground. The power spectrum may peak at frequencies as low as 0.05 MHz. We 
study also energy spectra and find results similar to those obtained for 


the power spectra. 
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We approximate an average energy spectrum by a rather simple 
mathematical formula that, at the high frequencies, reduces to that found 
by Wild in 1950. We discuss possible implications of our results on the 
theory of the mechanisms of Type III burst emission. 

In the fourth part of Chapter 4 we analyze the times of decay. We 
explain the technique used and its limitations. Combining our observations 
with all of those we could find at other frequencies we obtain a simple 
formula to express decay times versus frequency. We compare the observations 
with the Coulomb-collisional decay times computed from theoretical solar 
wind models and conclude the complete failure of the simple collisional 
damping hypothesis to explain the decay time of Type III bursts. Combining 
observations and theory we derive an approximate relationship for the ratio 
of observed to computed decay times as a function of distance. Because this 
function is derived from radio observations it is valid in a continuous 
range of distances. This relationship is equivalent to the one that solar 
wind theoreticians found while fitting the theory to temperature observations 
near the Earth. When these relationships are expressed in terms of the 
same parameters they become practically equal. 

In the last part of Chapter 4 we examine the association of the 
observed kilome trie-wave bursts with other solar phenomena: Ha and 

X-ray flares, energetic particles and radio bursts detected from ground. 

By assuming that the origin of the bursts coincides with the center of the 
Ha flares we find that the distribution of the kilometric-wave bursts 
is heavily biased toward the western solar hemisphere while the distribution 
found by ground based observers at higher frequencies is the opposite. Using 
unpublished data by Lin (1971) we find that all the kilometric-wave 
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bursts associated with high-energy .. solar electrons also occurred on 
the western solar hemisphere as do our Type III bursts. 

Chapter 5 contains a summary of the results achieved in this work. 
In Chapter 6 we make recommendations for future experiments. 

In the course of this research we concentrated not only on the 
properties of the kilometric-wave . bursts but also we tried to combine 
them with the results of other workers in as many frequencies as it was 
possible to find. This turned out to be a rewarding approach because it 
unveiled general trends across four frequency decades. Often this allowed 
us to write fairly approximate but simple formulas to describe the 
characteristics in question. 

Numbers of equations , tables or figures included in the appendices 
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CHAPTER II. THE OGO-V RADIO ASTRONOMY EXPERIMENT 

2.1 Description of the Instrumentation 

The radiometer consists of a 9 » 15 -m monopole antenna and a 
stepping- frequency superheterodyne receiver tunable through the eight 
frequencies of 0.050, 0.100, 0.200, 0.350, 0.600, 0.9 00, 1.800 and 
3.500 MHz in 9.2 seconds. This instrument has been described in detail 
by MacRae, (1968) and will only be described briefly below. 

The Radio Astronomy Experiment package is located in the Solar 
Oriented Experiment Package No. 1 (SOEP-l) of the Orbital Geophysical 
Observatory V ( 0 G 0 -V) spacecraft with the monopole antenna extended in the 
-X direction that is perpendicular to the Earth-space craft -Sun plane. 

(See Fig. 2 . 1 . ) 

2 . 1.1 The Antenna . The antenna unit consists of a 5 * 08 -cm wide 
and 9.15-m long beryllium- copper ribbon wound on a spool before deployment. 
The beryllium- copper ribbon is heat-treated in such a way that in the 
free state it forms a cylindrical tube of approximately 1.27 cm in 
diameter. A ground command releases a restraining strap and allows 
the . ribbon to unwind (by virtue of the stored energy in the wound 
condition) and to form a cylindrical boom. A plating of silver on top 
of gold is deposited on the outer surface of the beryllium-copper strip 
to improve its thermal properties. The inboard 3.05 m. of the ribbon are 
coated with Dupont "Kapton" to insulate the boom from the plasma sheath 
surrounding the spacecraft. 
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The electrical length of the boom ranges between 0.017 and 
0.0015 wavelengths at the frequencies of operation. 

The antenna equivalent circuit used in the engineering phase of 

the work is shown in Fig. 2.2 where 

= antenna capacitance, 

R = antenna radiation resistance, 

C = equivalent stray capacitance, 

B 

V 2 = 1+ k T R Af = mean squared noise voltage, 

T = antenna temperature , 

Af = bandwidth and 

k = Boltzmann's constant, 1.38 x 10 23 joule K 1 . 

The determination of the exact values of these parameters is 
not possible because of the complex configuration of the monopole- 
spacecraft structure, and also because the spacecraft is not in vacuum. 
However estimates can be obtained by adopting models for the structure 
and by assuming it to be in vacuum. 

Because of the shape of the spacecraft we considered two simple 
models. In the first we assumed that the Radio Astronomy antenna is 
fed against an infinite and perfectly conducting plane. In the second 
we considered that our antenna and the rest of the spacecraft behave 
like a center- fed short dipole. In spite of the fact that neither 
structure is an exact representation we felt that the above approxi- 
mations probably represent limits to the actual situation. From them 
we computed the parameters involved and found that the monopole and 
dipole results differ at most by a factor of two. We believe that the 
dipole model is better therefore we adopted It. 
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The radiation resistance of an electrically short dipole is 
given by the well known formula: 



where h = dipole half-length and 

A = free space wavelength. 

The capacitance derived from transmission line analogy and valid 
for an electrically short antenna is (Jordan, 195Q, p. 1*65): 

( 2 . 2 ) 

where £ q = permittivity of free space and 

a = antenna radius . 

The other parameters are the same as in Eq. (2.1 ). 

In Fig. 2.2, dg is the stray capacitance due to the proximity 
to the antenna of grounded surfaces and due to the connecting leads. 
Careful preflight measurements of the configuration gave a value of 
7-5 P F for the base capacitance of the antenna unit itself. The 
contribution of the connecting leads was estimated as 7.5 pF 
that gives =15 pF. 

E • A * 2 The He cei ver ■ As mentioned earlier, the receiver is of 
the superheterodyne type, with a center frequency tunable in eight 
steps. The discrete frequencies are 0 , 05 , 0.1, 0 . 2 , 0.35, 0 . 6 , 0.9, 
1.8, and 3 . 5 MHz. 


e Trh 
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The system block diagram is shown in Fig. 2.3. The antenna signal 
is amplified by a low-noise, broadband preamplifier. A low-pass filter 
in the preamplifier rejects undesired frequencies above 3.5 MHz and 
additional notch filters are located at the preamplifier input to 
prevent telemetry signals from entering the radiometer. The prea mp lifier 
output is mixed with signals from a crystal-controlled local oscillator 
that is stepped in frequency to achieve the eight discrete frequencies 
across the band. 

The 6-dB bandwidth of the I.F. amplifier is 10 KHz, and is 
determined by the crystal filter in the mixer output. Frequency stepping 
is in synchronism with the spacecraft data system. The I.F. stage 
has automatic gain control with a useful dynamic range of UU dB. The 
output of the detector is filtered with a time constant of 0.21 seconds. 

An analog voltage is generated in the oscillator circuitry 
which identifies the operating frequency. 

A solid-state four-level noise calibrator is connected in place 
of the antenna for inflight calibration every 9*85 minutes. Approxi- 
mately 37 seconds are required for a calibration, after which the 
antenna is restored to the radiometer input. 

2.1.3 The Timing . The system can operate in two modes. In the 
normal mode of operation the preamplifier output is stepped through 
the eight discrete frequencies already mentioned and is called the 
stepping mode. In the other mode, nonstepping, the local oscillator 
frequency is selected to give an output at only one of the eight 
frequencies . 
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FIG. 2.3 SYSTEM BLOCK DIAGRAM. 
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The sequence of events in the radiometer when it operates in 
the normal mode is illustrated in Fig. 2.4. The train of timing 
pulses furnished by the spacecraft to the receiver at intervals of 
1.152 seconds is shown in Fig. 2.4a. Each pulse causes the radiometer 
frequency to change as shown in Fig. 2.4b. The eight frequency steps, 
f through f , constitute a subcycle of approximately 9.2-seconds 
duration . 

During each frequency step when the spacecraft telemetry rate 
is 1 kb/sec (1000 bits per sec) three data samples are taken, the 
first is the sample containing the frequency identification voltage and 
the last two are the radiometer output data samples. The radiometer 
data samples are taken at timi intervals of 3*31 and 4.37 postdetector 
time constants, measured from the timing pulse. This is illustrated 
in Fig. 2.4c. See also Section 2.1.6. 

The subcycle repeats continuously except when a ground command 
is sent to disable the frequency stepping. 

The calibration cycle consists of four subcycles of eight 
frequency steps each as shown in Fig. 2.4d, similar to 2.4c, but drawn 
to a smaller time scale. The calibrator level, as shown in 2.4 e , is 
held constant during the first subcycle. Thus all eight frequencies are 
calibrated at the same levels. This is repeated for the remaining three 
subcycles at successively increased calibrator power levels. At the 
end of the 36.9-second calibration cycle , the calibrator is 
disconnected and the antenna is reconnected to the radiometer input. A 
calibration cycle occurs once during a main cycle which consists of 64 
subcycles of 9.83 min each. Four of them are devoted to calibration 



TIME 


-29- 



< 



-30- 


as shown in 2.4f. The antenna is connected to the receiver approximately 
9^ Percent of the time. The action of the antenna relay is illustrated 
in 2. kg. 

fLhLllt Preflight Calibration . The preflight calibration of the 
receiver involved several operations but only the final result concerns 
us here, i-e., the measurement of the overall system response over the 
anticipated ranges of operating temperatures and input signal levels. 

Plots of the radiometer output versus noise input signals were 
obtained at each frequency of operation. As is usual in dealing with 
noise power, the input is measured in units of the product RT (ohm K), 
where T is the effective receiver input noise temperature and R is the 
antenna radiation resistance. When the operating temperature was varied 
between -2o C and +^5 °C there was insignificant variation in receiver 
output. We adopted an average curve at each frequency. The corresponding 
values are shown in Table 2.1. As an example, the curve for 3.5 MHz 
is shown in Fig. 2 . 5 . Notice the saturation knee that has important 
implications and that will be discussed later. 

For frequencies from 3*5 MHz to 0.2 MHz the power transfer curves 
are practically the same over the entire range of signal levels. For 
0.1 and 0.05 MHz the curves differ from the others in the small-signal 
region, below 1.5 volts output, because of a higher receiver noise 
temperature. 

In the laboratory calibration the anfenna was replaced by a 
dummy equivalent circuit configured as in Fig. 2.2 using the parameters 
shown in Table 2.2 for all eight frequencies. 




Radiometer Output Voltages Preflight Calibration 


- 32 - 


o' 

0 

0 

00 

ON 

NO 

ON 

1 1 

ON 

O 

CO 

KN 

LTN 

ON 

OJ 

ITN 

n- 

00 

OJ 

O 

ON 

OJ 

LP\ 

NO 

O 

o\ 

d 

d 

OJ 

OJ 

0 

NO 

On 

d 

LT\ 

d 

d 

d 

d 

d 

d 

kn 

KN 

CnI 

OJ 

C\j 

rH 

1 — 1 

1 — 1 

O 


1 — ! 

0 

00 

NO 

1 — 1 

OJ 

CO 

NO 

NO 

ON 

OJ 

OJ 

O 

OJ 


NO 

O 

O 

0 

o\ 

CT\ 

ON 

UJ 

lf\ 

KN 

CO 

KN 

O 

NO 

1 1 

LTN 

OJ 

ON 

CO 


ITN 

d 

d 

d 

d 

d 

d 

KN 

KN 

KN 

OJ 

OJ 

1 — 1 

rH 

d 

d 

OJ 

O 

CO 

NO 

i — I 

OJ 

ON 

NO 

ITN 

O 

O 

0 

0 

N~N 


OJ 


0 

O 

CT\ 

ON 

ON 

CO 

n 

KN 

CO 

d 

O 

NO 

1 — 1 

-d- 

1 — 1 

c— 

ITN 


in 

d 

d 

d 

d 

d 

-d - 

KN 

KN 

KN 

OJ 

OJ 

1 — 1 

1 — 1 

d 

d 


KN 

O 

00 

NO 

1 — 1 

n 

1 — 1 

ON 

1 — 1 

NO 

ON 

O 

LTN 

IN- 

LTN 

OJ 

OJ 

d 1 

O 

ON 

On 

ON 

CO 

NO- 

KN 

ON 

d 

O 

IN- 

1 — 1 

d 

rH 

N- 

n 


ITN 

d 

d 

d 

d 

d 

-d - 

KN 

KN 

KN 

Oj’ 

oj 

1 — 1 

1 — 1 

d 

d 


NO I 

O 

CO 

no 

OJ 

OJ 

0 

irN 

n 

OJ 

Ol 

Ol 


O 

LTN 

ITN 

• I 
O | 

O 

• 

ON 

• 

ON 

ON 

• 

OJ 

• 

NO 

KN 

00 

• 

d 

O 

• 

VO 

• 

1 — 1 d 

1 1 

NO 

d 


n 

d 

d 

d 

d 

-d- 

d 

KN 

KN 

KN 

OJ 

OJ H 

1 — ! 

d 

d 


ON 

d 

0 

0 

CO 

ON 


OJ 

ON 

OJ 

00 

O 

NO 

00 

KN 

CO 

00 

ITN 

d 

o' 

n 

NO 

OJ 
1 — 1 

LTN 

-d- 

OJ 

H 

KN 

IN- 

OJ 

n 

n 

d 

d 

d 

-d 

d 

-d - 

KN 

KN 

KN 

oj 

OJ 

1 — 1 

1 1 

d 

d 


1.8 

0 

0 

CO 

ON 

NO 

ON 

OJ 

On 

co 

0 

NO 

O 

-d 

0 

On 

ir\ 

d 

O 

rH 

IN- 

NO 

NO 
1 1 

0 

n 

LTN 
1 — 1 

00 

IN- 

0 

NO 


ITN 

-d" 

-d 

-d 

-d 

-d 

-d 

KN 

KN 

KN 

OJ 

oj 

1 — 1 

1 — 1 

d 

d 

n 

KN 

0 

0 

CO 

ON 

NO 

On 

r— 1 

ON 

1 — 1 

CO 

n 

n 

n 

KN 

O 

00 

n 

KN 


LT\ 

LT\ 

S' 

O 

d 

10 

n 

N- 

0 

NO 

LTN 

-d 

-d 

-d 

-d 

-d 

d 

KN 

KN 

OJ 

OJ 

OJ 

1 — 1 

1 — 1 

d 

d 


"S' 

O 

H 

O 

1 — 1 

<0 

CM 

H 

0 

1 — 1 

CM 

H 
H 
O 
1 — l 

H 

0 

H 

0 

t — ! 

O 

0 ) 

0 

1 — 1 

0 ) 
0 
1 — 1 


CO 

0 

1 — I 


N 

O 

rH 


co 

0 

1 — 1 


M 

X 

H 

O 

X 

H 

0 

X 

H 

0 

X 

H 

O 

X 

X 

0 ) 

0 

X 

co 

0 

X 

N 

O 

X 

CD 

O 

0 

E-i 

K 

d 

1 — 1 

d 

1 — 1 

d 

1 — 1 

d 

1 — 1 


OJ 

rH 

d 

1 — 1 

d 

rH 

-d- 

rH 


0.45 o. 4 o 0.30 0.27 0.32 0.42 0.65 0.92 


-33- 


Table 2.2 

Different Parameters of the Antenna Equivalent Circuit Shown in Fig. 2.2 

Circuit R (ohm) C^ (pF) 

Laboratory calibration 25 l60 

This study Eq. (2.1) Eq. (2.2) 

2.1.5 Interference , Radiometer Stability and Sensitivity . 

Inspection of the data shows the presence of spurious signals which can 
be classified as impulsive or nonimpulsive interference. Impulsive 
interference is manifested simultaneously in the eight channels as a 
sudden increase in the output level followed by equally sudden decrease 
to the original level. The signals are very strong and last from 
minutes to hours. The four lower frequency channels are usually affected 
by this type of interference. They are not observed in the other 
channels. Since we found no coincidence with spacecraft ground commands 
we concluded that it was due to other experiments onboard. No obvious 
repetitive pattern was found. The fact that this type of interference 
occurs simultaneously in different channels provides a means to discriminate 
it from the Type III solar bursts which show a drift in frequency with 
time. 

The nonimpulsive interference is manifested as permanent high 
levels of the eight outputs, higher than the preflight receiver noise. 

This is especially true in the case of the 1.8 and 0.35 MH Z channels 
as can be seen in Fig. 2.6 where the typical system noise level of the 
radiometer output is shown. The output is measured in terms of the 
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FIG. 2.6 TYPICAL AVERAGE NOISE LEVELS OF THE RADIOMETER. 
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RT product of the antenna radiation resistance and the antenna temperature. 
Figure 2.7 shows that the amplitude of the noise fluctuations in the 
same two channels are much larger than those of the adjacent channels. 

The stability of the radiometer gain was checked by examining 
the inflight calibration approximately every two months between 
March 1968 and December 1989* We found that the output voltage 
corresponding to the four calibration levels stayed constant in all 
channels within a few percent during this period. (Alvarez, 1971) 

To study the sensitivity of the system, we examined the radio- 
meter output when the spacecraft was close to apogee, near l^-7>000 
km, at times in which the baselines of the eight frequency channels 
seemed to be fairly steady and free from solar bursts or impulsive 
interference. Figure 2.6 shows typical values of the baseline 
levels .Their peak-to-peak fluctuations are shown in Fig. 2.7. The 
measurements were made during the first 103 orbital revolutions of the 
spacecraft. We observe that the baseline levels of the system noise 
increase steadily as the frequency decreases, except for the 1.8 and 
0.35 MHz channels that obviously are abnormally high. The arrows 
indicate typical impulsive interference levels that most of the time 
affect the low frequency channels. 

Transformation of RT Product into Flux Density . In Section 2.1.1 
we considered a balanced dipole as a simple approximation of the monopole- 
spacecraft system. By using this model we can estimate the order of 
magnitude of the received flux densities. 

The collecting area A and the gain G of a loaded antenna are 
related by the formula, (Tai, 1961): 
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A 


q i q 2° 


i+7T 


(2.3) 


where q^ is the impedance mismatching factor defined as 

power delivered to the load 

^•1 power delivered to a matched load * 


l+R.Ry 
A L 


(B A + E L )2 + < X A + X L )2 


(2.4) 


also 

Z = R + JX = load impedance;, 

L 

Z^ = R^ + JX = antenna impedance defined at the 

load terminals and 

q^ = polarization mismatch factor. 

If the incident field is randomly polarized, the average value of 


is 


= 0.5 . (2.5) 

In what follows we will assume this condition and in Section 1+.3 
we will discuss further the state of polarization of Type III bursts. 

We will further assume that the antenna terminals are directly connected 
to the receiver input therefore R^ = R and X^ = l/2irfC^. 

The total flux density of randomly polarized radiation incident 
on an antenna is related to the antenna temperature T by the formula: 


kT 

A 


w m 


-2 -1 


Hz 


( 2 . 6 ) 


where k = Boltzmann's constant. 
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Inserting Eqs . (2.1), (2.3) and (2.5) into Eq. (2.6) we get 


k(RT) 

10 q 1 G7rh 2 


(2.7) 


where the product (RT) is the measured quantity. This formula is 
numerically the same for the antenna model consisting of a monopole 
perpendicular to an infinite and perfectly conducting plane. 

Equation (2.7) assumes that the antenna is a balanced dipole but 
the (RT) scale in which our data are expressed were obtained from a 
calibration using a different antenna model. See Table 2.2. Therefore, 
to be consistent, we should define a new calibration curve appropriate 
for the dipole parameters . To compute the correction factor that will take 
this into account let us call the power measured at the receiver 
output. Also let and P.^ be the power collected by the antenna as 
deduced by assuming the laboratory calibration model and by a balanced 
dipole model, respectively. Finally, let qj and q be the corresponding 
mismatching factors. If G R is the receiver power gain and there are 
no ohmic losses then 


and 


Therefore : 
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It is straightforward to compute q’ and q from Eq. ( 2 . 4 ). In 

1 l 

computing Z T ,C^ should be added to C c . The radiometer input shunt 
parameters were measured by Yorks (1968) and are shown in Table 2 . 3 * 

The results of the calculations are shown in Table 2 . 4 . Practically the 
same quantities are obtained for the monopole-infinite plane model. 

We observe that the corrections are large. 

Table 2.3 

Parameters of the Radiometer Input Impedance 
Frequency (MHz) Rg (ohm) Cg (pF) 
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Table 2.4 


Correction Factor for Antenna-Receiver Mismatching 
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We can write Eq. ( 2 . 7 ) as 


S 


f 



( 2 . 8 ) 


for the dipole G = 1 . 5 * The length of the OGO-V boom is h = 9*15 m. 
Then we can write Eq. ( 2 . 8 ) as 


S f = 3-5 x 10 " 27 ( ^) (RT) , w m" 2 Hz -1 . ( 2 . 8 a) 

l 

We will use this equation in Section 4 . 3 . 

The mean sky brightness temperature over the sky is equal to the 
antenna temperature T. The Reyleigh-Jeans approximation, valid in 
our cases, gives the specific intensity 


- 1 + 1 - 


2kT 

k 2 


( 2 - 9 ) 


Using Eq. (2.1) and introducing the correction factor (q ’ /q. ), 

li 


we have 



(2.9a) 


or, for our antenna, 

^.15 x 10 -24 ( i ) (ET) . 

i 

Using (RT) = 9*5 x 10 7 ohm-K from Fig. 2.6 we obtain at 3*5 MHz 
the following specific intensity: 



= 3*0 x 10 19 , w m 2 Hz -1 sr -2 . 

Values obtained by other observers range between 8 x 10~ 21 and 
6 x 10 2 °, R. R. Weber et al. (1969) ; therefore the OGO-V values are 
higher than the measured upper limit by a factor of 5. This means 
that eijther the system noise has contributions other than the cosmic 
noise or that our representation of the antenna by a dipole is incorrect. 
Figure 2.6 shows that the system noise increases with decreasing frequency 
This leads us to assume that the high values of I obtained are due to 
a contribution in excess of the cosmic background. 

2.1.6 Data Acquisition . The data received from the spacecraft 
telemetry system consists not only of the radiometer output but also of 
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other parameters related to the radiometer operation, for example, 
currents, voltages, temperatures, etc. The information is collected 
by the operation of a main commutator (main frame) and two sub-commutators. 
The main commutator has words numbered from 1 to 128. Each consists 
of 9 bits. Eight bits are used for analog data. The lowest repetition 
rate of the main commutator is 1.152 sec which corresponds to 1000 
bits/sec (1 kb/sec.). Of special interest is the word No. 15, that 
corresponds to the frequency identification voltage and Nos. 77 and 
125 which are the two samples of the radiometer output. Since the 
sampling the radiometer output at each of the eight frequencies lasts 
1.152 sec, consecutive samples at any frequency are separated by 9.2 sec. 

At the 1 kb/sec rate this interval constitutes the limitation in time 
resolution of the instrument, not the postdetector time constant of 
0.21 sec. This fact accounts for one of the reasons that our data time 
profiles are not as precise as we would like. 

The repetition rate of the main commutator can be increased by 
a factor of 8 or 64. At the rates of 8 and 6k kb/sec the time between 
samplings, at each frequency, in consecutive main cycles is only O-69 
and 0.08 post-detection time constants. However the frequency-stepping 
rate is always 1.152 sec per channel so the time resolution is not 
improved by using higher bit rates. 

The data from the OGO-V satellites are collected by Goddard Space 
Flight Center (GSFC) ground stations and sent to GSFC where the 
decommutated data tapes are prepared for each experimenter. The data 
are of two types, real time and playback. Real time transmissions 
are at 1, 8 and 64 kb/sec data rates whereas the playback data obtained 
from onboard tape recorders are recorded at 1 kb /sec. 



2.1.7 Telemetry Quantization . The telemetry uses 8 bits in 
a binary system allowing 256 discrete levels in the analog-to-digital 
conversion of an analog voltage between 0 and 5*12 volts. Since the 
output of the radiometer ranges between 0.0 and 5*1 volts, it is 
divided into 256 steps of 0.020 volts each. The design of the system 
is such that an incoming signal that has a value lying between discrete 
steps, is truncated to the adjacent lower discrete voltage level. This 
produces telemetry quantization that combined with radiometer saturation 
results in further limitation in time resolution and in signal amplitude 
resolution. 

The simple diagram of Fig. 2.8 shows the history of a fictitious 
input signal. The strength of the input signal is expressed in (RT) units 
and the receiver output in volts. The telemetry voltage steps V T 
are indicated. Figure 2.8 applies to any of the eight frequencies. 

The input signal (a) drives the receiver (b) to saturation. The dashed 
line follows a particular point (4). The signal is sampled, (c) 
and then the samples are transmitted by telemetry to ground. These 
samples are corrected by a computer using the radiometer response curve 
(d) and are finally plotted (e). it is seen that the incoming signal 
is changed by the telemetry quantization, especially in the nonlinear 
portions of the response curve. The time-of-start of the signal is 
modified and the amplitude is reduced. 

To assess the effect of this problem on the actual observations 
we measured the amplitude of the discrete levels in the range of 
strong signals. The range of interest in the (RT) product was above 
IO 10 ohm-K. We measured many strong signals of different amplitudes so 
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ON LARGE INPUT SIGNALS. 





that we believe that the list of levels shown in Table 2.5 is fairly 
complete. We prepared a diagram of discrete output levels ) for each 
frequency^ whose practical objective was to determine the uncertainty 
in the amplitude of very strong solar bursts. The effect on time 
is negligible in the four upper channels because the system noise in 
flight is in the linear part of the power transfer curve of the 
receiver. The combined effects of saturation and telemetry quantization 
can be appreciated in Fig. 3*2a (Section 3*1) for the event of Nov. 18, 1968. 
Saturation produces the horizontal upper parts of the profiles. The 
staircase-shaped discontinuities in the upper parts are due to the 
discrete voltage steps of the telemetry system. Tick marks on the 
right- and left-hand sides indicate the zero levels. 

2.2 The Spacecraft 

2.2.1 History and Orbit . The OGO-V satellite was launched 
on March 4 , 1968 from the Eastern Test Range at Cape Kennedy, Florida. 

The initial orbital parameters were approximately the following: height 

of perigee 292 km, height of apogee 147,000 km, inclination to the 
equator 51°, and period of 63 hr 25 min. We notice that the orbit is 
highly elliptical and that the height of apogee is O.38 of the Earth- 
Moon distance. The Radio Astronomy Experiment was turned on 19 hours 
after liftoff. The antenna was deployed on the fourth revolution at 
1825 U.T., l4 March 1968. During the time covered in this study 
(March 1968 to February 1970), the radiometer was in the normal stepping 
mode except for the following three periods: 25 April to 18 June 1968, 

12 -l4 September 1969, and 15-17 December 1969* The modes were nonstepping 
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Table 2.5 


Radiometer Discrete Outputs for Large Signals. Units (ohm-K) Should Be 




Multiplied by 10 9 
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at 3 - 5 , 0.6 and 0.6 MHz, respectively. During the last two periods 
the spacecraft was intentionally rotated around the Y-Y axis in a spin 
mode with an average period of 3^ minutes. 

Between March 1968 and February 1970 , the height of perigee 
increased from 6,800 to 21,800 km and the height of apogee decreased 
from 153^000 to 138,000 (geocentric distances). 

Since the spacecraft orbit is approximately fixed with respect 
to the stars it rotates relative to the Sun-Earth direction as the Earth 
moves around the Sun. 

Predictions of the spacecraft position and other related parameters 
are published in advance by GSFC to the experimenters in separate 
magnetic tapes along with the data tapes. The commands sent to the 
spacecraft are published in the Revolution Report Bulletins. Other 
information pertaining to the spacecraft is also published monthly in 
the Operation Summary Reports. The UCLA Fluxgate Magnetometer Group 
has generated 0 G 0 -V Orbital Plots available in microfilm form. 

2 . 2.2 Relative Motion of the Antenna Pattern . Referring to 
Fig. 2 . 1 , the radio astronomy antenna is contained in the plane of 
the solar paddles which are always normal to the Sun-spacecraft line, 
therefore, the antenna is normal to that direction. On the other hand, 
the Z-axis always points to the Earth, therefore the antenna is normal 
also to the Earth-spacecraft line. Thus the antenna is always 
perpendicular to the plane Sun-Earth-spacecraft. 

The plane of the orbit of the spacecraft is inclined about 31 ° 
with respect to the Earth's equator, and it is also inclined with 
respect to the ecliptic. In Fig. 2.9, (P) is the plane of the ecliptic 
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FIG. 2.9 ORIENTATION OF THE OGO-V ORBIT. 



and (R) is the plane of the spacecraft orbit i Plane (q) is normal to 

(P). S and E are the positions of the Sun and Earth, respectively. 

T represents the spacecraft. Thus, the antenna at T is always 

perpendicular to the plane (SET). As the spacecraft moves along its 

orbit the plane (SET) hinges about the intersection SE and it goes 

from a position in the plane (P) when T is at T , to a position normal 

1 

to (P) when T is at T . If we accept that the antenna pattern is that 

of a short dipole, then when the spacecraft is in T the direction of 

1 

the nulls is normal to the plane of the ecliptic (P), and when the 
spacecraft is in T o the direction of the nulls is parallel to plane (P) 
and therefore is perpendicular to (Q). For practical purposes it can 
be considered that the Sun-Earth direction coincides with the Sun- 
spacecraft line because the size of the OGO-V orbit is small compared 
with 1 Astronomical Unit (A.U. ) Then we can say that the direction of 
the nulls rotates in a plane perpendicular to the Sun-Earth line and 
that it coincides with the plane of the ecliptic twice during each 
spacecraft revolution of 65 . 4 hours. The time between successive 
coincidences are not equal because the speed of the spacecraft changes 
along the orbit. Due to the motion of the Earth the spacecraft orbit 
plane (R) rotates about an axis perpendicular to (P) at E completing a 
rotation in about a year. 



CHAPTER III. THE DATA 


3 »1 The Presentation of the Data 

The data can be presented in several forms, the two most useful 
being 35-mm Monitoring film and CalComp plots. The purpose of the 
Monitoring film is to display as much data as possible in a semi-processed 
stage in order to look for special events that can then be fully processed 
in the form of CalComp plots. All of the processing of the data was done 
at the Radio Astronomy Laboratory on the XDS-930 Computer and the associated 
special equipment designed and built by the UM/RAO Laboratory staff. A 
detailed description of the data processing has been given elsewhere by 
Breckenridge (1969) . 

The Monitoring Film . In this format the output voltages of the 
eight frequency channels are displayed simultaneously in 35-mm film. See 
Fig. 3«1* The time axis runs from left to right along the length of the 
film. The outputs of the eight channels are staggered across the width 
of the film. The 3-5-MHz channel is at the top and the seven remaining 
channels are displayed below it in order of decreasing frequency. The 
eight tracings have the same time reference, that is, a line perpendicular 
to the length of the film corresponds to the same time in all of the 
channels . 

The zero voltage levels of the channels are equally spaced across 
the width of the film except for the 1.8-MHz channel. In this case, 
and for unknown reasons, the background noise level is very high; presumably 
because of spacecraft interference. In order to avoid overlap with the 

tracing of the 3-5-MHz output, the 1.8 MHz output level was lowered by 1.0 
volt. 
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FIG. 3.1 SAMPLE OF 35-mm MONITORING FILM. 
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The time axis (in U.T.) runs continuously. During the periods of 
data gaps the computer generates fill data equivalent to a i+O.96 volts 
signal; these are not plotted. Each frame of film contains 39.0 minutes 
of data. It is important to notice that the blank spaces between frames 
are not time gaps but are due to the advancing mechanism of the film. 

The top of each frame contains the following labels: (a) the data and 

time for the start of the frame, (b) the number of the input data tape 
and (c) the number of the file in the input tape. Tic marks along the 
left side across the film indicate the zero voltage level for each of the 
eight outputs; tic marks along the time axis indicate half hours. Because 
of the compressed time scale, the inflight calibration signals are seen 
as vertical short lines about 10 minutes apart. 

If the Attitude Orbit (a/o) tape was available at the time of 
processing the Monitoring Film, some of this information was incorporated 
into the film. This included, among others, time marks for apogee and 
perigee passes. 

The software to handle the data is called the Monitoring Program. 
There are programs to process 1 kb /sec and 8 kb/sec data. The program 
reads the GSFC data tape, unpacks each of the eight outputs, converts 
them into display coordinates, transfers the data to the buffer core and 
then to the photographic cathode ray tube where all eight output voltages 
are plotted simultaneously and are photographed on-line by a 35-mm 
automatic camera. 

Of the two output samples taken during each main commutator cycle, 
only the second one (word 125) is used in 1 kb/sec data. In the 8 kb/sec 
data we use the word 125 of the eighth main commutator cycle. In both 
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cases the sample is taken 5-57 post-detector time constants after the 
frequency is changed. This criterion was adopted in order to arrive 
at Vj 0 of the steady-state value. There is no program for the 64 kb/sec 
data because the effort involved is not compensated by the small amount 
of information recorded at this data rate. 

The CalComp Plots . The periods of time of special interest found 
in the Monitoring Film, or otherwise, are processed in a CalComp Model 
565 plotter. Separate plots are made for each frequency channel. The 
computer corrects the output voltage for the transfer curve corresponding 
to the frequency channel being processed, and then displays it in units 
of RT (ohm - K) product. See Fig. 5*2. Because of the high intensity of 
the solar bursts selected in this study, scales of about IO 10 or 10 11 
(ohm - K) per inch are necessary in order to have the burst time profile 
within the limits of the paper, 8 inches maximum. The time scale used 
was 3 minutes per inch because in so doing, 0.1 inch contain two consecutive 
sample points of the radiometer output and so the samples are readily 
identifiable. This is relevant in measuring the time of start of a burst. 
The pen of the plotter connects the points by straight line segments in 
such a way that the curve intensity versus time is continuous. The CalComp 
plots are conveniently labeled with the pertinent information. The 
calibration cycles are not plotted and their time is filled with a straight 
segment . 

3.2 Data Coverage in Time 

The Monitoring film that was scanned in search of solar bursts 
included data between 14 March 1968 (the day the Radio Astronomy antenna 
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was deployed) and February 1970. For the period from March 1968 through 
June 1969, the data time coverage was fairly complete with a total of 
about 7212 hours. Table 3.1. The data between 24 April 1968 - 18 June 1968, 
during which time the radiometer was in the non-stepping mode, was not 
processed. When it was decided to stop searching for bursts with the 
purpose of analysis, only sparse data had been received for the period 
between July 1969 - February 1970. The data for this time period probably 
did not cover more than 2000 hours. 

3 .3 Characteristics of the Data 

The data may have been marred by any of the following reasons : 

(a) Inflight Calibrations. These last for 36.9 sec and repeat 
every 9.83 minutes and affect all frequency channels simultaneously. 

Sometimes they would occur at the beginning or the middle of an important 
burst. 

(b) Data Gaps in Time Coverage. These had the same effect as (a) 
except that they can be longer but they rarely coincided with important 
bursts . 

(c) Sporadic Spacecraft Interference. This permanently affects 
the 0.35-MHz channel and, most of the time, the rest of the lower ones. 

(d) Ionospheric noise. Strong noise is detected, especially in 
the higher frequency channels, when the spacecraft goes through perigee. 

It may be due to drastic changes in the antenna impedance or to noise 
produced by the ionosphere itself, or induced by the spacecraft. During 
the period scanned for this study the geocentric height of perigee increased 
from 6800 to 21800 km and the geocentric height of apogee decreased from 
153000 to 138000 km. 
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Table 3.1 

Total Amount of OGO-V Data Searched 



Hours 
of Data 

Total Coverage 
in Month ($) 

1968 



March 

609.3* 

95.20* 

April 

371.0 

51.53 

May 

10.3 

1.38 

June 

290.6 

40.36 

July 

711.0 

95-57 

August 

691.7 

92.98 

September 

511.3 

71.02 

October 

552.0 

74.19 

November 

552.0 

76.67 

December 

507.0 

68.15 

1969 



January 

519-0 

69.76 

February 

348.3 

51.83 

March 

476.5 

64.05 

April" 

4 11. 5 

57.15 

May 

4oi.3 

55.95 

June 

460.7 

65.99 

July 

' 364.0 

48.92 

August 

445.6' 

59.90 

September 

255.0 

55-14 

October 

503.0 

67.68 

November 

223.8 

31.08 

December 

55-8 

4.55 

1970 

January 

365.0 

49.1 

February 

336 ** 

50 

* 

This is 
received 

measured from the 5th of the month, 
data. The antenna was deployed on 

the time at which we : 

the l4th. 


** 


Estimated. 
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We may expect some ionospheric frequency cut-off when the spacecraft 
is near perigee. Typical values of electron densities in the ionosphere 
are shown in Table 3-2 (Bowles, 1963; Slysh, 1965). Also shown are the 
corresponding plasma frequencies. The distances are measured from the 
center of the Earth. 

With good approximation, we can compute the effect of the Earth's 
ionosphere, assumed spherically symmetric, on our data. Let us consider 
the revolution No. 117 (31 December 1968 - 2 January 1969) as an example. 
Taking the data of the OR 3/A predictions, mentioned in Section 2.2.1, 
and that of Table 3.2, we found that the spacecraft spent 96.9/0 of the 
period (62.4 hours) outside the 0.1-MHz level and 87.2$ of the period 
outside the 0.05-MHz level. Since the 0.05-MHz bursts occur rarely, we 
concluded that the ionosphere had very little effect on the observing 
capability of the radiometer. In fact, we found only one case of a 
strong event in which the 0.6-MHz burst must have occurred just after the 
spacecraft penetrated the corresponding critical level in the ionosphere 
because there was no indication of it in the 0.6-MHz cha nn el 
while the 0.9-MHz burst was very strong. (3 April 1969 ) 

The amount of time during which the spacecraft was eclipsed by 
the Earth from the Sun is negligible. 

3^4 Selection of the Data 

It became clear at an early stage in the scanning of the monitoring 
film that the quality of the data, as determined by the impulsive inter- 
ference, would impose the following logical division: the data from the 

highest four frequency channels (the high-frequency range) constitute 



-59- 


Typical 

Table 3.2 

Ionospheric Electron 

Densities 

N 

f 

R 

(cm -3 ) 

(MHz) 

(km) 

1.5 x 10 5 

3-5 

7200 

4 x 10 4 

1.8 

8200 

10 4 

0.9 

9800 

4.4 x 10 4 

0.6 

11000 

1.5 x 10 3 

0.35 

13500 

5 x 10 2 

0.2 

14900 

1.2 x 10 2 

0.1 

21400 

31' J 

0.05 

53700 


(*) 


Extrapolated from Slysh (1965, Fig. 2) 
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one group and the data from the lowest four frequencies (the low-frequency 
range) constitute the other. As we saw in paragraph 2 . 1 . 5 , only those 
bursts having a flux density of the order of 10 10 ohm-K (see Pig. 2.6), at 
and below 0.2 MHz, would exceed the limits set by the noise background 
of the system. In the same range of requencies, the requirement on 
minimum flux density rises to about 5 x 10 1 ' L ohm-K when the burst 
occurs during the presence of impulsive interference. Since the four 
low-frequency channels are affected by this type of noise most of the 
time, it turns out that only a few extremely strong solar bursts were 
observed by our radiometer in that range. There is a large number of 
strong bursts that are just below the intensity limit. These bursts are 
usually very well developed in the range of 3 *5-0. 6 MHz. We estimate 
that their number may be about ten times larger than that of the very 
low frequency bursts. 

Since the 0G0-V radiometer has at least two channels below the 
lowest frequency of observation of Type III solar bursts reported in 
the literature, we decided that the most important area of study was that 
of the very low frequencies, down to and below 0-35 MHz. As a direct 
consequence the observed events would be affected, in general, by the 
problem of saturation combined with telemetry quantization. 

It is important to realize that the way in which the data have 
been chosen may introduce a selection effect on the final conclusions 
drawn from this study. A study of the hundreds of bursts in the range 
between 3.5 and 0.6 MHz would shed light on this matter. The amount of 
data scanned is approximately that shown in Table 3.1 and it covers 
about 9000 hours. In these data we found 64 events that satisfied the 
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selection criterion. Inspection of data processed after the 64 bursts 
were selected indicates that there were 15 more bursts in the period 
extending from July 1969 through middle February 1970. These 15 bursts 
were not included in this analysis. The 64 selected events that we could 
properly call kilometric ~ wave bursts are listed in Appendix 1 , Table A.l. 



CHAPTER IV. ANALYSIS OF THE SELECTED EVENTS 


■ 1 Measurable Parameters 

In the last section of the previous chapter we stated the criterion 
adopted to select the radio events for further study. They should be 
observed down to or below 0.35 MHz. Since the noise system level of the 
four lowest channels is usually high because of radio interference, a 
burst to be observed in these channels needed to be extremely strong. 

We found that a characteristic inherent to all these events was a very 
complex time structure. We have implied that a very strong event is also 
very complex. 

For discussion it is appropriate to establish some definitions. 

We will call EVENT the complete radio phenomenon from the time it starts 
in our records until it ends. This definition may be used in a more 
restrictive sense when it refers to only one frequency. The objective 
of the definition is to distinguish the whole complex time profile, EVENT, 
from its individual COMPONENTS, or BURSTS components. Because of 
limitations in time resolution due to sampling rate or to saturation 
effects already discussed, the number of identifiable components 
constitutes a lower limit of the actual number. 

The profiles of intensity versus time that the radiometer provides 
at each frequency allows us to measure the following parameters of interest 

(a) Time of start of an event and of its components at the different 
frequencies, 

(b) Rise time, defined as the time it takes a burst to rise within 
certain fractions of its peak value. 
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(c) Decay time, defined as the time it takes a burst to decrease 
from its peak to 37 percent of that value (e 1 ), 

(d) Duration of an event or its components, 

(e) Peak intensity. 

4.2 Time of Start 

As it has been explained in Section 2.1.6, the maximum accuracy 
with which we can measure a starting time is 9*2 sec. We can reach 
this accuracy in the 3*5, 0.9, and 0.6 MHz channels. Since the 1.8-MHz 
channel is usually noisy the accuracy is reduced. The accuracy is 
always poor in the 0.35-MHz channel and it varies between good and poor 
in the 0.2-, 0.1- and 0.05-MHz channels, depending on the presence of 
impulsive interference. This is usually present at the times of 
occurrence of the selected bursts. However, this does not render our 
data useless because the delay in the times of arrival of the bursts 
increases rapidly with decreasing frequency and the accuracy in the 
measurement becomes less critical. 

The times of start of the components of an event were plotted in 
graphs displaying frequency versus time. The scale used was 2 inches/MHz 
in the ordinate and 3 min/inch in the abscissa. We plotted not only the 
time of start of well identifiable components but also those times at 
which a discontinuity in either the rising or decaying part of a burst 
profile suggested a buried component. Very often in the case of strong 
events and because of the quantization-saturation effect, the peak 
becomes a horizontal trace of up to several minutes of duration. In 
those cases it is impossible to guess how many components are hidden. 
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The uncertainties in the times of arrival were represented in 
the plots by horizontal segments. The inaccuracies in frequency are 
not noticeable in the chosen scale. The time of start of the 3-5-MHz 
leading component was used as reference. An example is shown in Fig. 4.1. 

Sometimes a few or all points corresponding to a frequency channel 
may be missing because of gaps in the data, or interference, or just 
because the burst was weak. 

k.2.1 Observed Frequency Drift Rates . A direct result from the 
measurements of the arrival time of a burst is its frequency drift rate (FDR). 
The first step in the technique used to determine the FDR consisted in 
drawing a smooth curve through the observed points corresponding to a 
certain component in the plot of frequency versus time of arrival. The 
second step was to measure the slope of the curve thus obtained at a 
certain number of points. The advantage of the method is that it handles 
graphically the uncertainties in the time of arrival. It is equivalent 
to a weighted graphical averaging. Whenever it was possible two extreme 
smooth curves compatible with the data would be drawn. The slope measured 
at a certain frequency in these two curves would give the upper and lower 
limits of the FDR at that frequency. Because we had obtained a continuous 
curve in frequency versus time we could measure the slope at any frequency. 

In each case we selected the two extreme frequencies in the following way: 
the highest frequency was midway between the first and the second present, 
the lowest frequency was midway between the last one present and the one 
before the last. The intermediate points were usually selected at the 
frequencies of observation except in special cases where more points 
were needed. The values ; of FDR so obtained were plotted on log-log paper 
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FIG. 4.1 EXAMPLE OF FREQUENCY VERSUS OBSERVED ARRIVAL TIMES. 



- 66 - 


with df/dt as ordinate and f as abscissa. The points were then connected 
by straight segments with an approximate best fit. 

The validity of the method depends on the unmistakable identi- 
fication of one and the same component at each of the frequencies of 
observation. This limited the applicability of the technique because 
of the extreme complexity of the selected events, plus the lack of time 
resolution. We thought of two cases in which the method would be 
applicable. In the first case the criterion adopted was to use the time 
of start of the event and assume that at each frequency this corresponds 
to the time of start of the first component of the complex. In the 
second case, whenever an event contained a component several times 
stronger than the rest it was assumed that this component could be 
followed without too much ambiguity at lower frequencies. As we will 
see shortly the assumption involved in the two cases were not always 
found to be valid. Those events in which the assumptions were correct 
permitted us to extend the knowledge of FDR of Type III bursts down to 
0.075 MHz. The others led to the discovery of the unexpected prevalence 
of second harmonic emission below about 1 MHz. (Section 4.2.5). 

By using the above technique we obtained FDR curves for the 
first 52 selected events. (Of the rest, 11 were selected from the raw 
data after we found problems in interpreting the results and therefore 
we did not process them in this way, and one was useless for this purpose. ) 
We found that, based on their shapes, we could classify the curves into 
several groups that we designated as A, B, C, D, E, and Unclassified. See 
F'ig. 4.2. Haddock and Alvarez, (1970). The FDR's of the bursts in the 
last group were not amenable to the drawing of any meaningful curve. It 
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turns out, as we will see later, that the only cases for which the 
curves really represent FDR are those in group A. The peculiar shapes 
of the other groups are artifacts produced by a technique that does 
not discriminate between fundamental and second harmonic emission. 

Fortunately the first bursts selected for analysis were those 
in group A. They were chosen because of the simplicity of their FDR 
curves: a straight line. The technique to be described in Section 4.2.5 
revealed that the emission of these bursts at the different frequencies 
was mostly at the fundamental or at the second harmonic, to the 
extent that this could be determined. The FDR curves were obtained for 
17 of the 53 bursts. (The number is 53 and not 52 because in event 
No. 62 we could measure two components.) The numbers of the corresponding 
events, as given in Table A.l of Appendix 1, are: 6, 13, 16, 18, 19, 

20, 22, 24, 25, 27, 29, 30, 31, 34, 62. 

A straight line can be characterized by its slope, a, and by 
its intercept with the ordinate axis, a. We can thus write the equation 
of a straight line in a log-log plot as : 

d f* 

log — jj- = a + a log f . (4.1) 

In our plots we expressed f in MHz and t in seconds. The 
parameter "a" is the logarithm of the FDR at 1 MHz. The slope of the 
line is positive, a > 0. 

Equation (4.1) can be written: 


df 

dt 



(4.2) 



-69- 


where the negative sign has been added, because as the time increases 
the frequency drifts from high to low values. The form of this equation 
suggests defining a as the FDR Index . Within the accuracy of our 
measurement this relationship is valid at least between 2.65 MHz, the 
mid frequency between the two highest channels, and 0.075 MHz, the 
mid frequency between the two lowest channels. We will see very 'soon 
that its validity holds up to 550 MHz. (Haddock and Alvarez, 1969)- 

The measured values of a and a are given in Table 4.1 and their 
frequency distributions are presented in Fig. 4.3. 

Figure 4.3a shows that the spread in the distribution of "a" is 
not too large while Fig. 4.3b shows that the spread in a is larger. 

It is interesting to observe that a does not become less than 1.0. 

We will refer to this point in Sections 4.2.3 and 4.2.4. The distri- 
butions have the following values of average and standard deviation: 

a = I.92 , a = -1.93 , a = 0.38 , a = 0.13 • (4.3) 

LX a. 

Table 4.1 does not contain observations below 0.1 MHz. In fact, 
the only event observed down to 0.05 MHz is the event No. 21. The 
characteristics of its FDR curve are those of a straight line except 
close to the higher frequencies where the points lie somewhat above the 
extrapolated straight line. Because of this closeness and because of the 
uniqueness of this event it was included in our calculations but it 
was treated in a different approach. This consisted in averaging, at 
each of the frequencies selected for FDR measurements, the values of 
FDR obtained for the 18 different bursts.; These averages 
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Table 4.1 

Parameters of Frequency Drift Rate Curves for Bursts of Group A(*) 


Event No. 

a 

a 

6 

- 1.92 

1.42 

13 

- 1-75 

2.05 

16 

- 1.90 

1.98 

18 

- 1.90 

1-31 

19 

- 2.18 

2.24 

20 

- 1.96 

2.01 

22 

- i.84<**> 

2.27 

2 4 

- 1.83 

1-97 

25 

- 2.02 

1.83 

27 

- 2.01 

1.59 

29 

- 1.83 

2.48 

30 

- 2.o4 

2.03 

31 

- 1.90 

2.3^ 

34 

- 1.65 

1.21 

42 

- 2.03 

1.93 

62 

- 2.03 

1.60 

62 

- 2.08 

2.32 


(*) 


These measurements were made from graphs in which f was 
in MHz and time in seconds. 


expressed 


(**) 


Averages. 
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were plotted against frequency in a log-log graph and a straight line 
fitted by the least-squares method. The frequencies used and the 
averaged FDR obtained are presented in Table 4.2. The least-square 
fitting of a straight line gave the following parameters: 

a = 1.93 , a = -1.89 f a = 0.08 , (4.4) 

where a is the standard deviation, of the dispersion in ordinates of 
the data points with respect to the fitted line. The line and the 
data of Table 4.2 are shown in Fig. 4.4. 

We see that the values of a and a found by the two methods are 
not significantly different. For describing an average FDR curve we 
have chosen the second procedure because it is the one consistent with 
the form in which the data at other frequencies, that we will soon use, 
are given in the literature. 

The next step in the study of FDR of Type III solar bursts was 
to find the extent of the validity of the empirical law described by 
Eq. (4.2). For this we collected some of the data on FDR published in 
the literature. These data are shown in Fig. 4.5 that also shows the 
OGO-V points listed in Table 4.2. The least squares fit of a straight 
line through the points other than OGO-V 's gave the following parameters: 

a = 1.88 , a = - 2.07 , a = 0.13 . (4.5) 

The least square fit to a straight line of all the points 
contained in Fig. 4.5 except Young £t al_. , gave the following results: 
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Table 4.2 


Average Frequency Drift Rate for 18 Bursts 




df 


df 


-1 

fiTT 17 non 

f (MHz) 

log 

dt 

a ( } 

dt 

r 

in 4 bet 

2.70 

- 

1.07 

0.16 

8.5 

X 

lO" 2 

1.80 

- 

1.44 

0.10 

3.7 

X 

10" 2 

o'. 90 

- 

1.90 

0.19 

1.3 

X 

10“ 2 

0.75 

- 

1.98 

0.23 

1.0 

X 

10“ 2 

0.60 

- 

2.36 

0.17 

4.4 

X 

10" 3 

0.48 

- 

2.48 

0.22 

3*3 

X 

10" 3 

0.35 

- 

2.86 

0.11 

1.4 

X 

10" 3 

O.275 

- 

3.08 

0.15 

8.4 

X 

io“ 4 

0.20 

- 

3.28 

0.12 

5-3 

X 

10“ 4 

0.15 

- 

3.47 

0.l4 

3-4 

X 

10" 4 

0.10 

- 

3-79 

0.13 

1.6 

X 

10- 4 (**> 

0.075 

- 

4.01 

0.15 

9-8 

X 

10 -5 (.*) 


' ' This is the standard deviation corresponding to the values in 
the second column. 

' ' These two values come from only one observation. The error 

quoted corresponds to the maximum and minimum values compatible 
with the observation. 


FREQUENCY DRIFT RATE (MHz sec') 
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FIG. 4.4 AVERAGE FREQUENCY DRIFT RATES FOR 18 BURSTS OBSERVED 


10 


BY OGO-V. 


OGO X 

Haddock 8s Graedel (1969) 
Hughes 8 Harkness (1963) 
Dulk 8 Clark (Hartz,l969) 
Young et al. (1961) 
Alexander et al. (1969) 
Maxwell et al. (I960) 
Boischot (1967) 

Hartz (1969) 


FREQUENCY (MHz) 


FIG. 4.5 FREQUENCY DRIFT RATES FOR TYPE III BURSTS. 
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a = 1.84 , a = -2.00 , a = 0.13 • (4.6) 

A line with the parameters of Eq. (4.6) is shown in the same figure. 
The summary of all these data is presented in Table 4.3. 

Table 4.3 

j -P 

Least Square Fit of the Line log — = a + a log f 

& u 


Data 

a 

a 

a 

OGO-V. Table 4.2 

1-93 

- 1.89 

0.08 

Other than OGO-V' s 

1.88 

- 2.06 

0.13 

All 

1.84 

- 2.00 

0.13 


In order to represent the variation of the frequency drift rate 

with frequency we will adopt the parameters derived from "all" the 

data. The error that we will make by using this line to represent the 

OGO-V data is several times less than a and a . We can conservatively 

a a J 

state that Eq. (4.2) is valid between 550 and 0.075 MHz and can be 
written 

Tj-j- = - 0.01 f MHz sec 1 

The data of Young et al. (1961) between 500 and 94o MHz were not used 
because it is not certain that the bursts are of the classical Type III, 
(Kundu et al. 1961). However they agree with the values expected from 
an extrapolation of the computed line. Young et al. found that the frequency 
drxlt rute is = 2000 MJIz sec . 
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The line that fits all the observations in Fig. ^.5 deserves 
special attention because it is remarkable that such a good linear fit 
can be obtained over seven decades in frequency drift rate with 
observations made in such a variety of circumstances. For example, the 
earliest observations included are from Maxwell et_ al_. (i960) and were 
taken close to the sunspot maximum of 1958.25* On the other hand, the 
OGO-V data was obtained at the next sunspot maximum, in 1968. The data 
contained in Fig. ^.5 span 80 percent of an 11-year cycle, from peak to 
peak. The positions of the associated flares are not known, except for 
OGO-V data, and we can expect them to cover a wide range of values. The 
velocities of the exciter particles also are known to have an appreciable 
distribution of values. Finally, the coronal electron densities as well 
as the emission mode (fundamental or second harmonic) associated with 
each measurement may have been different. In spite of all this, the 
standard deviation in the spread of the points in log |dt/dt| from the 
best fit line is only 0.13 . 

Attempts have been made in the past to find a mathematical 
expression for the FDR as a function of frequency. Wild (1950) found 
in the range 70 to 130 MHz : 

df f - 1 

^ , MHz sec . 

In fact, he started with the same relationship as Eq. (^-.2) and he then 
used oc = 1. 

There was early recognition that the plot of log|df/dtj versus 
log f was a straight line (for example, Hartz 1969)* Graedel (1969) 
wrote an expression in the form of Eq. (U.l) that he used to' compute 
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electron density gradients in the solar corona. The form Eq. (4.2) 
simplifies mathematical manipulation and retains the sign. 

We can summarize this section as follows: 

(a) We have developed a technique for measuring the FDR of Type 
III bursts. Its indiscriminate use in cases where it is not applicable 
led to the disclosure of an unexpected prevalence of second harmonic 
components, as will be explained in Section 4.2.5. 

(b) We established that a plot of frequency drift rate versus 
frequency in a log-log scale can be well represented by a straight 
line down to 0.075 MHz. 

— Basic Assump tions . The analysis of the frequency drift 
rates of Type III bursts depends on a knowledge about the following 
subjects : 

(a) Electron density of the solar corona, 

(b) Mechanisms of emission of Type III bursts, 

(c) Trajectory of the exciter particles, 

(d) Velocity of the exciter particles, 

(e) Propagation of the electromagnetic radiation through the 
solar corona to the Earth’s environs, 

(f) Other subjects. 

The purpose of this section is to obtain models of the electron 
density distribution that are consistent primarily with observation and, 
if possible, with theory. These models will be based on assumptions 
regarding the last five subjects. We will start by discussing these 


subjects. 
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(a) Electron Density Distribution of the Solar Corona 

The only ranges of distance studied experimentally before the 
OGO-V observations were near the Sun and near the Earth. 

Optical observations have provided information to distances 
out to 20 solar radii (R @ ). Some of the optical studies were done by 
Baumbach (1937), van de Hulst (1950), Michard (195*0, Blackwell (1956), 
Blackwell and Petford (1966), Gillet _et al. (196^), Newkirk et al. (1969), 
Hansen et al. (1969), Newkirk et al. (1970) and Schmidt (1953)* 

The electron density N of the solar wind near the Earth has 
been determined mainly by experiments onboard spacecraft. Montgomery 
et al. (1968), Koehler (1968), Neugebauer and Snyder (1966a, 1966b), 

Wolfe et al. (1967), Bonetti et al. (1969), Stelzried et al. (1970), 
and by the group at Stanford University headed by V. R. Eshleman (see 
Staff, 1966). 

A determination of N at 1 A.U. using observations of the zodiacal 
light was done by Blackwell and Ingham (1961), by Ney et al . (1961) 
and also by Elsasser (1963). Radar techniques have also been used 
involving reflections from Venus (Campbell and Muhleman, 1969, Brandt, 
196^; Smith et al. 1968), from the Moon (Howard et al. 1965) and from 
several spacecrafts (Howard and Koehler, 1967; Howard, 1968). 

Measurements near Venus were made by Mariner V and were reported 
by Bridge et al. (1967). It is not clear however if the measurement 
actually represents the solar wind electron density because of the 
nearness of the spacecraft to the Venus shock front. 

Other radio techniques to study the solar corona are based on 
the observation of the scintillation of radio sources, on the occulation 
of radio sources by the solar corona and on the observation of pulsars. 
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Care has to be exercised however in selecting electron density- 
values near the Earth because many of them are inferred through the use 
o-f-" particular model. This is the case of the results, for example, 
of Stelzried et al. (1970), Blackwell and Ingham (1961), Campbell and 
Muhleman (1969), Brandt (1964) and Smith et al. (1968). 

Very near the sun N varies typically between 10 s and 10 9 cm" 3 for 
the corona and up to more than 10 9 cm 3 for coronal streamers. 

Near the Earth there is also a wide range of values that go 
from 1 cm 3 to 100 cm 3 , or more. However the average is generally 
accepted as lying between 5 and 30 cm" 3 . 

We have recognized the fact that the electron densities depend 
on the phase of the solar cycle in which they were measured, however we 
did not consider this effect in our models. 

(b) Mechanism of Emission 

There is little doubt that the radio bursts originate in plasma 
oscillations set up by the passage of exciter particles through the solar 
corona. These oscillations represent longitudinal waves which can couple 
into transverse electromagnetic waves by mechanisms which are not well 
understood. This is the so-called Plasma Hypothesis. As was mentioned 
in Section 1.1.2 the original idea was suggested independently by 
Shklovskii (194-6) and by Payne-Scott et al. (1947) but it was developed 
primarily by Wild (1950). 

Our observations and others have revealed the existence of the 
fundamental emission as well as the second harmonic, at twice the fundamental 
frequency. The identification of harmonics in Type III bursts was made 
at least in 1954 by Wild Murray, and Rowe (1954). 
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The well-known plasma frequency equation is : 

f = j 9 vf , (1+.7) 

where f = frequency, Hz, 

N = electron density, m 3 , 
j =1, for the fundamental, 
j = 2, for the second harmonic. 

The solar corona has a negative electron density gradient. Any 

of our single-frequency channels, let us say centered at frequency f , 

can receive indiscriminately hurst radiation emitted as fundamental 

from a coronal level whose plasma frequency is f , or/and emitted as 

o 

second harmonic from an outer level whose plasma frequency is 1/2 f . 

o 

It will he important then to distinguish between fundamental emission 

at a frequency f , its second harmonic at frequency 2f , and the emission 
o o 

at frequency f Q generated as second harmonic. 

In this study, we have assumed that a fundamental frequency f and its 
second harmonic at 2f are produced simultaneously and in the same general 
region. 

(c) Trajectory of Exciter Particles 

In the past it has been sufficient for observers to consider that 
the exciter particles were ejected along a straight path. The distances 
from the Sun to the plasma levels corresponding to the lowest frequencies 
of observation.', then available were not large enough to require a 
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discrimination between the shapes of the trajectories. We realized that 

the source of our 0.05 MHz burst could be placed at 100 R A , or beyond. 

u 

This made it necessary to choose a path for the particles. 

The two obvious choices were a straight line coming radially out 
from the center of the Sun and an Archimedes spiral as suggested by the 
solar wind studies. In the first case the particles would not be 
influenced by the solar magnetic field, except perhaps very close to the 
Sun, therefore the trajectory would be presumably along the radius joining 
the center of the Sun with the center of the associated optical flare. 

This would lead to a simple three dimensional geometry. In the second case, 
the exciter particles would be guided by the interplanetary magnetic field 
along an Archimedes spiral. 

We expected that the observations would permit us to decide which 
trajectory was a better assumption. 

The Spiral Trajectory . Parker (1958) based on studies made by 
Alfven proposed that the interplanetary magnetic field would be frozen-in 
to the solar plasma and dragged by this plasma as it streams continuously 
out from the Sun. The magnetic field lines would be shaped into 
Archimedes spirals whose characteristics would depend on'the rhte 6f 
rotation of the Sun and on the velocity of flow of the solar plasma. 
Experimental studies made first by McCracken (1962) proved Parker's 
prediction. Observations performed later have confirmed McCracken's 
findings. (Dessler, 1967; Kavanagh et al. , 1970). Because the plasma 
outflow is mainly in the equatorial plane of the Sun the magnetic field 
may be expected to be parallel and close to it except perhaps near. the 
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Sun, (Gold, 1958; Jager, i 960 , Sturrock and Smith ,‘ 1968 ) . For simplicity 
we assumed the spiral trajectories to he contained in the plane of the 
ecliptic in spite of the fact that it is actually inclined about 7*5 
degrees with respect to the solar equator. In general this will result 
in an underestimation of the exciter velocity because the introduction 
of the third dimension makes the distances larger. Figure h.6 shows the 
geometry involved. The plane of the ecliptic coincides with the plane 
of the paper. S represents the Sun and E the Earth. 

0^ * heliocentric longitude of the associated flare and 

& = one astronomical unit (A.U.) 

The equation of the spiral shown in Fig. k.6 is 

0 = 0-. -e - , (fc.8) 

F a 

where e is the tightness of the spiral which can be shown to be: 

e = tan </> = SL , (1+.9) 

E V s 

where ^ is the angle at which the spiral intersects the radius vector 

Jcj 

at 1 A.U. Also, e is directly proportional to the Sun's angular rotation 
and inversely proportional to the solar wind velocity V . At the 

O 

solar equator Q = 2.7 x 10 6 rad sec 1 . Typical values of V at the 

D 

Earth range between 300 and U00 km sec 1 . (Dessler, 1967 ; Parker, 1969 ). 
In this study we assumed V to be constant and independent of distance. 

D 

Some values of ^ measured and predicted are presented in Table U.ii. 

After approximating 9 = 0^, at r = 1, we adopted in this study e = tt/ 2 or 
= 57-5 degrees. Thus, exciter particles originating on the west limb 
will pass through the Earth. 
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4.6 GEOMETRY OF THE TRAJECTORY OF THE EXCITER PARTICLES. 



-85- 


Table 4.4 

Characteristics of the Solar Wind Spiral Structure 


4>° 

(D 
’t) O 

Author 

45 

57 

Parker (1964) 

~55 

.90 

McCracken (1962). 

30-40 

~50 

Fan et_ al. (1968) 

55 


Wolff et al. (1971) 

46 

60 

Wilcox (1968) 

57-5 

90 

This Study 


The Radial Trajectory . The three-dimensional form of this geometry 
was not used. Instead we made some use of radial paths contained in 
the plane of the ecliptic. These results were obtained directly from 
Eq. (4.8), making e = 0. 

(d) The Velocity of the Exciter Particles 

Following the local plasma hypothesis, we have accepted that the 
bursts are produced by a group of very fast particles traversing the 
solar corona with certain velocity v. The original hypothesis assumed 
that the velocity was constant during the motion of the particles. This 
matter has not been settled. Stewart (1965) has accepted the constancy 
of v. Malville (1972a) observed a burst whose velocity seems to decrease 
from 0.45c in the 180-50 MHz range to 0.l8c in the 19-9 MHz range, where 
c is the velocity of light in vacuum. However we believe that his 
evidence is inconclusive because in the light of our findings, discussed 



- 86 - 


in Section 4.2.5, the same velocity pattern could be given by a burst 
with a velocity 0.40c radiating at the fundamental frequency in the 
high-frequency range but then radiating, at the same frequency but as 
second harmonic, in the low-frequency range. The "switch" in mode 
prevalence from fundamental to second harmonic produced by a constant 
exciter velocity could be observed as an apparent reduction in velocity 
by a factor slightly less than two. 

Boischot (1967b) working in the range from 40 to 8 MHz found that 
the velocities derived from the bursts radiating in the fundamental mode 
decreases with increasing distance from the Sun and he also found the 
same trend for the velocities derived from the second harmonic. 

Close to the Sun we could think that some ion-electron collisional 
effects and magnetic effects are strong enough to produce deceleration, 
for example, by energy losses, however at distance larger than a few 
tens of solar radii we can suggest no obvious cause to reduce the velocity 
of the particles. Small-scale inhomogeneities in the electron density 
and small-scale structure in the interplanetary magnetic field may 
certainly cause fluctuations in the velocity, however the timing accuracy 
of our observations may be unable to resolve them. 

Considering the above discussion we decided to use a constant 
velocity. Referring to Fig. 4.6 we have: 

arc PP^ = ds 

and 


v 


ds 

dt 


(4.io) 
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We will also use the velocity parameters 8 defined by 

6 = - , (U.ll) 

c 

where c is the velocity of light in vacuum. The observations will give 
information about the velocity of the particles associated with each 
burst . 

(e) Propagation of the Electromagnetic Radiation 

In a plasma the trajectory of a ray traveling at an angle to a 

direction in which there is a density gradient is not a straight line. 

\ 

In connection with the solar corona this problem has been treated by 
several authors (Smerd and Westfold, 19^9; Jaeger and Westfold, 1950; 

Link, 1951; Reule, 1952; Bracewell and Preston, 1956; J. C. James, 1966). 

The results achieved by these authors are valid only in the high frequency 
range where the distances of the plasma levels to the Sun are small. 

For example, Jaeger and Westfold, and also Bracewell and Preston used 
Baumbach's model of the solar corona modified by Allen; Link used an 
inverse power law distribution and Reule used the original Baumbach model. 

. Since they considered rays originating (or reflecting) at points 
close to the Sun, the work of these authors was simplified because they 
could assume that all the rays observed at the Earth came parallel to the 
Sun-Earth direction. (See Fig. 1^.7a). But for a burst at a very low 
frequency originating at a distance that is an appreciable fraction of 
1 A.U. the situation is different, as we can see in Fig. L.7b. Accordingly, 
the approach has to be different also. 
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FIG. 4.7 RAY PATHS IN THE CORONA FOR HIGH AND FOR LOW FREQUENCIES. 
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At the time we did this work our expression for the electron 
density distribution was complicated (see Eq. A.l^) and it was not 
appropriate for analytical manipulation that could lead to a simple 
answer. We considered it premature to use more involved numerical 
methods . 

Because of the complications mentioned above we decided to 
assume in this study that the electromagnetic energy propagates along 
a straight line to the Earth, travelling with the velocity of light in 
vacuum. 

There are several implications involved in this assumption. First, 
if we assume a spherically symmetric corona with a negative electron 
density gradient then, strictly speaking, we could not see the bursts 
ocurring within the surface of escape of radiation defined for the 
frequencies of observation and for the Sun-Earth direction. (See Fig. 
l».7a). At lower frequencies this restriction is more limiting because 
the lower the frequency the larger the proportion of the ray path that 
will be in a region where the index of refraction is near zero, and 
this is the region of appreciable refractive effects. The geometry of 
the trajectory of the exciter particles also has a bearing on this 
problem, as it can be realized by thinking of particles ejected from a 
west limb flare along either a straight line or a spiral. The nature 
of the radiation, whether fundamental or second harmonic, is also related 
to this problem because a frequency emitted as second harmonic originates 
at a level where its index of refraction is * 1. (Wild, Sheridan and 
Neylan, 1959)- 
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The problem of escape of the radiation has received attention by 
Shain and Higgins (1959) > who proposed a solution based on the existence 
of small-scale irregularities in the solar atmosphere from where the 
radiation would be scattered to the Earth. 

The alternative to a spherically symmetric corona is a corona 
containing large scale enhancements like long streamers. The concept of 
solar streamers became necessary in radio work because since the early 
studies, the models of electron densities that were required to fit the 
observations were found to be about 10 times as dense as the typical 
coronal models (Shain and Higgins, 1959; Wild, Sheridan and Neylan, 

1959). The critical level of a certain plasma frequency is higher in 
the streamer than in the homogeneous corona. We will come back to this 
subject when we discuss the results. 

Second, the assumption that the electromagnetic waves travel in 
a straight path implies that the velocities obtained for the exciter 
particles are underestimated. 

Third, the assumption of propagation as in vacuum conditions ignores 
group retardation effects which will also result in an underestimation 
of the velocities (Wild, Sheridan and Neylan, 1959). 

The problem of the effect of a plasma on the propagation of 
electromagnetic waves is closely related to the mode of emission. Stewart 
(1962) found that the second harmonic precedes its fundamental and he 
attributed it to group retardation effects. Smerd et al . (1962) observed 
that m some cases the second harmonic arrives from a direction closer 
to the center of the disk than that of the fundamental at the same 
frequency. Elgaroy and Rodberg (1963) have also dealt with the problem 
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of the delay between fundamentals and second harmonics at the same 
frequency. The theory of Type III bursts developed by Smith (1969) 
and based on these observations succeeds in predicting the early arrival 
of the second harmonic with respect to its fundamental. 

With the simple equation discussed later for the electron density 
distribution it should be easier to refine our calculations by introducing 
some of the effects discussed in this paragraph. 

(f ) Other Factors 

In the derivation of the electron density distribution no 
temperature effects were considered; no magnetic field influence was 
considered either except for the guiding of the exciter particles along 
a spiral. 

4.2.3 Electron Density Distribution — Form 1_. In our study of 
the distribution of electron density in the corona we arrived at two 
mathematical expressions that we will call Forms 1 and 2. In spite of 
the fact that we consider only Form 2 as being important we will discuss 
briefly Form 1 because it led to Form 2 and also because some results 
presented later were based on it. 

The Eq. (4.2) represents a differential equation that together with 
the relationship between the plasma frequency and electron density suggested 
that we could go one step further beyond obtaining electron density 
gradients, and by integration solve for the electron density distribution 
itself. This was done using Eqs . 4.2, 4.7, 4.8 and 4.10 plus the 
assumptions discussed in the previous section. The details are given in 
the Appendix No. 2. The final result is Eq. (A.l4) that represents Form 1. 
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Conceptually there are objections to this expression because the electron 
density distribution depends on the velocity parameter 8 and on the 
position of the associated optical flare. In order to determine the 
constant of integration we used the value of the electron density at an 
arbitrary distance. Equation (A.l4) permits us to impose arbitrary 
conditions at the base of the corona but not at the orbit of the Earth 
because the function may become ill-behaved at shorter distances. 

On different grounds we may expect the Eq. (A. 1*0 to be invalid 
for large distances, say r ^ 1 A.U. We recall that by hypothesis 
particles travelling through a density distribution given by Eq. (A. 14), 
under the assumptions made, give a frequency drift rate curve expressed by 
Eq. (4.1) , that is, a straight line in a log-log plot that extends from 
an infinite frequency down to zero frequency. This is not physically 
reasonable . 

In spite of the above deficiencies we found that, using some 
reasonable values for the parameters involved, Eq. (A.l4) gives 
quantitative results that are consistent with the observations (electron 
densities and frequency drift rates) within the orbit of the Earth. 

Figures 4.15 and 4.17 in Section 4.2.5 show examples of models based on Form 1. 

Electron Density Distribution — Form 2. These deficiencies 
of the function, Eq. (A.l4), led to search for another one that would 
eliminate the unwanted features while keeping the desirable ones. This 
function should give results in agreement with the observations of 
frequency drift rates and with observations of electron densities. If 
possible it should also be simpler. The form of the new expression, 
suggested by Eq. (A.l4) itself, is: 
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N = 


(r - b y 


(4.12) 


We will refer to this as electron density distribution Form 2. 

In order to arrive at acceptable models based on this expression 
we need to adopt values for the three parameters and to test the models 
against observations. 

We will express r and b in units of solar radii and A in cm 
The radial distance r will be measured from the center of the Sun. 

We will adopt the exponent p found in Eq. (A.l4) and whose 
expression is given by Eq. (A. 11), reproduced below: 


p - 5-rr • (a. ii) 

As before a is the frequency drift rate index. Using the value of 
a = 1. 81+ obtained from data over all the range of frequencies as shown in 
Table 4.3 we arrive at 


P = 2.38 . (4.13) 

It is interesting to compare the values of p derived from the 
radio data in this manner with those obtained either from observations 
or theories. Table 4.5 presents some values at distances at which the 
curve log N versus log r is very close to a straight line and therefore 
p could be measured easily. 

We notice that the radio determinations agree well with the 
observations and with the theory of either solar wind or streamers models. 
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Table 4.5 

Values of p Measured from the Work of Several Authors 


p 

Range of 

Source 

2.38 

-1 - -200 

0G0-V observations plus others (Table 4.3, third 
row) 

2.16 

zb - -200 

0G0-V Group A. Eq. (4.3) 

2.3 

5 - 16 

Blackwell and Pet ford's Observations (1966) 

2.57 

-1*0 

Newkirk's presentative corona (1967) 

2.28 

~70 

Whang et_ al . , solar wind model (1966) 

2.09 

-100 

Cuperman and Harten's solar wind model. 
(1971 Fig- 2,) 

2.12 

( V 
LO 
O 
O 

Hartle and Sturrock's solar wind model. (1968, 
Fig. 7) 

2.25 

l V 
H 
O 
O 

Noble and Scarf's solar wind model (1963) 

2.09 

-200 

Neuman and Kopps ' streamer model (1969?, 1970?) 

~ 3 

153 - 214 

Neugebauer and Snyder's Observations (1966b) 

7 

1.2 - 2.2 

Gillet et_ al. , Observations (1964) 

1.94 

5-70 

Erickson's Observations (1964) 



-95- 


Most theories, if not all, predict something close to an inverse power 
law at distances larger than about 100 R . 

From Table 4.1 we can compute the extreme values of p deduced from 
0G0-V observations. They are 1-35 and 9-52 derived from events No. 29 
and No. 34, respectively. We consider that the last is too high, 
possibly due to an error in the determination of a. More meaningful are 
the p values corresponding to the values of Ct that are one standard 
deviation from the mean index given in the relations, Eq. (4.3); the 
resulting values of p are 1-55 and 3-70. We see that they cover the 
range shown in Table 4.5. The exception is given by the observations of 
Gillet et al., probably because they correspond to distances very close 
to the Sun* 

As we noticed before, the measured a is never = 1. (See Fig. 4.3). 
This should be expected because, examining the Eqs. (A. 11) and (A. 12), we 
see that p and D approach infinity as a approaches 1 and that Eq. (A.l4) 
is undeterminate at that value. For a = 1 Eq. (4.12) gives N = 0 and 
for a = 1 the electron density increases with increasing distance, all of 
which is absurd. 

In order to study further Form 1 we attempted to fit Eq. (4.12) 
to some observations and models by selecting three points from the data and 
determining the parameters A, b, p such that the curve would pass through 
those three points. 

We obtained the solution by computing the functions Q (b) and 

l 

Q (b) derived from Eq. (4.12) and defined below: 
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where N is the density at distance r , etc. 

1 i 

The solution for both b and p was found graphically by plotting 
and as functions of b, and by finding the intersection: 


Q (b) = Q (b) = p . 

1 2 


Table 4.6 shows sets of values of A , b and p that fit the results 
of several authors with different degrees of precision. Except for 
Whang et al. model the rest were done by the method explained above. 

The selection of the three points was by visual inspection and we did 
not attempt to optimize the fitting. In order to have an estimate of the 
quality of the fitting we calculated at different distances the ratio 
n/n', of the electron density obtained from the fitted model to that of 
the original data. Table 4.7 shows these ratios for some of the data 
included in Table 4 . 6 . The fitting to observational data is usually good, 
for example, that of Newkirk (1969), Newkirk (1970), Blackwell and 
Petford (19 66 ), van de Hulst (1950). The fitting to models in the range 
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Table 4.6 

Fitting of N = A/(r - b)^ 


p 

b 

log A 

2.55 

O.89 

6-33 

4.66 

O.52 

8.34 

2.86 

0.60 

6.87 

5-31 

0.4-1 

7.32 

5-35 

0.06 

7.83 

2.23 

0.02 

6.11 

2.27 

1.72 

5-77 

2.33 

0.93 

6.27 

3-51 

0.68 

6.87 

2.25 

0-79 

6.33 


Authors 

Newkirk (1967) 

Newkirk (1969) 

Newkirk (1970) 

Newkirk (1970) 

Newkirk (1970) 

Blackwell and Petford 
(1966) 

Pneuman and Kopp 
(1969?, 1970?) 

Whang et al . (1966) 

van de Hulst (1950) 

Hartle and Sturrock 
(1968, Fig. 1) 


Remarks 

Composite representative 
corona. Equator min. 

Average active enhancement 
plus quiet component. 

Southwest streamer. 

Southeast streamer. 

Equatorial excluding 
streamer. 

Corona, r < 20 R. 

© 

Streamer core model 

r > 2 R . 

- © 

Solar wind model. 

Equator max. 


Solar wind model. 
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Table 4.7 

Values of n/n ' Representing the Quality of Fitting ^ 


Models from Observations Models from Theory 




Newkirk 

Blackwell- 

van de Whang 

Pneuman- 

Hartle- 

r 

Newkirk 

Newkirk (1970) 

Petford 

Hulst et al. 

Kopp (1969,? Sturrock 

<V 

(1967) 

(1969) S.W.Str. 

(1966) 

(1950) J 19 SS) 

1970 ?) 

( 1968 ) 

1.00 




1.00 * 


2.34 

1.02 

1.00 >* 



1.15 



1.12 


1.00 '* 





1.25 


1.00 





1.37 


0.95 





1.50 


1.00 * 





1.6 

0.51 






1.66 


1.00 * 





1-75 


0.95 





1.98 


1.03 





2.00 


1.00 * 


1.00 '* o. 48 ' 

10.00 

1.00 * 

2.36 


1.00 * 





2.56 


1.10 





3.00 




1.12 

1.4l 


3.16 






1.00 

3*17 


1.00 





4 

1.00 * 



1.00 0.89 

1.00* 


5 



1.00 * 

1.00 * 



6 

1.09 



o.73 



7 



O.98 




10 

1.02 


1.00 * 

1.12 

1.00* 

1.00 * 

13 



1.01 




16 



1.00 r* 




31.6 






1.00 

4 o 

0.84 




1.16 


100 

0.86 





0.99 

200 




0.89 

1.00* 

1.00 * 

215 

1.00 * 






1000 




0.63 


0.63 


N - electron densities computed from models using A, b and p listed in 
Table 4.6 

N ' = electron densities from authors listed in Table 4.6 
* = Indicates the points used in the fitting 
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from 1 to 215 is usually poor at small distances but, in general, for 
r greater than 2 or 4 R^ the agreement is satisfactory. For theoretical 
models, such as those of Whang ejt al. (1966) and Hartle and Sturrock 
(1968), no attempt was made to have a good fit beyond 215 Rq* In fact, 
we observed that the corresponding values of p in Table 4.5 and those 
used in the fitting model, Table 4.6, are not the same, therefore we do 
not expect a good fit at large distances. This suggests that an alternative 
way to find the three parameters that will give a good fit to a theoretical 
model is to adopt p as predicted from the model and then determine A 
and b by a fit at two points. The agreement of Eq. (4.12) with obser- 
vations and theory over a large part of the distance range, and its 
simplicity makes it useful for analytical work. 

Several authors have in the past fitted their observations to the 
first few terms of a sequence of inverse power of r (Stelzried et al., 1970; 
Baumbach, 1937 5 Hollweg, 1968; Blackwell and Petford, 1966; 

Ingham, 1961; Saito, 1950; Newkirk et al., 1970; Newkirk et al., 1969; 

Allen, 1947a)or else they have used other functions that are susceptible 
of being expressed in that form (Bohlin and Simon, 1969; Newkirk, 196l; 
Bohlin, 1970; Chapman, 1957; Wild et al. , 1959)' It seems that the 
simplicity of the formula Eq. (4.12) and of the procedure for determining 
the three parameters is advantageous even if the fitting is not so precise. 

Electron Density Models . The next step in the study of Form 2 is 
to adopt values for A and b that will lead to the formulation of some 
numerical models . 

Since we want to compare the computed values of drift rates with 
the observed values, it seems natural to select an exponent p obtained 
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from the radio data# Sines we also want our density model "to cover "the 
widest range of distance we should adopt the p derived from the largest 
range of frequency. Earlier we found it to be p = 2.38, Eq. (U. 13 ). 

After p has been chosen the values of A and b can be obtained 
by assigning electron densities at two distances. The obvious choices 
are 1 and 215 solar radii. From Table U .6 we see that b is expected to 
be of the order of unity, therefore at 1 A.U. r » b and A will be 
determined by N at the Earth, N(215). Consequently b will be effectively 
determined by N near the Sun, N(l). 

Observational facts and the acceptance of the local plasma 
hypothesis impose some restrictions on the numerical values of N(l) and 
of N(215 ) . Because Type III bursts have been observed at frequencies 
as high as 550 MHz (Maxwell et al . , i 960 ) the electron density near the 
Sun cannot be smaller than 3-73 x 10 9 cnf 3 . (See Table 4.8). The 
corresponding plasma level is not known and for simplicity we will take 
it at 1 R 0 . Within the wide range of values of N measured near the Earth 
we choose two alternatives: one with N = 7.7 cm " 3 and the other with 

N - 30.9 cm 3 . The values are within accepted observational limits and they 
correspond to a plasma frequency of 0.05 MHz as second harmonic and 
fundamental, respectively. 

The choice of the electron density at the Earth has implications 
regarding the homogeneity of the corona that we want to model. If we 
accept an homogeneous corona then the maximum eledtron density at the 
Earth compatible with the observation of a burst at 0.05 MHz is 30.9 cm -3 
because the Earth would be precisely in the plasma level. This is true 
in the case in which the burst is supposed to originate in the homogeneous 
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Table U.8 

Plasma Frequencies and Electron Densities 
(f = 9 x 10 -3 JN^ 2 ) 


P 37 (cm 3 ) 

\ 



J = 1 

.1 = 2 

550 

3.73 x 10 9 

9.3*J x 10 

3.5 

1.15 x 10 s 

3-78 x 10 

1.8 

k.O x 10 4 

1.00 x 10 

0.90 

1.00 x 10 4 

2.5 x 10 3 

0.6o 

*kUi. x 10 3 

1.11 x 10 

0.35 

1.51 x 10 3 

3.78 x 10 : 

0.20 

U.93 x 10 2 

1.23 x 10 : 

0.10 

123.5 

30.9 

0.05 

30.9 

7.7 
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maximum-density corona itself or in the case of a denser streamer 

embedded in that corona. Evidently in the first case there are more 

stringent limitations in directionality due to refractive effects. 

Some models based on the above conditions are shown in Table 4.9. 

Model T (thin) is not dense enough at the base of the corona to give 

a plasma frequency of 550 MHz, (see Table 4.6). The designation D stands 

for Dense. Model C (Combination) is as dense as D at r = 1 and as 

T at r = 215. D*2 corresponds to a model twice as dense as D. Because 

of the above discussion model D*2 probably represent a denser 

streamer. Figure 4.8 shows some of these models and some observational 

data. On the frequency scale f ? represents a fundamental at frequency f 

emitted at level r(f) and f , a second harmonic at frequency f emitted 

n 

at the level r(f/2). 

Table 4.9 

Characteristics of Some Models. N(r) = A/(r - b ) P 


Model 

N(l) 

(cm 3 ) 

N(215) (cm -3 ) 

log A 

b 

p 

T 

4.0 

X 

H 

0 

CD 

7-7 

6.44 

0.88 

2.38 

C 

4.0 

X 

10 9 

7-7 

6.44 

0.95 

2.38 

D 

4.0 

X 

10 9 

31 

7.05 

0.91 

2.38 

C*2 

8.0 

X 

10 9 

15.4 

6.75 

0.95 

2.38 

D*2 

8.0 

X 

10 9 

62 

7.35 

0.91 

2.38 


Direct determinations of the heights of the Type III sources have 
been made by several authors. Table 4.10 shows the results found in the 
literature, and Fig. 4.9 compares them with the predictions from several 


of our models . 



(crrT 3 ) 
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FIG. 4.9 OBSERVED AND PREDICTED HEIGHTS OF TYPE III SOURCES. 
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Table 4.10 

Direct Measurements of the Height of Type III Solar Bursts Sources 


f (MHz) 

R <V 

Author 

Instrument^ 

200 

1.4 

Morimoto (1961) 

I 

169 

1.3 - 1.6 

Bourgeret et al . (1970) 

R 

52.5 

2.07 ~ 0.10 

Morimoto (1964) 

SFI 

50 v ' 

~1.7 

Kundu et al . (1970) 

SFI 

60 

-1.66 

Wild et al. (1959) 

SFI 

45 

~2 .0 

Wild et al. (1959) 

SFI 

(*) 

30 

~2 . 3 

Kundu et al. (1970 ) 

SFI 

26.3 

OJ 

0 

+1 

OJ 

Malitson & Erickson (1966) 

I 

26.3 

2.0 - 2.5 

Erickson (1963) 

I 

19- 7 

2.9 

Shain & Higgins (1959) 

I 


Fundamental. 

I means interferometer, R means radio heliograph and SFI means sweep- 
frequency interferometer. 

We can see that direct observations require a model with a distri- 
bution of density between two and four times that of model D. The last 
model corresponds probably to streamers . At the lower frequencies the 
observations fit the less dense model D*2 better, and we suspect that the 
higher frequency bursts may originate in local condensations which may 
disperse and blend into a homogeneous corona at greater distances . 

It has been customary to refer coronal electron density distributions 
to a model considered as standard. Some of the reference models have 
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been Baumbach-Allen , van de Hulst (1950) and Newkirk (1967). We will 
adopt the last one. In Table 4.11 we give an approximate equivalence 
expressed in terms of the ratio of electron densities of the different 
models (N) to that of Newkirk's (N^). 


Table 4.11 

Approximate Equivalent of Some Electron Density Models 


Model 


Oj/N ) Between 
10-40 (R Q ) 


N (215) Plasma Level of 0.05 MHz 
(cm' 3 ) Fundamental. (R ) 


Newkirk (1967) 

2 

T 

4 

C*2 

4 

D 

8 

D*2 

16 

Fainberg & 
Stone (1970) 

16 

Penman & Kopp 
(1969?, 1970?) 

( Interstreamer) 

0.20 

Pneuman & Kopp 
(1979?, 1970?) 

0.20 

(Streamer's core) 


( 1" ) 

Extrapolated. 



2.5 

85 

7.7 

115 

15.4 

155 

31 

215 

62 

300 

to (t) 

230 

1.1 

44 


3.1 71 


Prediction of Frequency Drift Rates ( FDR ) . The next step in the 
study of Form 2 is to investigate if it predicts results compatible with 
radio observations. For this we will predict FDR and we will compare 
them with the observations. 
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Due to the simplicity of Eq. (It. 12) it is easy to derive an 
expression for the frequency drift rate observed at the Earth. The details 
are given in Appendix 3. The final expression for the drift rate given 
by Eq. (A. 26) is : 


log 


df 


peg 

dt 

- log 

_ 2R 0 (8lj 2 10" 6 A) 1/p _ 


p + 2 


log f - log [H(f, e, e p ) + 6 G(f, £, 0 p )] 


(A. 26) 


The terminology has the same meaning as before except for the 
functions H and G which are defined in Appendix 3. In general, H is 
related to the trajectory of the exciter particles and G is related to 
the trajectory of the electromagnetic waves. The frequency is expressed 
in MHz, the time in seconds and the distance is solar radii. 

In order to compute numerical values for the frequency drift rates 
we need to adopt values for the rest of the parameters involved in 
Eq. (A. 26). 

The velocity 6 and the mode of emission j will be kept as variable 
parameters . 

For the geometry of the plane trajectory we will consider two 
alternatives: a radial line and an Archimedes spiral with e equal to 0 

and tt/ 2, respectively, in Eq. (4.8). The values of £ can be substituted 

directly into the Eq. (A. 22) to compute the functions H and G. The 

> 

position 9 of the associated flare will be kept as an adjustable 


parameter. 
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To illustrate the effect of the orientation of the trajectory on 
the FDR curves we have chosen the following values for Q : 90% 30% 0°, 
-90°. We recall that 90° corresponds to a west-limb flare. The angle 
of 30° was selected from the study of the distribution of the helio- 
centric longitudes of the flares associated with our bursts. As we 
will see in Section 4.5-l } the 9^, distribution of 0G0-V bursts has an 
average of 30° with a standard deviation of 51°. We shall assume that 
this is the case also for the bursts outside the frequency coverage of 
0G0-V in order to use all the observations. 

To summarize, the following elements were used to compute 
log | df / dt I versus f from Eq. (A. 26): 

Models: C, D and D*2 (see Table 4.9) 

Trajectory: radial (e = 0) and spiral (e = n/2) 

Mode: Fundamental (j = 1) and second harmonic (j =2) 

Position of associated flare: 9_ = 30°, also 0, 90° and -90°. 

r 

Figures from Figs. 4.10 to 4.l4 show 'Curves representative of 
the results. We should compare them with Fig. 4.5. in general we 
observe that the resulting curves are fairly straight down to about 
0.2 MHz. In this range their slopes are equal or close to what is 
predicted, this means that the term A in Eq. (A. 27) varies slowly. 

Special effects originating in A are predicted at frequencies below 
approximately 0.2 MHz and we will discuss them shortly. 

In order to illustrate the effect of the different parameters 
on the FDR curves, only one parameter was changed in each of the figures. 

Figure 4.10 shows the effect of varying 9 , i.e., the orientation 

r 

of the trajectory with respect to the Sun-Earth direction. The curves 
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FIG. 4.11 EFFECT OF THE SHAPE OF THE TRAJECTORY ON THE PREDICTED 
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FREQUENCY (MHz) 

FIG. b.lk EFFECT OF THE ELECTRON DENSITY MODEL ON THE PREDICTED 


FREQUENCY DRIFT RATE CURVES. 


are fairly straight except at the lower frequencies where a change 
in shape occurs. This effect is due partly to the proximity of the 
source of the burst to the Earth. For the case illustrated in this 
figure the trajectory for 9^ = jt/2 passes, by construction, through 
the Earth and because of the model chosen the 0.05 MHz as fundamental 
occurs at 215 Rq • For this particular case there is a sudden discontin- 
uity in the FDR because when the exciter particles pass through the 
Earth the term A in Eq. (A. 27) changes sign. Physically this is due 
to the fact that before the particles go through the Earth with a 
vectorial velocity ~v the radio waves that we observe propagate in a 

j 

direction very close to v* but just after the particles pass the Earth 

the radio waves that we observe are those propagating in the direction 

-v. This is discussed in the Appendix No. 5- The sharpness of the 

discontinuity depend on how close to the Earth the trajectory passes. 

The possibility of observing this effect is low because the' large size of the 

source, in which we would be immersed, would smear out the time of 

arrival. It is of interest to notice that the separation of the 

curves shown is small considering the wide range of values of 9 , 

F 

between -jt/2 and -bt/2. We recall that for OGO-V bursts one standard 

deviation in the distribution of © in heliocentric longitude is only 51 ° • 

r 

Figure 4.11 illustrates the effect of varying the shape of the 
trajectory. Under equality of the other conditions the radial and 
spiral paths are coincident down to approximately 0.2 MHz. If we 
extrapolate the straight line defined in the high frequency range and 
use it as a reference then we can have a quantitative idea of the 
amount of bending of the computed curves. We will define f^/2 as 
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frequency at which the drift rate drops by a factor of 1/2 from the 
straight line extrapolation. For the conditions of the figure the 
frequencies f ^ for the spiral and radial models are O.083 MHz and 

< 0.025 MHz, respectively. The drift rates predicted at 0.05 MHz 
differ by a factor of 2.2. This effect might be detected by observations 
obtained with high time resolution. 

Figure 4.12 illustrates how the velocity affects the drift-rate 
curves. Above 0.2 MHz they are fairly straight and parallel. As we 
could expect the larger the velocity the higher the curve. We observe 
that the amount of curvature in the low-frequency range depends also 
on the velocity. For the conditions of the figure the frequencies f-^ j g 
are, in order of increasing velocities, 0.078, O.O85, and 0.100 MHz. 

Figure 4.13 shows the effect of the mode of the emission. We 
should note that the frequency scale is the same for both curves but 
the corresponding distance scale is not. In fact, for the conditions 
of the figure the fundamental radiation at 0.05 MHz occurs at 215 R 0 
while the same frequency as second harmonic occurs at about 390 R Q . The 
amount of departure from a straight line is about the same for both 
except that the curve corresponding to the second harmonic starts bending 
down at a higher frequency. For the conditions of the figure the ^^_/2 
values for the fundamental and second harmonic are O.O85 and 0.180 MHz, 
respectively. It is of interest to notice that the two curves are 
practically parallel, the lower one corresponding to t he second harmonic^ 
consequently, in fitting observations this mode leads to higher velocities, 
as expected. We will come back to this point when we discuss the 
phenomenon of the prevalence of the second harmonic radiation. 
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Finally, Figure 4.l4 illustrates how the overall density affects 
the FDR curves. As expected, the denser the model the higher the 
frequency at which the curves bend down. The f values for the 
conditions of the figure, in order of increasing densities, are 
0 • 045 , 0.084 and 0.126 MHz. The curves are fairly parallel and are all 
practically straight lines above 0-300 MHz. The figure shows the effect 
of doubling the electron density of model D. For reference we recall 
that for distances corresponding to frequencies less than 25 MHz 
model C is four times less dense than Model D. 

In Section 4.2.1 (Table 4. 3) we found that the scatter of the 
observational points in Fig. 4.5 from the least-square fit line has a 
standard deviation a =0.13 in log j df/dt | (a factor of I.35 in the 
drift rate). It is interesting to express the effects on the FDR curves 
resulting from the variation of the different parameters in terms of 
this standard deviation. This information is contained in Table 4.12. 

is the coefficient A in Eq. (4.12) corresponding to Model D. We 
observe that the position of the curves in the plot is more sensitive 
to changes in the velocity than in the other parameters. For comparison 
we note that the parameters j, A and P were doubled. 

When we adopt a density model we can select combinations of 
trajectories, orientations and velocities to compute the FDR curves for 
the fundamental and second harmonic. Since we generally know the position 
of the associated flare of our bursts, we could find by trial and error 
the set of other parameters for which a computed curve best fits the 
observations. If the observations are adequate we could' obtain the velocity 
the trajectory and the modes. However, this method would apply only to 
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! Table ^1.12 

Effects of the Variation of Parameters on the FDR Curves Expressed 

Approximately in Terms of a 

Approximate 

Range Approximate Frequency 

of Separation Range of Validity 


Parameter 

Parameters 

of Curves 

(MHz) 

Figure 

9p (Orientation) 

it it 

2 > ~ 2 

2a 

all 

4.10 

e (Trajectory) 

0 - — 

u 2 

< 2a 

> 0.10 

4.11 


(radial - 

spiral) 



j (Mode) 

1-2 

2a 

> 0.50 

4.15 

A (Model) 

A D ' 2 Aj 

la 

> 0.50 

4.14 

(3 (Velocity) 

0.2 - 0.4 

< 4 a 

all 

4.12 


- 118 - 


those bursts for which there is no uncertainty about the measured drift 
rates. This restriction limits considerably its applicability because 
it is indirect; it involves derived quantities (FDR) obtained from 
measurements (times of arrival). A more direct method is to fit directly 
the times of arrival. This has the advantage that it can be applied to 
complicated cases containing mixtures of fundamental and second harmonic 
emission modes. This approach was used and it will be discussed in 
Section ^.2.5. 

Summary 

We have found a mathematical expression to represent a coronal 
electron density model that has the following properties: 

(a) It is of simple form, 

(b) It contains only three parameters that can be determined 
from the density gradient and the number density at any two distances, 

(c) It can be fitted to individual sets of optical observations 
by an adequate choice of parameters, 

(d) A selection of the appropriate parameters produced models 
that are consistent with optical observations near the Sun, with 
measurements in the Earth neighborhood, with radio observations of 
Type III bursts, and with direct measurements of the height of Type III 
sources, 

(e) Its gradient at large distances agrees with those predicted 
by theoretical models, and 

(f) A derivation of the frequency drift rates of Type III bursts 

predicts some effects that may allow testing of some of the assumptions 
made in this study. 
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Electron Density Di.stri but ion - - Other Forms . A function that 
resembles Eq. (4. 12) in its behavior is 



(It. 14) 


A comparison of the series expansion of Eqs. (4.12) and (4.l4) 
reveals that under the normal conditions r^ 3 > b and p > 1, the latter 
function decreases faster at small r, otherwise the behavior of Eq. (4.l4) 
and (4.12) are very much the same. The function Eq. (4.l4) was not used 
in this work. 

In general the theoretical models predict a decrease in density 
slower than that given by Eq. (4.12). (See for example Table 4.7.) 

A modification of Eq. (4.12) to obtain this effect, at the expense of 
increasing the complexity, is: 


P q , 

! r - mr - b 

where m > 0 and p > q. This form was not used either. 

4. 2. ^ Determination of Velocities and the Second Harmonic 
Phenomenon . 

Introduction . In the preceding two sections we have discussed 
the electron density models evolved from the present study. We have 
explained the deficiencies of Form 1 and how this led to Form 2 which 
was finally accepted. 

Because it was developed first and because under appropriate 
conditions it provided reasonable electron density distributions. 
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Form 1 was used in the next chronological stage of our research. This 
consisted in predicting the times of arrival of the bursts at the 
different frequencies and comparing them with the observations. Later, 
after Form 2 was developed, we found that the early model used in the 
prediction of the times of arrival and based on Form 1 was numerically 
close to Model C (Form 2), at least within 1 astronomical unit. We 
estimated then that the use of any of the two models would not modify 
appreciably the deduced velocities and that the presence and character- 
istics of the second harmonic phenomenon, described shortly, would be 
unchanged. Based on these considerations we decided that it was not 
necessary to repeat the work using a new model. We will start then 
by discussing briefly the model used. 

The Electron Density Model . We used Eq. (A.l4). The selection of 
the parameters was done as follows: The parameters a and a, related to 

p and D by Eqs . (A.ll) and (A. 12), were chosen from an earlier fit of a 
straight line through the observational points in Fig. 4.5. These points 
consisted of the November 18, 1968 event observed by OGO-V, and of the 
observations of other bursts made by various workers . The values adopted 
were a = 1. 89 ..and a = -I.97. 

Comparing them with those in Table 4.3, we see that they are close 
to values determined subsequently in a more precise way. The adopted 
a leads to p = 2.25 that should be compared with p = 2.38, Eq. (4.13). 

We computed models for j = 1 and j = 2. 

The trajectories selected were spiral and radial with e equal to 
n/2 and 0, respectively, in Eq. (A.5). Equation (A. 10) gives the 
explicit form of s for e > 0. For 9 , the orientation of the trajectories. 
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which appears in the term x in Eq. (A A), we selected jt/2 and 0. We 
applied the boundary condition at r Q = 1 assuming that N q = 3 x 10 s cm 3 . 

The fact that this density is not high enough to produce a plasma 
frequency of 550 MHz does not invalidate our results because we use 
the model only in connection with the OGO-V data corresponding to much 
lower frequencies. 

Finally we let P be a running parameter. We proceeded as follows: 

We computed N = N(r) for different P but for a given set of N , r , 

j and trajectory. Since N q was fixed the result was a set of curves 

converging at r = 1 and that could be pivoted about this point by 

varying p. By controlling the position of the curve with p in this 

fashion we selected a value of P such that the corresponding curve gave 

a good fit to the optical observation out to 20 fh and to the spacecraft 

© 

observations near the Earth. By hypothesis the model is consistent 
with the drift rates of radio observations. Figure 4.15 shows curves 
for N q = N(l) = 3 x 10 s cm 3 , j = 2 , P = 0.30 and a spiral trajectory. 

This figure illustrates also the effect of the orientation of the spiral. 

In order to select one model we made a rough average of the two curves. 

This is equivalent to choosing a 9 probably somewhere between 30° and 

r 

50°. We will see in Section 4.5.1 that 65 percent of our events were 
associated with solar flares on the western hemisphere. The electron 
distribution adopted is shown in Fig. 4.l6. For comparison we include 
some observations. The dashed part of it was not used because it 
corresponds to plasma frequencies > 3-5 MHz. Due to the relatively low 
density this distribution corresponds to a background corona rather than to a 
streamer. This model is comparable to model C of Fig. 4.8, based on Form 2. 
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FIG. 4.3.6 ELECTRON DENSITY MODEL ADOPTED TO COMPUTE TIMES OF 


ARRIVAL 



We also computed electron distributions assuming a radial trajectory 
but we made little use of it. For the sake of completeness we show in 
Fig. 4.17 models corresponding to the same conditions as those in 
Fig. 4.15, except for the trajectory. 

Prediction of the Times of Arrival . Once we adopted an electron 
density distribution it was straightforward to compute the times of 
arrivals of the bursts at the Earth. For this we used the same. basic 
assumptions discussed in Section 4.2.2. We assumed that the particles 
travel with constant speed along a spiral trajectory (with e = n/2 in 
Eq. 4.8) contained in the plane of the ecliptic. We also assumed that 
at a particular level the plasma frequency and its second harmonic are 
emitted simultaneously from the same point with the same propagation 
velocity in a straight line toward the Earth. The expected sequence of 
events is as follows: The first burst to arrive at the Earth is the one 

emitted at 5.5 MHz as a fundamental; its second harmonic at 7.0 MHz falls 
outside the range of our radiometer. As the particles proceed in their 
motion they reach the 1.8 MHz fundamental level; its second harmonic at 
5.6 MHz is not observed because the bandwidth of the radiometer channel 
centered at 3*5 MHz is only 10 kHz. The process continues similarly at 
higher levels. Table 4.13 shows the levels that, according to the adopted 
model, contribute with emission at frequencies observable with our 
instrument. Subscript F refers to fundamental and H to second harmonic 


emission. 
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FIG. 4.17 ELECTRON DENSITY MODEL BASED ON FORM 1. RADIAL 


TRAJECTORY 
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Table U.13 

Levels that Contribute with Radiation 
Observable in OGO-V Radiometer 


r( V 

(MHz) 

f (MHz 

n 

3.75 

3.5 

— 

5-95 

1.8 

— 

6.00 

— 

3.5 

10.3 

0.9 

1.8 

Ik. 9 

0.6 

— 

19.0 

— 

0.9 

23.6 

0.35 

— 

27-0 

■ — 

0.6 

38.1 

0.2 

— 

k2.k 

— 

0.35 

68 

0.1 

0.2 

11 L 

0.05 

0.1 

178 

— 

0.05 


According to our definition we see that the only levels from 
which we can observe a fundamental and a second harmonic at the same 
frequency are at distances of 10. 3 , 68 and llU R^ corresponding to a 
fundamental of 0.9, 0.1 and 0.05 MHz, respectively. 

Let us refer to Fig. 1*.6. Equation (A. 10) gives the distances 
along the spiral between the center of the Sun and a point at a radial 
distance r. The distance x is given by Eq. (A. 4). Then the time 
elapsed since the exciters leave the surface of the Sun with constant 
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velocity 3 until the fundamental emission at frequency f , excited at 

F 

the level r, is received at the Earth is: 

\ R © 

+ » (4.15) 

where t(f^,) is in sec, R^ in Km and c in Km sec 1 . The relationship 
between r and f is given in Table 4.13. A similar expression can be 
obtained for the second harmonic at the same frequency by replacing F 
by H. Since the first frequency observed in our radiometer under the 
above conditions is 3.5 MHz emitted as fundamental, we will take its time 
of arrival as reference and work with time delays. 

We computed and plotted time delays versus frequency for different 
e F and 3. We varied 0 F between -90° and +130° in steps of 10°. For 
each 0 F we varied 3 between 0.10 and 0.60 in steps of 0.02 and sometimes 
0.01. An example of predicted times of arrival is shown in Fig. 4.l8. 

In what follows we will use the notation fF or fH to identify 
points in this plot, f will be the frequency in MHz of a channel of our 
radiometer. F and H have the same meaning as before. 

The Second Harmonic Phenomenon . In Section 4.2.1 we described the 
technique used first in order to determine the frequency drift rates. This 
consisted of drawing a smooth curve through the observed points in a 
frequency versus time of arrival plot and then measuring the slop of this 
curve. We explained that this method produced FDR curves that could be 
separated into several groups according to their shapes (see Fig. 4.2). 
Except for those curves that were a straight line, the shapes were 
difficult to explain in terms of reasonable physical arguments. 
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FIG. 4.18 PREDICTED ARRIVAL TIMES. 
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The second technique that we devised consisted of comparing directly 
the observed times of arrival with the predicted ones. For this we 
selected the set of predicted plots corresponding to the © of the 

r 

associated flare of the burst under study and by overlapping with the 
observational plots (Fig. 4.1) we selected the $ within the set that would 
give the best fit. Note that in this way we determine the modulus of the 
velocity and not the radial component. 

We first used this technique by applying it to the time of start of 
the first burst of each event. In some cases the observed points would fall 
on the positions predicted for either the fundamental (F) or second 
harmonic (H). In other cases the observations would fall on the (F) above 
a certain frequency and on the (H) points below that frequency. Figure 4.19 
shows a typical example in which curves have been drawn through the predicted 
points. We can see that the observations at 3-5» 1.8, 0.9 and 0.6 MHz 
fit the predicted arrival of fundamentals while the 0.35 MHz fits the 
arrival of the second harmonic at this frequency. This pattern gave the 
clue to the peculiar shapes of the spurious FDR curves obtained by the 
first technique, that is, that we were indiscriminately passing a curve 
through points that corresponded to emission in two different modes and 
produced at two different levels in the corona: a fundamental and a second 
harmonic at the same frequency. The disappearance of the observed fundamental 
and the emergence of the second harmonic emission at the same frequency has 
been called by us "the second harmonic phenomenon." 

A completely developed single burst observed in our eight channels 
would give a l6-point plot. In general, however, some of them may be 
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FIG. 4.19 COMPARISON OF THEORETICAL AND OBSERVATIONAL ARRIVAL 


TIMES OF A TYPE III BURST. 
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missing either because the corresponding emission was too weak for 
physical reasons, or because of instrumental reasons (interference, 
data gaps, etc.) Figure it. 20 shows a case in which the 1.8-MHz 
fundamental point is missing and so are the other fundamentals below 
0.6 MHz. The combination of a few missing points would evidently 
produce very peculiar FDR curves by the first technique. In fact, we 
recall that there was a group of bursts that could not be classified 
because no simple curve could be passed through the observed arrival 
times of the leading edge of the time profile at the different frequencies 
The bursts that produced a fairly complete set of either F or H points 
would have a FDR that is a straight line : those are the ones that we 
classified earlier at Type A. Now, we find that the drift rates derived 
from the fundamental and from the second harmonic curves are practically 
straight lines with the same slope in the classical log-log plot (Fig. 
^.13), regardless of the value of 9 p . However, the value of the apparent 
drift rate A.#At is smaller for a west limb flare than for a center disc 
flare because of the propagation time of the electromagnetic radiation. 
This means that when we draw a straight line segment through the drift 
rate measured in this range of frequencies it is going to be steeper 
for the larger 9 p . This indeed was found. 

We have mentioned in Section h.2 of this work that most of the 
events studied (if not all) were extremely complex. At each frequency 
they were made up of numerous burst components that in general could not 
be separated mainly because of lack of time resolution. We explained then 
that for each event we made a chart displaying frequency versus time of 
arrival for all the components that were unequivocally present or whose 




existence was suspected on good grounds. When using the second technique 
there exists then the possibility of confusion in the sense that at a 
certain frequency the fundamental of a late component can be misinterpreted 
as the harmonic of an earlier component. The situation is complicated 
when some observational points are missing. 

The first step in the use of the new technique was to apply it 
only to the beginning of the events. The second step was to apply it to 
all the components of each event, as already explained. We identified 
the components of an event by Roman numerals according to the chronological 
order of arrival. The operating procedure (based on the overlapping) 
was as follows : We first fit the earliest times of arrival of an .event 

at the different frequencies to predicted F or/and H points. For example, 
we labeled points as IF and IH to indicate that they belonged to the 
first component and were identified, respectively, as fundamental and 
second harmonic. We wrote down the corresponding 6. We moved in time 
and took the next set of points to arrive within the same event; we 
repeated the same fitting operation except that these points were labeled 
II. The process was continued until we reached the last points to arrive 
in the events . 

In an ideal situation all the observed points of an event would 
be accounted for in the fitting and a velocity would be determined for 
each component. However in practice several things can happen to introduce 
confusion : 

(a) Some points of a component may be missing, 

(b) There may be data gaps. 
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(c) A point can be assigned to one or more components as either 
F or H if the velocities of adjacent components are so different that 
the curves cross each other, 

(d) When the uncertainties in time of arrival are large the 
predicted points for different velocities and components fall within the 
error har and therefore are indistinguishable. 

Table l».lU shows an actual example of the way in which we tabulated 
our results. For this event we determined 9 components. The first is 
uncertain and it was assigned only one point at 3*5 MHz; this could not 
be fitted with the rest into a single component:. Component II starts 
as F and then goes to H; it is not seen below 0.9 MHz. For large velocities 
the separation in time between 3-5 F and 3. 5 H is of the same order of 
magnitude or less than our uncertainty in time resolution, therefore, the 
predicted F and H may fall within the time error bar. This happened 
to component III at 3.5 MHz. At 1.8 MHz the same observed point can be 
accounted for as III 1.8 H or IV 1.8 F. This is indicated by encircling 
and connecting them as shown. In component V we can fit a point to 
3.5 F and another to 3.5 H. The predicted F at 1.8 MHz falls in a 
calibration cycle therefore its existence cannot be either accepted or 
rejected. 

What we have called the structure of the F-H pattern is a 
combination of the sequence, number and mode of the observational points 
assigned to a component. If we do not impose some criteria on this we 
may have a large amount of possible fits for the points of an event. 

The main criterion was to find fits using the least number of components 



Example of Tabulation Resulting from the Fitting of Times of Arrival (Event No. 1) 
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with the smallest spread in the velocities. This last requirements seems 
reasonable if the components come from the same flare within a few 
minutes of time; but there is no guarantee that it is always true however. 
In deciding to what component to assign the observation corresponding to 
ambiguous points (those encircled and connected in Table lull*) we tried 
to avoid those resulting sequences in which the modes alternates between 
fundamental and harmonic as the frequency decreases, and also those that 
were too incomplete. In selecting a representative velocity for the 
whole event we gave more weight to the velocities of the components fitting 
the lower frequencies. For the event No. 1' illustrated in Table k.lk we 
selected 3 = 0.36. 

The above conditions usually lead to a more or less unique solution 
for each event. However there are some cases that admit two solutions. 

This is due to the fact that the predicted times of arrival of the second 
harmonics for a burst with an associated velocity 3 are very close to 
the times of arrival of the fundamentals with a velocity approximately 
equal to 3/2, the other parameters being the same. An example of this 
is shown in Fig. A. 21, where at the higher frequencies there is practically 
an overlap between the H curve for 3 = 0. 31 and the F curve for 3 = 0.18 
leading to ambiguity in the velocity. At the lower frequencies where the 
coincidence starts to fail, the lack of time resolution prevents us from 
resolving the ambiguity. The structure of these cases is a sequence of F 
points with an assigned velocity 3/2 or a sequence of H points with a 
velocity close to 3. We called these cases "ambiguous" and we will discuss 

them shortly. Not surprisingly we found that most of the ambiguous cases 
had been classified earlier as Type A. 



137 ' 



(ZHI/M) ADN3n03id J 



- 138 - 


The method explained above is not completely objective. In order 
to reduce subjectivity, the fits were performed by two persons. The 
results were not significantly different initially and the cases of 
discrepancy were discussed until an agreement was reached. 

In order to have quantitative information with respect to the 
second harmonic phenomenon we treated the F-H distribution of the events 
on statistical bases. For this we counted at each frequency of obser- 
vation, the number of points in all the unambiguous events that were 
fitted either as F or H, then we analyzed the result according to 

two criteria. Following the first and less rigorous criterion we 
counted only the F and H points that presented no major uncertainties 
in their identification. The count was done regardless of the F-H structure 
except subjected to the condition that a point, to be counted, should 
belong to a component with at least 3 points. The second criterion 
consisted of counting the points of only well defined components with 
at least 4 points . 

Calling F and H the number points counted as fundamental or harmonic, 
respectively, the results were expressed in terms of the fractional content 
of components fitted as second harmonic: 

H 

H + F ’ 

Figure k.22 shows this quantity as a function of frequency for 
the two criteria (see note No. 7 of Table A. 3). The data on second 
harmonic obtained from ground-based observations are very limited and we 
estimate that the fractional harmonic content in the range about 1*0 to 



H/(H + F) (%) H/(H + F) TO 
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fig. ^ .22 FRACTION OF SECOND HARMONIC COMPONENTS IN UNAMBIGUOUS 


TYPE III EVENTS 



2h0 MHz is not larger than 10 % (Wild et al. 195M or less (Smerd et al. 
1962). With regard to this subject we think that Kundu's interpretation 
of the paper by Wild, Murray and Rowe (195*0 data is misleading. His 
statement is: "The harmonics occur in about 60 percent of the bursts 

observed in the U0-2U0 Mc/sec range." (Kundu, 1965, p. 282). The 
statement of Wild et_ al . is : "Of the several hundred Type III bursts 
observed during the year, the recognition of harmonics was considered 
certain on 20 ocassions. In addition there were many doubtful cases. 
Certain recognition is often made difficult owing to the combination of 
wide bandwidth and rapid frequency drift, which may result in the merging 
of the two bands. In cases where the bursts show tails, however, it is 
possible to decide without ambiguity whether or not the two bands are 
present. Of 19 such cases recorded,.' 12 were found to have the two bands. 
If it is assumed that the cause of harmonic production is independent of 
the cause of tail production, this result suggests that the proportion of 
bursts in which fundamental and harmonic are both of detectable magnitude 
may exceed 50 per cent." Kundu's statement apparently refers to the 12 
out of 19 (or 20) special cases. We can estimate approximate upper and 
lower bounds to the ratio H/(H+F) from the above statement by Wild et al. 
It thus appears that during special rare periods the fraction of harmonics 
to total Type III bursts components H/(H+F) can approach 57 <J 0 but over a 
year the average H/(H+F) ~ 20/several hundreds, or < 10%. A statement 
similar to Kundu's has been made by Boischot (1967a). On Fig. h.22 the 
number in parentheses indicates the number of counts. We observe a clear 
tendency of increase in the harmonic fractional contents as the 
frequency decreases. This is the quantitative expression of the second 
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harmonic phenomenon. The 50$ frequency seems to be around 1 MHz. We 
observe that the use of different criteria, does not 
change the results appreciably. It should be emphasized that these 
results were derived from very intense and complex radio events. 

In order to explain this phenomenon we cannot use current theories 
of the emission mechanism of Type III bursts because they have been 
derived for the metric range of wavelengths that is 10 2 and 10 3 times 
shorter than OGO-V's. (Smith, 1969; Melrose, 1970; Zheleznyakov and 
Zaitsev, 1970; Yip, 1970). These theories require observational data 
on the size of the source, on the density of the exciter particles, on 
the distribution of their energies , on the r&tio of intensities of the 
fundamental emission and the second harmonic on the temperature of 
the ambient coronal plasma, etc. Most of these parameters are partially 
or completely unknown at very low frequencies. 

To explain the second harmonic phenomenon several 
possibilities could be considered: 

(a) Directionality effect due to the radiation pattern of the 
source, or to focusing by the ambient plasma or both. In any of these 
cases there should be a dependence on the position of the associated 
flare. In Fig. A. 2 in Appendix No. k, we plotted the distribution in 
position of the associated flares for ambiguous and unambiguous cases. 
Since the last ones were in general those in which the fundamentals 
disappear within the frequency range of our radiometer we could expect 
them to have a different distribution than the ambiguous. As we will 
show shortly there is no evidence that the two groups differ with respect 
to the positions of the associated flares. 
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(b) The fundamental radiation at level r(f) may be actually much 
weaker than the second harmonic at level r(f/2). 

(c) The two radiations in (b) may be of the same order of magnitude 
but scattering and refraction effects of the intervening medium would 
tend to attenuate the fundamental more than the second harmonic at the 
same frequency because the first is emitted at a point further away from 
the Earth and because the ambient conditions at the two levels are 
different . 

As we will discuss in Sections 4.3 and 4.1* the second harmonic 
phenomenon will have important implications in the interpretation of 
the bursts 1 spectra and in the determination of the electron temperature 
of the ambient corona. 

The Distribution of Velocities . The determination of the modulus 
of the velocity of the exciter particles associated with a burst is a 
direct consequence of the fitting of the observed data to the predicted 
times of arrival, together with the recognition of the second harmonic 
phenomenon. As we have explained, we selected a representative value 
for the whole event from the velocities of the burst components. For 
some events we could get a unique velocity but for others we could obtain 
a fit to two solutions with velocities related by a factor near two: 

One solution corresponding to second harmonic radiation associated with 
certain velocity S and the other corresponding to fundamental radiation 
associated with a velocity close to S/2. These cases were called 
ambiguous. Table A. 3 in Appendix No.. 1 contains the list of events and 
their corresponding velocities. S R and S L stand for the high and low 
alternatives of the ambiguous velocities. These values depend on several 
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assumptions, particularly on the electron density model shown in Fig. U.l6. 
There are 32 unambiguous and 3-1 ambiguous cases 9 that is , both types are 
common. 

In an attempt to resolve the ambiguity we plotted the distribution 
of velocities for the unambiguous and ambiguous cases. The first case 
(una mb iguous) is shown in Fig. 4.23 and we denote the velocity g without 
subscript. (See Note No. 7 of Table A. 3). In the second case for each 
event there is a choice between the high velocity and the low velocity that 

we will call g and g respectively. We decided arbitrarily to consider 

H L | 
only the alternatives of all velocities high or all low. The resulting 

distributions are shown in Fig. 4.24. The average velocity and standard 

deviation a of each of the three distributions are shown in Table U . 15 

where n is the number of cases. (See Note No. 7 of Table A. 3). 


Table U. 15 

Parameters of the Velocity Distributions 


Case 

n 

8 

CT 

g i a 

Unambiguous 

32 

0.32 

0.09 

0.32 (1 t 0.3) 

Ambiguous High 

31 

0.33 

0.05 

0.33 (1 ± 0.2) 

Ambiguous Low 

31 

0. 18 

0 . 03 : 

0.18 (1 t 0.2) 


It is interesting that the averages of 3 and are practically 

the same however the spread in individual values is larger in 3 than in 

g by a factor of 1.8. Also the average of 8 is about twice that of 
H 

8 T s in fact the ratio is B /g = 1.83 but the spread in percent is the same. 
L xl v 

This could have been expected because of the following: In Eqs . (A. 26) 
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FIG. 4.24- DISTRIBUTION OF VELOCITIES FOR 51 AMBIGUOUS CASES. 


and its equivalent Eq. (A. 27) let us neglect the correction term A 

that in the range of frequencies of our interest does not contribute 

with more than 20%. The only term left that depends on the velocity is 

the first term in the right-hand side of these equations. The velocity 8 

O / 

and the mode of emission j enter in this term in the ratio 8/j whose 

value is not altered when we take j = 1 with the corresponding velocity 

B or j = 2 with the corresponding 8 . In any case, 

F n 

?H . p2/p 

h 

which equals 1.85 with p = 2.25 (used in the model). This is close to 
1.83 obtained from Table U.15. 

It is interesting to notice in Fig. it. 23 that the unambiguous cases 
show a narrow distribution around 8 = 0.19 » which make the distribution 
look bimodal. The implication could be that there may exist two ranges 
of velocities that center approximately around 0.19c and 0.36c, and 
this would account in part for the distribution of the ambiguous cases. 
However we have not observational evidence to support this hypothesis and 
we know of no theoretical prediction of the phenomenon either. 

In trying to find a clue that would allow us to resolve the 
ambiguity in velocity we ran some statistical tests on the distributions 
obtained and also searched for possible correlations that would 
characterize the unambiguous cases. The details of this work are given 
in the Appendix No. k , here we will only mention it. 
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(a) We applied a y 2 test to the unambiguous, ambiguous -high and 
ambiguous-low distribution in order to estimate the likelihood of the 
hypothesis that each of them came from a normal population with the parameters 
shown in Table 4.15. The probabilities obtained range between 0.40 and 

0.64 and indicate that there is no reason to suspect the hypothesis. 

(b) We did a more detailed study of the unambiguous distribution 
by computing its moments up to the fourth order. The calculated skewness 
is -0.91; that it is negative can be seen in Fig. 4.23. The calculated 
amount of flattening is 2.22, that is, the distribution is more flattened 
than a normal distribution (3-00). 

(c) We applied an F-test to test the hypothesis that the unambiguous 
and ambiguous-high distributions came from two normal populations with 
the same variance (not necessarily the same mean). The probability 

is 0.01 therefore the hypothesis is unlikely. 

(d) We studied the frequency of occurrence versus 9 in unambiguous 

Jr 

and ambiguous cases. We applied a Student's Test and concluded that it 
is likely that the two distributions found are samples of a same normal 
parent population. The probability at the 20% level is 0.10. 

(e) We tried to find correlations involving the ambiguous and 

unambiguous cases. For example, we studied the relationship B versus 

9 for each distribution. These are shown in Fig. A. 3 and A. 4. There 
F 

may be a weak indication of absence of low B at 9 > 60 ° in the unambiguous 

F 

cases, however, the results are considered inconclusive. 

In conclusion, the tests performed gave no statistical evidence 
that would help us to resolve the ambiguity. 


Before we discuss the velocities determined by us we will review 
the results of other workers. Most of them obtained their values by 
assuming different coronal electron density models and other different 
assumptions; therefore it is almost meaningless to compare individual 
results. Table 4.l6 and Fig. 4.25 show data found in the literature. In 
this figure the numbers in parenthesis refer to the authors listed in 
Table 4.l6. In some cases the data has been interpreted by the authors 
in two different ways with completely different values of velocity. 

For those cases we kept the same author number. It is interesting to 
observe in Fig. 4.25 that the values of B concentrate in two regions: from 
0.10 to 0.20 and from 0.30 to 0.50. The fact that we have found the 
tendency to our velocities to concentrate about 0.l8c and 0.36c could be 
the same effect. Our velocities were obtained by assuming a model for the 
electron density distribution. Had we used a denser solar wind model 
within the limitations , we would have obtained higher average velocities 
and the form of the distributions would have been different. We will 
discuss this soon. 

An advantage of our method of velocity determination is that the 
fitting for each burst was done in a frequency range of more than a decade, 
corresponding to a range of distances of about 100 R . 

Other distributions of velocities have been given by Wild et_ al. 
(1959), Mai vi lie (1962a), Stewart (1965), Alexander et_ al_. (1969), and 
Haddock and Graedel (1970). The distribution observed by Wild et al. 

(1959) corresponds to radial velocities determined independently of a 
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coronal model. The distribution shows an abrupt decline in the number 
of bursts for velocities between about 0.35c and 0.2Qc and it shows no 
cases of velocities less than 0.20c. This apparent cutoff in velocity 
has been discussed by de Jaeger (1960a), Malvilie (1962a), Sturrock ( 196H ) , 
Takakura (1966), Melrose (1970). They offer explanations of cutoffs 
between 0.13c and 0.20c assuming that the exciter particles are 
electrons. The distributions obtained by Malvilie (1962a) and by 
Haddock and Graedel (1970) show a faster decline in the low frequencies 
and the cutoff occurs at 0.05c. Alexander et_ al. 's distributions do 
not show an abrupt decline; using their Models I and II they do not 
deduce velocities lower than 0.01c or 0.05c, respectively. 

Effects of the Assumptions on the Deduced Velocities . The 
velocities derived in this section rest upon the validity of the 
assumptions made in the computation of the times of arrival from Eq. 

(4.15). We will discuss them in the rest of this section. 

(ja) Electron Density 

In order to have an estimate of the effect of a denser solar 
wind model on the velocities determined by our method we computed 
times of arrival using a model about eight times denser than that of 
Fig. 4.l6, and fitted a few events. We selected the events at the extremes 
of the unambiguous distribution D f pig. 4.23. The results are shown 
in Table 4.17- 
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Table 4.17 

Velocities 3 Obtained from Different Models 


Event No . 

Model Fig. 4.l6 

8 x Model Fig 

8 

0. 18 

0.38 

36 

0.18 

0.38 

11 

0.19 

0.38 

23 

0.19 

0.38 

64 

0.19 

0.39 

28 

0.20 

0. 4l 

3 

0.43 

0.80 

32 

0.49 

0.70 


As expected the velocities corresponding to the denser model 
are higher. Except for event No. 32 the ratio of velocities is about 
two. The distribution is shifted up to range between 0.38 and 0.80, 
and the average could be expected to be somewhere around 0.60. These 
values appear too high. Also it is suspicious not to obtain velocities 
lower than 0.38. At the higher velocities the errors in fitting are 
larger because the time delays are smaller than the time resolution 
of the radiometer. 

(b_) Trajectory 

In our treatment we assumed that the trajectory of the particles 
lies in the plane of the ecliptic. A three dimensional trajectory 
needs to take into account not only the heliographic latitude of the 
associated optical flare but also the position of the Earth with respect 
to the solar equator. The geometry is shown in Fig. 4.26 where: 
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QP = A = heliographic longitude of a flare at F, 

PF = 6 = heliographic latitude of a flare at F, 

QT = B = heliographic latitude of the solar disk center and 

FT = \p = angular distance "between the radius at the center of 

the flare and the Sun-Earth line. 

Approximating angles "by integer degrees , the amplitude of B q is 

never larger than 7° • From Table A. 2 we note that the absolute values 

of 6 are not larger than 25° • Because 6 and B^ are not too large we have 

made an approximation that will allow us to estimate the effect on 

arrival times for trajectories out-of-the-ecliptic. If* the exciter 

particles are ejected radially from F and if the Earth is on the plane 

of the solar equator then it is simple to calculate the time delays 

of the bursts observed at the Earth. The approximation we have made is 

to assume that ip equals F'Q rather than FT, that is, we have assumed 

the Earth on the solar equator and we have increased the declination PF 

by FF' = B . Let us estimate how good the approximation is. When A = 0 

the situation corresponds to the exact case. When A = 90° we have 

calculated, for 6 = +25° and B = -7° » that FT = 93° and F'Q = 90° • 

o 

For our purposes the difference is negligible, especially when we recall 
that the uncertainties in the flare position is of the same order or 
larger. 

Table 4.l8 shows delays based on the electron density model of 
Fig. k.l6 and assuming a radial path for the exciter particles. We 
recall that the time resolution of our radiometer is of 0.15 min. Let 
us take as an example a burst originating at a flare located at A = 0, 


Time Delays in Minutes Predicted at the Earth and Between the 
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6 = +25°, and assume that B q = -7°. Using our approximation we can 
consider an equivalent 6 equal to 32°. If the particles are traveling 
with a velocity of 0.3^c we observe a delay of approximately 15.89 min. 
between the arrivals of the 3-5- and 0.05-MHz bursts, assuming second 
harmonic emission. However, the original assumption that the trajectory 
lies in the plane of the ecliptic is equivalent to take 6=0. For 
this condition the delay closest to 15.89 min. in our table is of 
15.81 min. and it is associated with a velocity of approximately 0.30c 
only. Then according to the above conditions we have underestimated a 
velocity by about 12%. It may be expected that for a spiral trajectory 
the error is of the same order. 

In Fig. k.27 we have plotted 8 versus |6 - B q | for the unambiguous 
cases in an attempt to detect the effect of neglecting the third 
dimension. There is no obvious correlation that would indicate the 
existence of such an effect. 

( c.) Propagation of Radio Waves 

The neglect of the effect of the solar plasma on the propagation 
of the electromagnetic waves results in an underestimation of the 
velocities. The correction for the bending of the rays due to 
refraction would make the path longer than we assumed (see Section 
k.2.2) and the correction for group retardation would make the time 
delays shorter. It is difficult to evaluate these errors. From Fig. 15 
in the paper of Wild et al . , (1959), we can grossly estimate that an 
apparent velocity of 0.35c should be about O.kOc when this correction 
is made. The error is approximately of 12% in this velocity range. 



158 - 



I6 # 20° 24° 28° 32° 

|S-bJ 

.NE EFFECTS ON THE VELOCITIES OF 
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The estimated errors given above should be taken with caution 
■ and they could be considered perhaps as lower limits. They can be 

I tentatively used in discussing coronal models. For example, let us 

( 

take the event No. 32 that has the following parameters (Table A. 2): • 

| ^ = -2°, 6 = +18° , B q = -7°. From Table 4.17 we see that if the true 

distribution actually corresponds to the denser model then 8 * 0.70. 

! If we further assume that we have underestimated the velocity by 2k% 

\ 

because of effects (b) and (c) then the corrected value would rise to 
8 s; 0.87. This value is high but not impossible. Wild et al. (1959) 
reported bursts with velocities between 0.8c and 0.9c. However, the 
assumed correction factor of 24$ would exclude any coronal model that 
would give '8 >^0.8l because then the corrected velocities would exceed 
I the speed of light in vacuum. 

4. 3 Spectra 

4.3.1 Introduction . As was shown in Section 2.1.5, the 
output of the OGO-V radiometer gives a quantity RT that is proportional 
to the incoming flux density, see Eq. (2.8). 

The concept of the power spectrum of Type III solar bursts is 
different from that of galactic or extragalactic radio sources where, 
due to the fact that the source size scale is small compared to their 
distances from the Earth, it is assumed that the power emitted at the 
, different frequencies come from a common region. 

For solar Type III sources the situation is different for 


several, reasons.: 
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(a) The emission at different frequencies originates at different 
levels of the corona, 

(b ) The physical conditions of the ambient coronal plasma and 

of the exciter streams are different at different distances from the 
Sion, 

(c) Radiation at a particular frequency may originate from two 
different levels , from one level as fundamental and from an outer 

level as second harmonic. It is necessary then to deal with two spectra, 

(d) The sources at very low frequencies may originate close to 
the Earth , 

(e) The directivity pattern of the electromagnetic radiation 
at the source, for both fundamental and second harmonic, will affect 
the observed flux density depending on the shape and orientation of the 
trajectory of the exciter particles. 

(f) The effect of the coronal plasma on the propagation of the 
electromagnetic waves to the Earth can modify drastically the received 
radiation from its origin, 

(g) In general the position of the source changes during the 
lifetime of the burst. 

The factors listed above added to the complexity of our data 
makes it very difficult to interpret the measured intensities especially 
since little is known about these factors. The identification of 
harmonic and fundamental emission is possible. However the presence 
of both modes may require very strong events , but this condition then 
makes the separation of components more difficult. In fact, we have 
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seen that the time profiles of strong events are highly complex and 
that the components usually merge into one big burst. 

The anisotropy of the plasma waves and the directivity of 
electromagnetic waves at the source of Type III bursts has been 
theoretically investigated, for example by Smith (19^9) > t>y Melrose 
(1970), and especially by Smerd, Wild, Sheridan (1962), by Zheleznyakov 
and Zaitsev ( 1970a, b) and by Yip (l970a,b). 

The patterns given by Zheleznyakov and Zaitsev for the 
fundamental and second harmonic emission are quite directive . There 
is some loss of directivity in the second harmonic when the reflected 
emission is taken into consideration, but in general the orientation 
of the radiation diagram with respect to the Earth should influence the 
observed intensity, particularly in the fundamental mode. 

All the theories however have been developed for a range of 

frequencies extending down only to about 20 MHz, and it is not 

3 

obvious that they are still valid at frequencies 10 to 10 times lower. 

The effects of the refractive, reflecting and absorbing 
properties of the corona on the propagation of type III bursts is not 
fully known. These effects and the directivity of the emission are 
related to the physical properties of the exciter streamer and of the 
local coronal plasma in which the burst originates. 

In order to search for directivity effects we grouped the OGO-V 
kilometric wavelength events according to the lowest frequency (f^) 
which they were observed. We defined four groups corresponding to 
f equal to 0.35> 0.2, 0.1, and 0.05 MHz. Within each group we plotted 

J_l 


-162- 


the distribution of the bursts in angular distance (9_) of the associated 

r 

flare from the Sim's central meridian. The results are shown in 
Fig. it. 28. They indicate that the eastern tail of the distribution is 
displaced to the west as the lowest frequency decreases. This can be 
attributed to directivity effects at the emission and during propagation. 

The effect of the OGO-V radiometer system on the measured flux 
densities can be attributed to the antenna radiation pattern and to 
the receiver response profile. The receiver response and mismatch 
at each of the observing frequencies is taken into account by the 
preflight calibration curves and by the correction factors computed 
in Section 2.1.5; the antenna system can be assumed to be that of a 
short dipole, as was discussed in Chapter II. In Section 2 . 2.2 we saw 
that the direction of the nulls of the antenna rotates slowly in a 
plane that is approximately normal to the Sun-Earth line. If a burst 
occurs whenthe axis of the malls is almost contained in the ecliptic 
then the source of the bursts at about 0.05 MHz or below may be close 
to the direction of a null. This depends on the trajectory of the 
exciter particles, on the mode of the emission and on the electron 
density distribution. Figure 4.29 is an idealized picture of two 
trajectories of exciter particles viewed from the north ecliptic pole. 

The circles shown are the fundamental plasma level corresponding to 

p 

Model D in Table 4.9. The figure shows also a cos 9 dipole pattern 
centered at the Earth when the nulls axis is passing through the 
ecliptic. The direction of 0.5 and 0.1 of the maximum gain are also 
indicated. The trajectories SA and SB correspond to a west-limb 
and center-limb flare, respectively. 




0.35 MHz 



FREQUENCIES OF OBSERVATION. 






0.025 MHz 
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FIG. 4.29 TWO TRAJECTORIES OF EXCITER PARTICLES AND A cos 2 9 DIPOLE PATTERN 


WHEN THE NULLS DIRECTION LIES ON THE ECLIPTIC. 
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The position of the antenna with respect to the planes (F) and 
(Q) in Fig. 2.9 is known at any time from the values of the angles $ TEU 
and SEU, listed as 9 and 0, respectively, in the OR 3/A predictions. 
When 0 = 0° or l80° the direction of the nulls is contained in the 
ecliptic; when 9=0° the direction is normal to that plane. The values 
of 9 and 0 at the times of the events analyzed in this work are listed 
in Table A.l of Appendix 1. 

**•3.2 Power Spectra . The fact that the source of a Type III 
burst moves outward through the corona introduces the concept of a 
composite spectrum. When excited a uniform and homogeneous plasma emits 
a spectrum of radio waves at and above the corresponding plasma frequency. 
This is called by us local spectrum . Assume that the passage of the 
exciter particles through the corona produces a local spectrum at 
each level. Then as the exciter particles travel outwards they produce 
a continuous succession of local spectra. Because of the negative 
electron density gradient in the corona the successive local spectra 
are shifted toward lower frequencies. At any instant there will be 
many local spectra excited at a different stage of development because 
of the finite velocity of the exciter particles. The local spectra are 
also different in form because of the different local conditions, 
expecially the electron-ion collision frequency. The incident spectrum 
at the observer at a given instant is formed by the superposition of 
the incoherent radiation from a continuous sequence of local spectra. 

This superposition of local spectra is a function of the exciter 
particles trajectory orientation with respect to the observer because 
of the different distances from the source to the observer and because 
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of propagation effects, such as scattering and refraction. The main 
fact known about the local spectrum is that the radiation at the 
fundamental and the second harmonic of the plasma frequency are 
prominent. In the earliest attempt to study theoretically the local 
spectrum Jaeger and Westfold (19U9) considered the problem of transients 
in an ionized medium. Their treatment did not consider the second 
harmonic. Modern theories have been developed that start by investi- 
gating the spectrum of the emitted plasma waves, Smith (1969), Melrose, 
(1970), Zheleznyakov and Zaitsev (l970a,b), Yip (l970a,b). 

In connection with the above discussion it is interesting to 
mention the so-called "instantaneous" bandwidth "which is defined as the 
instantaneous frequency extent at any given time," Kundu (1965 , p. 280). 
The instantaneous bandwidth is measured directly on dynamic spectro- 
grams. It lies between 10 and 100 MHz in the metric range, Wild, (1959), 
Elgaroy and Rodberg ( 19 63). Hughes and Harkness (1963) measured 100 MHz 
bandwidth at 100 MHz; using the same definition we can estimate from 
the event shown by Graedel (1969) in his Fig. 7 that in the range between 
2 and k MHz the instantaneous bandwidth is about 2 or 3 MHz. If we 
call Af the instantaneous frequency bandwidth these numbers suggest a 
relationship of the form: 


Af 

f 


* 1 


The 0G0-V radiometer obtains only limited samples of the incident 
spectrum at the 8 channel frequencies, each with a 10 KHz band, whose 
output levels are only recorded sequentially every 9.2 seconds. From 



these samples a plot is made of peak: power level against frequency. 

These will he referred to hereafter simply as power spectra of Type III 
hursts . 

When dealing with lengths at the different plasma frequency levels , 
especially if the frequencies are very low, it is convenient to have 
a reference length scale. If we assume that the electron density 
varies as r then it can he shown that the distance Ar between levels 
with plasma frequencies f and Pf t where P < 1, is 

Ar = — (l - P)r 
P 

Adopting p = 2.38, see Table b.9 , and a 10 % drop in frequency we 

obtain: 

Ar = 0.08 r 

where Ar and r are in solar radii. 

Perhaps a quantity more significant for the interpretation of 
our results is the distance Ar between plasma frequency levels 
corresponding to the extremes of the receiver bandpass (10 KHz). Using 
model D in Table 1+.9 it can be shown, for the frequency range of the 
0G0-V observations, that: 

Ar = 2.8 x 10 _i * r 2 ’ 19 , 

where Ar and r are in solar radii. We will call Ar the "frequency 
length-scale." Table U.19 contains the numerical results. 
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Table 4.19 

Frequency Length-Scale for Model D 
and for the OGO-V Radiometer 


f (MHz) 

3.5 

1.8 

0.9 

0 . 6 

0.35 

0.2 

0.1 

Ar (R ) 
0 

* 0.014 

0.055 

0.19 

0.39 

2.0 

2.8 

9-9 


The energy emitted at any one level in the corona is a function 
of the size and shape of the packet of exciters and of the time elapsed 
since their passage through the level. We will make the important 
assumption that the extension of the source at any frequency is very small, 
compared to the frequency length-scale. Preliminary report of the 
IMP-VI satellite data indicates that the apparent size of the sources 
at low frequencies are large if they are close to the ecliptic. 

Haddock (1971). 

In spite of the fact that the fundamental and second harmonic 
emissions are not separable in our spectra the discussion that follows 
will only be for the fundamental; a similar treatment can be made for 
the second harmonic. 

We can express the flux density by the functional dependence: 

S = S(f,t) . (4.l6) 

For a fixed frequency this equation represents the time profile 
of the burst at that frequency. To investigate the amplitude of time 
profiles we can plot quantities of interest versus frequency. A one 
such plot is the maximum flux density versus frequency. To correct 
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for the changing distance of the source from the Sun we can reduce 
the observed flux densities to a common distance such as 1 astronomical 
unit. To compensate for the inverse distance-squared reduction, the 
flux density of radiation observed at a frequency f and emitted at a 
distance x from the Earth should be multiplied by a factor: 


q. 



> 


u r 

where £ is the Sun-Earth distance. See Fig. 4.6. This factor ignores 
the effects of the corona on the direction of propagation of the 
radio waves and also depends on the electron density model assumed, 
q can be computed from the Eqs . (A. it) and (A. IT). Table 4.20 shows 
values of q for different conditions. It is observed that the values 
on some columns reverse their trend as the frequency decreases. This 
happens when the sources start receding from the Earth. 


Table 4.20 

Values of the Correction Factor q for a Spiral Trajectory 



Model C 
Fundamental 

Second 

Harmonic 

Model 

Fundamental 

D 

Second 

Harmonic 

X 9 F 

(MHz)\ 

*/2 

0 

jt/2 

0 

*/2 

0 

*/2 

0 

3.5 

1.00 

0.96 

1.00 

0.94 

1.00 

0.94 

.99 

.89 

1.8 

1.00 

0.94 

0.99 

0.90 

• 99 

.90 

.98 

.82 

0.9 

0.99 

0.90 

0.98 

0.83 

• 98 

.82 

.94 

• 71 

0.6 

0.99 

0.86 

0.97 

0.76 

.96 

.76 

• 89 

.63 

0.35 

0.97 

0.79 

0.92 

0.66 

• 92 

. 66 

• 75 

.52 

0.2 

0.93 

0.69 

0.80 

0.54 

.80 

.54 

.44 

.60 

0.1 

0.80 

0.54 

0.45 

0.60 

.44 

.60 

.00 

2.03 

0.05 

0.45 

0.60 

0.00 

2.01 

.00 

2.03 

3.08 

7.63 
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^•3-3 Energy Spectra . In deriving equation (2.8) that relates 
the measured RT product and incident flux density we assumed that the 
radiation was randomly polarized. The assumption was implicit in the 
definition of the flux density S f that refers only to one state of 
polarization. This assumption is not always valid in the case of 
Type III hursts. Komesaroff (1958) observing in the range 1+0 to 2*+0 MHz 
found that half of the 500 bursts exhibited the characteristics of 
fairly strong circular polarization. Takakura ( 196U ) states that "most 
of Type III bursts at 200 Mc/s seem to be unpolarized or weakly 
circularly polarized." Kai (1970) observing at 80 MHz found that the 
bursts showed weak circular polarization or were unpolarized while 
Harvey and McNarry (1970) examples illustrates the full range of variation 
in ellipticity in the polarization of bursts observed at 7U MHz. In 
view of this we will assume that the radiation of the bursts under 
study was unpolarized, especially since the linear dipole will receive 
one-half of the incident flux density for circularly polarized signals. 

If we integrate S(f,t) over the duration A of a burst we obtain 
the energy density received per unit bandwidth centered at the frequency 
f. Assuming that the radiation is emitted into it steradiaft's and that 
it comes from a small source, the total energy emitted by the burst in 
one polarization and per unit bandwidth is: 

E(f) = ™ 2 (f\/ A(f) S(f,t) dt , joule m" 2 Hz" 1 , (1+.17) 


where x is the distance from the source to the Earth. 



Finally the total energy emitted by the radio event in the 


range of the frequencies of observation is 

E t = f 2 E(f) df , joule m" 2 . (fc.l8) 

1 

This result applies to the spectrum of either the fundamental or 

second harmonic emission. The addition of their contributions gives the 

total of energy of the burst in one state of polarization. This is in 

theory because in practice we did not separate them. 

4 . 3 .^ The OGO-V Observations . As we have seen there are many 

factors that contribute to the spectrum of the observed flux densities 

and any interpretation has to be considered very provisional. Before 

we realized the prevalence ' of the second harmonic emission below 1 MHz 

we selected those events having a strong component that was recorded on 

the eight frequencies. We assumed that these components corresponded to 

the same burst and to the fundamental. Their peak-flux density was 

computed from Eq. ( 2 . 8 a). The correction factor q for inverse square of 

distance effect was based on a coronal model about 20% less dense than 

that of Fig. I+.16. As a consequence q was close to ■unity in most cases. 

In the following flux densities the contribution from the 

background noise level has been subtracted. Considering the shape of 

these peak-power spectra we can distinguish two basic types: single 

and double peaked. They are illustrated in Fig. k. 30 . Usually the 

-13 -12 

amplitude of the maximum fluctuates between 3 x 10 and 3 x 10 

-2 -1 

w m Hz These measurements are about seven orders of magnitude 


ifiO 
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large r than the ones obtained from ground-based observations at 
frequencies above 10 MHz. This could be partly due to the closeness 
of the burst sources to the Earth but from Table 4.20 we see that the 
correction, even at 0.2 MHz emitted as second harmonic, is small and 
practically independent of the orientation of the trajectory of the 
exciter particles. However for the combination of model D, second 

harmonic and e_ = tt/ 2 the 0.1-MHz flux density should be much affected. 

. r 

In Fig. it. 30 we include the orientation © F of the trajectory. We 
observed that the flux density at 0.1 MHz of event No. 2 6 is very large 
in spite of the fact that it originated 18° East of the Sun's central 

meridian. 

An unequivocal identification of modes and possible correlations 
of spectral shapes with the position of the associated flares would 
help to learn more about double peaked spectra. In particular the 
trough could be produced by the onset of the second harmonic with 
increasing intensity. 

We me asured the energy spectrum E(f) for eight events (Nos. 13* 
16 , 24, 31, 32, 33, 34, and 46) which were selected because of the 
simplicity of their time profiles. To obtain E(f) we measured the 
area under the time profile at each frequency without regard to mode 
separation . The results are shown in Fig. 4.31 where the symbol x 
in our 0G0-V data indicates the average over the different events. The 
error bars represent the extreme values measured. The number of 
events taken at 3.5, 1.8, 0.9, 0.6, 0.35» 0.2 and 0.1 MHz were 8, 8, 

8, 8, 3, 3 and 2, respectively. At 0.1 MHz the bar represents one case 
and the arrow the lower bound of another case. We notice that the 
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FIG. 4.31 AVERAGED ENERGY SPECTRUM. ALL VALUES EXCEPT WILD'S AND 
OGO-V's WERE DEDUCED FROM THE REFERENCES. 
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2 

extremes can fluctuate by a factor of 10 , however, they indicate a 
definite trend with frequency. In particular the average follows a 
smooth curve that seems to peak about 0.3 MHz. In the same figure 
we show similar data obtained by Wild ( 1950b) between 70 and 130 MHz. 

He found that a curve of the form f n well fitted his observations 
and that the best fit is obtained with n = 3.25* Plausible values of 
n lie between 2.8 and 3.6. We found that the equation of his best 
fit curve is : 

E(f) = 10 -11 * 65 f -3,25 

where f is in MHz. 

The exponent in Wild's functional form was misprinted on page 
279 of Kundu's book ( 1965 ) as -3-50 instead of -3.25, as given by 
Wild. This error has been propagated by Melrose (1970) and by 
Graedel ( 1969 ). Graedel also misinterprets Wild's equation as referring 
to flux density rather than to energy. In Fig. k.31 we plotted Wild's 
curve for n = 2.8 and n = 3.6. We found that the agreement for n = 3.6 
with the 0G0-V values is good but n = 3.8 would be better. The rest 
of the data included in Fig. I*. 31 was deduced from the references 
indicated and it should be taken cautiously. The increase in flux 
density with decreasing frequency observed by Wild was also noticed 
by Slysh ( 1967 b) who compared a burst at O .985 MHz with the associated 
Type III burst detected from ground. Values as low or lower than that 
deduced from Slysh 's publication can be expected since it is a matter 
of the weakness of the event and the capability of the radiometer to 
detect it. We wish to remark again that the 0G0-V data corresponds to 
a selected group of events characterized by high intesnity. 
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In order to estimate the total energy in the spectrum defined 
in the low frequency range by the OGO-V average values , we fitted it 
to a function of the form: 


E(f) = E m (r) 

where n and p are constants, 

= energy of the peak of the spectrum, and 
fjyj = frequency at which occurs. 

For large values of f the function behaves as f n , as required 
by Wild's observations. Our observations of the eight events plotted 
in Fig. 4 . 31 suggest the following approximate values: 

Em = 3 x 10 10 , joule m ^ Hz ^ 

f M = 0.3 , MHz . 

We found that with p = 1 and n = 3.8 to 4.2, the function 
Eq. (4.19) gives a good fit to the data. The curves are shown in 
Fig. 4.32. For the purpose of computing the total energy, both curves 
are satisfactory. 

The definite integration of Eq. (4.l8) using (4.19) can be 
achieved by the following change of variable: 




(4.19) 


r 



5 


which transforms the integrand into the form of an incomplete gamma 


function. The difference between two incomplete gamma functions gives 
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FIG. 4.32 FIT OF EQUATION 4.19 TO AVERAGED ENERGY SPECTRUM. 
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the result. We found that the spectral energy density received at the 
Earth in the band between 0.1 and 300 MHz, using n -k.2, is of 
1.5 x 10-^ joule m - 2. 

We can deduce the total energy produced by this "average" radio 
event, assuming that the source is at 1 A.U. , that it is small and that 

x 8 

it radiates into TT stereoradian. The result is E = 7 x 10 joule. 

This value could be used to estimate the number of particles involved if 
we knew the kind of particles, the efficiency of the transformation of 
the kinetic energy of the exciter into plasma waves and the efficiency 
of the coupling between plasma and electromagnetic waves. The nature 
of the exciters is not known and electrons and protons have been 
proposed. Since the average speed of the particles is 0.32c this 
corresponds to 30-KeV electrons or 5^-MeV protons. We have not found 
explicit information in the literature about the loss of kinetic energy 
to plasma waves. 

Little is known about the effectiveness of the coupling of 

plasma waves into electromagnetic waves; to complicate the problem, 

it is necessary to differentiate between coherent and noncoherent 

processes. The results also depend on the mechanism invoked and on 

whether the emission is at the fundamental or the second harmonic 

frequency. Ginzburg and Zheleznyakov (1958a, 1958b) have derived an 

approximate expression for the efficiency Q of the conversion of 

incoherent plasma waves into electromagnetic waves . An upper value 

”8 T 

of Q is about 3 x 10 , while a mean value is < 10 . This subject 

is complex and uncertain and we consider any attempt we make to estimate 
the number of particles involved almost meaningless. This situation 
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is further aggravated by the fact that we do not know how applicable 

are the current theories in the kiloinetric range of wavelength. 

To estimate the number of electrons involved Melrose (1970) 

computed the total energy E^ from Wild's (1950b) data, in a way similar 

to ours. His calculations assumed however a low- frequency cutoff near 

11 12 

35 MHz. He arrived at a total energy of 10 - 10 joules. (The 

13 l4 

actual values he gives are 10 - 10 joules due to an error in the 

exponent of f in his equation 8, Paper II. The error originates from 
a misprint in Kundu's book ( 1965 ) where the exponent appears as 
-3-5 when it should have been -3.25.) Notice that the total energy 
derived by us is about 10 times larger. We present in Table 4.21 
the total number N of particles estimated by different authors. The 
subscript e and p refer to electrons and protons, respectively. 


Table 4.21 

Number of Exciter Particles for Type III Bursts Estimated by Various 

Authors 

N N Energy /Particle Authors 

g P 


10 


35 


9 x 10 


29 


5 x 10 


25 


33 34 

10 ^ - 10 J 


is3 x 10 


33 


30 KeV 
50 MeV 

> 40 KeV 


Wild (1964,1969) 
Smith ( 1969 ) 

Smith (1970) 
Melrose (1970) 
Lin (1968) 
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Th_e 18 November 1968 Event. This event, No. 21, deserves special 
consideration because it is the only one analyzed that was observed down 
to 0.05 MHz, see Fig. 3.2. It was the radio component of a complex 
solar proton event associated with a west-limb flare. We recall that 
in our analysis of velocities we used a coronal model close to Model C, 
and that the event No. 21 is ambiguous, having the alternatives 
3 = 0.l8 radiating at the fundamental or 8 = 0.36 radiating at the 
second harmonic. (Table A. 3 Appendix No. l). We assume that the 
second harmonic is most likely because of the results shown in Fig. 1*.22. 
Figure ^ . 33 shows the peak-power spectrum of this event without the 
distance correction. From Table 1*. 20 we see that for Model C a correction 
due to the inverse square of the distance would change the 0.1-MHz 
flux density by a factor of 0.80 or 0.1+5, depending on the mode. This 
is not much considering the uncertainties and the high power level of 
the whole spectrum. We decided to present the spectrum without the 
distance correction although this correction is not negligible at 0.05 
MHz. In the same figure we include also the available ground-based 
data. The event occurred at a time, 1026.5 U.T. , in which the 
major solar radio observatories were at night. 

1 

A comprehensive report on the phenomena associated with this 
important event has been published by the World Data Center A, 

Report UAG-9 , 1970. The observations of the Weissenau Observatory 
are not included there. 

The event was observed over 6 decades between 36 GHz and 0.05 MHz, 
which makes it the most extensive frequency coverage of any radio event 
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FIG. ^-33 PEAK-POWER SPECTRUM OF THE l 8 th NOVEMBER 1968 EVENT. 
BURSTS GROUP STARTING AT 10^0 U.T. NO DISTANCE 
CORRECTION APPLIED TO OGO-V DATA. 
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so far reported. The data shown in Fig. 4.33 correspond to the group 
of hursts that started at 1040 U.T. at 3.5 MHz. Previous to this 
strong group there was a weak group starting at 1026.5 U.T. at 3.5 MHz, 
that is, at the time of the optical flare. See Fig. 3.2b and Table A. 2. 
The radio event was extremely complex in the entire frequency range. 
Urbarz (1969) reports bursts of at least Types II, III and IV in his 
spectral observations between 30 and 965 MHz. Thus, it is difficult 
to establish what part of the single- frequency ground-based observations 
used in Fig. 4.33 corresponds to the Type III group. Above about 100 MHz 
the spectrum follows the same pattern found by Castelli et_ al . , (1967) 
in other proton events . 

Because of the typical starting frequency observed in the Type III 
bursts (Malville , 1962a, 1967 ) we could say that the spectrum below 
260 MHz (and certainly below 23 MHz) belong to the Type III group. The 
shape of the spectrum between 0.05 and 260 MHz suggests a double peak 
as those shown in Fig. 4,30b . It is observed that the main peak seems 
to occur, at 0.05 MHz or below. 

Figure 4. 34 shows the energy spectrum of this event, not corrected 
for distance, which is more striking than the power spectrum because, 
instead of turning over at 0.05 MHz as in Fig. 4.31, this spectrum 
continues to rise at a rate approximately proportional to f - ^. At 0.1 
MHz the spectral energy density received at the Earth is between two 
and three orders of magnitude higher than typical strong events observed 
by 0G0-V. An estimate indicates that the total integrated energy 
involved is about 40 times larger than that in a typical strong event. 
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Theoretical Implications . If we accept the flux densities from 

the OGO-V observations we find that they exceed by several orders of 

magnitude the values obtained at higher frequencies from ground-based 

stations. This should be useful information in any theory developed 

to explain the emission mechanism of Type III bursts. The fact that 

we have selected the more intense events is interesting because, as 

Smith ( 1969 ) recognizes, "the explanation of the most intense burst 

is likely to present the most stringent test of a model." Smith 

_19 _2 -1 

considered 3-2 x 10 w m Hz as a maximum possible peak value and 

applied his model to a frequency as low as 20 MHz. The theory developed 

by Zheleznyakov and Zaitsev (1970) is not applicable to bursts with a 

-17 -2 -1 

flux density exceeding 2 x 10 w m Hz at 100 MHz and, as they 
indicate, the generation of bursts that violated this condition 
deserves a special investigation. 

-12 -2 -1 

The average peak-value observed by us is 10 w m Hz . We 

conclude that our observations may represent a severe test to theories 

on the generation of Type III bursts although the highest frequency 

observed by 0G0-V is at least one decade below the lowest frequency 

to which theories have been applied. 

Another aspect in which our observations should influence theory 

refers to the second harmonic emission. The relative abundance of one 

of the two modes in Type III bursts as a function of the strength of 

the bursts has been pointed out by Zheleznyakov and Zaitsev ( 1970 ) 

and by Melrose (1970). Our results contradict the prediction of the 

t'Oviet authors that tin* "intensity ratio of the fundamental and second 

harmonic should increase as the concentration N in the stream decreases 

s 
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and thereby as the power of the burst decreases. One would expect in 
this connection that the percentage of bursts where the fundamental 
and second harmonics are simultaneously present would also increase 
as the intensity of the burst decreases." Our results, on the other 
hand, support Melrose's (1970) model that predicts that "the more 
intense the burst the higher should be the ratio of the intensity at 
the second harmonic to that at the fundamental." None of the proposed 
emission mechanisms predict the increase of the relative content of 
the harmonic with decreasing frequency found by us and illustrated in 
Fig. 4.22, however, the second harmonic prevalence may be due to 
phenomena affecting the radiation during propagation. 

Mai vi lie (1962b) predicted a turn-over in the energy spectrum, 
but he placed it at 50 MHz. Smith (1969) accounts for the increase 
of flux density at lower frequencies by invoking a decrease in the 
width of the velocity distribution function of the exciter stream. 

4. 4 Time of Decay of Time Profiles 

4.4.1 The Observations . The time profile of a burst at a 
frequency f is defined by Eq. (4‘.l6) with f constant. The output from 
each of the channels of the OGO-V Radio Astronoiry receiver constitutes 
a time profile. Figure 3.2a shows the eight profiles for the 18 November 
1968 event. 

The shape of a time profile depends on many factors, for example 
the size, shape and energy of the bunch of exciters particles, the 
temperature and density of the local plasma, and the modifying effect 
of the receiver bandwidth. In this respect it is interesting to calculate 
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the time of transit of the exciters particles between the two plasma 
levels corresponding to the frequency extremes of the bandpass of the 
OGO-V radiometer for the different frequencies. It depends on many 
uncertain factors (velocity of exciters particles, coronal model etc.,) 
however we can obtain representative values from the empirical drift- 
rate law, Eq. (h. 2)^ adopting the parameters corresponding to "all data" 
shown in Table k.3. If Af is the receiver bandwidth of 10 kHz then 
the transit times At are given by: 

At = f _1,81t sec , 

where f is in MHz. The At calculated for the OGO-V frequencies are 
shown in Table it. 22. 


Table it. 22 

Typical Transit Times Through the Eight OGO-V Bandpasses 

f (MHz) 3.5 1.8 0.9 0.6 0.35 0.2 0.1 0.05 

At (sec) 0.10 0.3it 1.21 2.56 6.9 19.2 69 2 b6 

As was stated in Chapter II, the receiver time constant is 0.21 
seconds, but each channel is sampled every 9.2 seconds, therefore down 
to about 0.35 MHz we are limited in time resolution by the sampling 
rate. On the other hand at frequencies equal or less than 0.2 MHz 
the transit time through the bandpass prevents measurements with an 
accuracy better than At . 

One of the most important characteristics of a time profile 
is the rate of amplitude decay . We observe the decay of the 
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electromagnetic radiation that, according to the local plasma hypothesis, 
depends on the decay of the plasma oscillations set up by the passage 
of the exciter particles through the corona. 

Jaeger and Westfold (1949) developed the thoery of transients 
in homogeneous ionized medium. They concluded that the oscillations set 
up hy a pulsed excitation decay with time according to an exponential 
law 

-vt 

e 

where v, the ion-electron collision frequency (Coulomb collisions) is 
given by: 


V = 1*. 2 x 10 5 NT 3/2 sec" 1 (4.20) 

where N = electron density, m 3 and 

T = electron kinetic temperature^K. 

They proposed that Eq. (4.20) could be applied to bursts of solar 
noise produced in the solar corona. Evidence found by Williams (1948) 
and Payne-Scott (1949) suggested that the decay in intensity of the 
bursts was approximately exponential; but they did not compute 
temperatures. Further study was reported by Jaeger and Van't Verr 
(1957) and Boischot et al. (i960). These last authors used Eq. (4.20) 
in the range of 15 to 38 MHz assuming the local plasma hypothesis and 
computed temperatures in agreement with that obtained from optical 
observations. Mai vi lie (1962a) using a similar formula given by Spitzer 
obtained good agreement also in the derived temperature at 25 MHz. 
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Boischot (1967) found that the temperatures obtained from observations 
at 36, 18 and 8 MHz showed a decrease with frequency, as expected from 
pure collisional damping. The first measurements from an artificial 
earth satellite reported by Hartz (196*0 also showed this trend between 
10 and 2 MHz that agreed with Slysh's value at O.985 MHz. However, 

Hartz (1969) first pointed out that the temperatures obtained in the 
range of 9 to 0.7 MHz were "substantially below those expected for the 
solar wind." Further studies of Type III bursts at lower frequencies 
have confirmed this disagreement. Haddock and Graedel (1970), 

Alexander et_ al. (1969 ) . 

Criticisms of this determination of temperature has been raised 
by Hughes and Harkness (1963). They argue that the time decay of a burst 
is also affected by the increase in the size of the exciter packet by 
dispersion and by contamination of the Type III burst with Type V and 
with second harmonic radiation of the same Type III. If these factors 
are taken into account the deduced terpieratures would be even lower thus 
increasing the discrepancy. 

We believe that it is necessary to ascertain that the burst decay 
is in fact exponential before using Eq. (*+.20) or its equivalent. How- 
ever, to present the 0G0-V results we have measured the time x that 
it takes a burst to drop from its peak to e of that value regardless 
of whether the decay is exponential or not. Besides the factors that 
affect the time profile mentioned by Hughes and Harkness we should add 
the impossibility of resolving individual components in the extremely 
complex time profiles of the events we have analyzed. All these effects 
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become worse in our data as the frequency of observation decreases. In 
many of those cases affected by saturation and telemetry quantization 
(see Chapter II ) we were reduced to drawing an idealized curve through 
the observed profile. Because of all the reasons discussed our 
measurements of T are upper limits of the times of decay expected from 
collisional damping. This is supported by the fact that in bursts that 
exhibit saturation the decaying portion, after the saturation terminates, 
decreases at a faster rate than the ideal profile. This decay could be 
closer to the free relaxation after the exciter particles have passed 
through the appropriate plasma level. Figure 1+.35 shows the measured 
decay times. The sybmols x indicate the averages, the length of the 
error bar is 2 a centered at the average. The number in parenthesis 
indicates the number of cases considered. At 0.05 MHz we have used the 
data of the radio burst that occurred on the 6th of April, 1971 at 
0915 U.T. , observed by 0G0-V but not included in our analysis. We see 
that in spite of the crudeness of the measurements the decay times follow 
a definite increase with decreasing frequency. In Fig. b .36 we have 
gathered all the information we could find in the literature. It is 
interesting to observe that the observations fall approximately on a 
straight line, except for a few points. The points obtained by 
Boischot et_ al. (i960) that depart the most were considered perhaps 
erroneously high by Kundu (1965, p. 329). 

In this plot we can fit approximately a straight line whose equation 
is 

Tf °-97 = 10 7.75 _ 


(4.21) 


DECAY TIME 
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FiG. 4.35 DECAY TIMES OF TYPE III BURSTS OBSERVED BY OGO-V. 
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where x is in seconds and f in Hz. It is interesting to note that 
Wild (1950b), observing between 80 and 120 MHz, fitted a line 

Tf = 10 8 , 

that is, practically the same. As Bracewell pointed out in the dis- 
cussion of the paper by Alvarez and Haddock (1971), this means that 

Q 

it requires 10 cycles of the radiated waves for damping to occur, 

independently of frequency or plasma density. 

In order to study the experimental results, we have obtained the 

decay times x^ due to Coulomb ion-electron collisional damping derived 

from some theoretical solar wind models. X is defined bv 

c J 

T c - £ <k.22) 

and v is given by Eq. (4.20). Figure It. 37 shows x for two solar 

c 

wind models and one solar streamer model; it also shows the measured x 

that, as we said, constitutes at best upper limits. The increase in 

the disagreement with diminishing frequency is obvious and at 1 MHz the 

difference is of one order of magnitude. It is interesting to observe 

that the variation of x^ with frequency is similar in the three 

theoretical models in spite of the fact that Pneuman and Kopp's streamer 

(core) is very thin. We realize that actual Type III bursts may not 

occur in media given by these models. However, at the lowest frequencies, 

where the discrepancy is larger the bursts may occur in the solar wind 

rather than in well-defined streamers. 

At 1 A.U. Whang et al. model gives x = 18.U hours and Hartle 

c 

and Sturrock's give an even longer time. This large disagreement could 
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be explained by invoking additional damping to the ion-electron collisions 
contribution (we will see shortly that this is what the solar wind 
theoreticians have assumed.) We can think for example, of the dissi- 
pative effect of turbulence in the local plasma of either a permanent 
nature or a transient character and perhaps produced by the passage 
of the exciter particles. The time profile may be sensitive to the 
efficiency of the coupling between plasma waves and electromagnetic 
waves. This efficiency could perhaps be altered by the passage of the 
exciters particles. Figure 4.37 indicates that whatever the source 
of the additional damping, it increases with decreasing frequency. It 
would be important to measure the times of decay with better time 
resolution in well identified fundamental and second harmonic emission. 

4. 4.2 The Temperature of the Corona . The problem of the 
determination of the coronal temperature has proven to be more complex 
than that of the coronal density. Like density, the temperature has 
only been measured near the Sun and near the Earth. It has been well 
established that the temperature near the base of the corona ranges 
between 7*5 x 10^ and 4.5 x 10^ K, and the values are evidently 
influenced by the method used for the determination, Billings ( 1966), 
Zirin (1966), Rosh (1967). 

The temperature near the Earth has been measured from 
artificial satellites and space probes. The information from these 
experiments has brought the necessity to modify some concepts that were 
hitherto applied to the solar wind. For example, the solar wind cannot 
be considered any more as a one-component fluid rotating and expanding 
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quietly but rather as a turbulent multiple- component fluid having 
important discontinuities and being continuously crossed by plasma 
waves . 

The temperature of the protons was measured first. Due to the 

interplanetary magnetic field their random motions are not isotropic 

and it is necessary to speak of "transverse" (T ) and "longitudinal" 

(T ) temperatures; the so-called thermal anisotropy. Its magnitude 

is expressed by the ratio T /T^ that is found to be near 2. Usually 

the proton temperature is given as a weighted mean T = 1/3 (T + 2T ) . 

1 

The generally accepted value for the proton temperature in the "quiet" 

... i 4 5 

solar wind is 4 x 10 K. Values as high as 8 x lCr K have been obtained 

at times when the solar wind velocity was above its quiet value, 

Hundhausen (1968). The important point is that the assumption of 

isotropic Maxwellian distributions is not valid for the solar wind 

protons, Hundhausen ( 1968 , 1970).. 

The problems involved in the observation of thermal electrons 

have been solved only recently. The electrical neutrality of the solar 

wind and the equality of flow speeds of electrons and protons are well 

established. The thermal anisotropy of the electrons is small, and 

much less than that of protons . Here again , the velocity distribution 

function of the electrons is not Maxwellian and seems to fit a bi- 

Maxwellian distribution except for a high-velocity tail that is usually 

present. The accepted mean electron temperature of the "quiet" solar 

wind is between 10 5 and 1.5 x 10 5 K, Hundhausen (1970). Variations 
4 5 

between 7 x 10 and 5 x 10 K have been observed in a two-month period^ 
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Kavanagh et al. (1970). Usually the temperature of the electrons is 
about four times higher than that of the protons. 

The theoretical treatment of the solar wind started with one- 
fluid models, for example. Noble and Scarf (1963) and Whang et al. 

(1966), and continued with two-fluid models, Hartle and Sturrock (1968), 
Cuperman and Harten (1971)* and Wolff et_ al . (1971). The last three 
predict, to a different extent, electron temperatures higher than 
proton temperatures . 

Because of these circumstances it seems logical to expect that 
the use of Eq. (k.2l) to determine electron temperatures from radio 
observations at low frequencies would not give accurate values because 
the derivation of Eq. (1+.20) assumes a Maxwellian distribution of 
velocities of electrons and colder protons. There are other arguments 
against this formula. For example, the ratio of the mean collision 
time for protons (or electrons) to the scale time for the coronal 
expansion, although very small near the Sun, becomes unity at about 
0.13 A. U. (28 Rq). The ratio continues to increase with distance and 
beyond approximately 0.13 A.U. the plasma is essentially collisionless, 

Hundhausen (1968). Thus near 1 A.U. the electrons make about one 

. 5 -1 5 

collision per 10 sec, for N « 5 cm and T = 10 K, while it takes the 

solar wind 3 x 10 5 sec to travel from about 0.1 to 1 A.U. The 

importance of the plasma -wave-particle interaction in the physical state 

of the coronal plasma has been discussed by solar wind theoreticians. 

This interaction produces plasma turbulence that in turn "provides a very 

local effective dissipation mechanism because the wave-particles 

scattering leads to an anomalous (low) electrical conductivity," 
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Scarf (1970). This may be the way to increase the Coulomb collision 
frequency. Our results support this suggestion. As can be seen in 
Fig. 4.37 the predicted decay times are longer than the observed ones, 
that is, the predicted Coulomb collision frequencies are too small. 

This effect was considered in the two-fluid model of Cuperman and 
Harten (1971), who introduced an enhanced noncollisional coupling 
between the proton and the electrons to bring the theory into agreement 
with the observations near the Earth. This was expressed by taking, 
ad hoc, an effective collision frequency 30 times larger than the 
classical Coulomb collision frequency. Nishida (1969) had earlier used 
a factor of 10. The latest solar wind model, by Wolff et_ al . (1971), 
introduces a distance dependent coefficient to increase the collisional 
energy exchange between electrons and protons. Figure 4.38 shows 
Hartle and Sturrock's temperature model based on collisional coupling 
and that of Wolffet al . including enhanced noncollisional coupling. 
Observations are also included; the temperature near the Sun was 
arbitrarily placed at r = 1. 

In Fig. 4.37 we can draw approximately a straight line through 

the low- frequency portion of the curve derived from the Hartle and 

12.71 -1.58 

Sturrock's model. Its equation is T =.10 f , where x is in 

C 

sec and f in Hz. If we call x the observed decay time given by 

Eq. (4.21) then we can write the discrepancy in the following functional 

form: 


T 


10 5.00 f -0.63 


(4.23) 
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Using the local plasma hypothesis, Eq. (^.T)» this equation 
can he rewritten as : 


Is. = xo 1 *- 140 H- 0 - 315 y 0 - 63 




where N is in cm 7 and j = 1,2 for fundamental and harmonic, 
respectively. The distance dependence is determined by the electron 
density model chosen. Using for N the fitting given in Table h.6 for 
Hartle and Sturrock's model we can write explicitly: 


= 3.3 r °' 71 , for j = 1 , (U.25a) 

and 

:j— = 2.1 r'°’ 71 , for j = 2 . (U.25b) 

In their model Wolff et al (l9Tl) increase the ion-electron 
energy exchange by assuming that the electrom mean free path is 
proportional to the reciprocal of the distance raised to an arbitrary 
power. This power was chosen so as to give the best agreement between 
the theory and the observations at the Earth, particularly of the 
proton temperature . In terms of our notation their law can be written 
as : 


Is. = r “°-T3 

T 


(It. 26 ) 


Considering the uncertainties in the measurements of temperatures 
and in our determination of decay times we believe that the agreement 
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between Eqs . (4.25) and (4.26) is good. The comparison of Equation 
(4.25) based on Hartle and Sturrock's model with Eq. (4.26) based on 
Wolff et_ al . is valid because Wolff et_ al. reduces practically to 
Hartle and Sturrock's when no corrections are made. 

We observe that the second harmonic gives better agreement. 

To summarize Section 4.4 we can say that we have shown in a 
conclusive way, the failure of the pure collisional damping hypothesis 
to explain the decay times of Type III bursts. Our results are consistent 
with those found by solar wind theoreticians independently and by 
completely different means. Both results complement each other because 
one is valid near the Sun and at the Earth while ours is valid in a 
continuous range of distance in between. 

4.5 The Association with Other Solar Phenomena 

The Type III solar bursts are usually associated, to a different 
degree, with other manifestations of the solar activity. We will 
briefly examine the association of the bursts analyzed in this work 
with optical flares (Ha), x-ray flares, solar cosmic rays and radio 
bursts observed from ground-based stations. Other correlations with 
terrestrial effects, i.e., ionospheric disturbances, geomagnetic storms, 
etc. , will not be discussed. 

**-5-l Optical Flares (Ha). The association of Type III bursts 
with Ha flares has special importance because, except for appropriate 
radio observations , it is the only means known so far to establish 
the approximate location of the exciter particles origin. 




The association of enhancements of solar radio emission with 
solar flares was established by Appleton and Hey (19^6). 

Specifically , an association of flares with Type III bursts was 
pointed out by Wild, Roberts and Murray (195^0; afterwards it has 
been studied by several authors. 

The subject can be studied in at least two respects: time of 

occurrence and position of the flare on the disk. 

Temporal Association . In the study of the temporal association 
it is important to define the criterion 0 f coincidence, especially the 
phase of the optical event considered as time reference. This is 
confusing in the literature because the different authors use a variety 
of phase references and not all flares exhibit the same phases. For 
example, the following terminology occurs in connection with this 
subject: puff, surge, spray, flash phase, explosive phase, flash- 

explosive phase , etc. . . . 

In order to study the correlation of our data with solar flares 
we used the tables of Confirmed Solar Flares compiled by the ESSA 
Research Laboratories (now NOAA) and published in the Solar-Geophysical 
Data. In only two cases (events Ho. i* and 1 6 ) we used the Unconfirmed 
List. Reports of the same flare from different observatories are 
listed together in a group and average values for the group are given. 

The time reference for the Type III bursts analyzed by us is the 
start time at 3-5 MHz. In most cases we found that a flare occurred 
close to our radio events. We compared the time of start of the radio 
event with two times associated with a flare: (a) with the closest 
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time of start within the group of observations and (b ) with the average 
time of start of the group. The times in the ESSA lists are given 
to within one minute. We rounded the OGO-V data (usually measured to 
one-tenth of a minute) to the closest minute. The differences of 
time of start t - t^ are then an integral number of minutes. In 
plotting their distribution we centered the one minute intervals of 
the histograms on integer numbers. Figures U. 39 and l+.i+O show the 


distribution of the differences in time of start according to criteria 
(a) and (b), respectively. 


Of the 6k radio events studied eight were not considered in 
Fig. ^ . 39 ; five occurred behind the limb (1+3, 53 , 5 !+, 55 , 56 ), two 
( 27 , 1 + 7 ) did not have a reported flare during these hours and one ( 63 ) 
occurred 38 minutes after the closest flare. Two events ( 30 , 1+6) were 
probably associated with Ha flares already in progress when they were 
first observed. Evidence from x-ray flares supports the association 
so we decided to keep these two cases adopting the minimum value of 
the time difference. In Fig. 1+.1+0 ten events were not included: The 

eight excluded from Fig. 1+.39 Bnd two (30, 1+6) for which the average 
of time of start of the Ha flares was not available. Event No. h 
is probably wrongly associated because it is the one that departs the 
most from the general trend both in temporal (Fig. 1+.39) as well as in 
positional association (Fig. l+.l+l). in addition, it is associated with 
an unconfirmed flare. 

The statistics obtained from these figures are summarized in 
Table 1+.23. We observe that most of the bursts 


occur after the optical 



UNCONFIRMED FLARE 
(EVENTS 4 a 16) 
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CRITERION (a). 


-20k- 
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flare and that a large percentage occur within 1*.5 minutes of the 
optical event. The use of criterion (a) or (b) does not affect 
appreciably the results. 


Table k.23 

Temporal Association of the Type III Bursts Observed 
by OGO-V with Ha Flares 

(Number of Occurrence Are Expressed in Percentage) 


Criterion 


a 

b 


*3.5 ~ Ha (min> 

-8 to 0 0 to 2$ -2 to 0 0 to +2 

31.3 68.7 2k. 1 38. k 

36.1 63-9 20.it 28.7 


-it. 5 to +it.5 


83.0 

76.0 


For comparison we gathered in Table it.2h the results of other 
observers. We should stress the fact that our results refer to very 
intense radio events. The data of at least three references relate 
also to major radio events (Swarup et al . , Boorman et_ al . , Riihimaa). 

Because of the different criteria used for coincidence it is 
difficult to compare our results with those of other investigators. 

Malvi lie ' s criterion of Start ±2 min seems to be the closest to ours. 

He found that 50$ of the bursts were associated with flares. Our 
results for the condition Start -2 min (Table h.23 columns i+ and 5) give 
62.5$ and U9.1$, depending on the criterion used. 

The observed proportion of bursts preceding or following the 
flare varies for different authors. From Malville's figure 5 (1962c) 
we see that k3% of the bursts preceded the onset of the flash phase 
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at most by 1.5 minutes while 37$ followed it within 1.5 minutes. This 
agrees with Grae del (1970) who found that on the average a burst 
preceded the reported flare maximum by 1 . 1+5 minutes. Our data indicate, 
that on the average a burst follows a flare (column 3 of Table 4.23, 
Figs. 4.39 and 4.40). This agrees with Swarup et_ al . (i960) who found 
that most of the bursts occur after the flare start, and also with 
Riihimaa (1963) who found that 7 out of 11 radio bursts occurred 
between the onset and maximum of the flare. 

Positional Association ■ Before solar bursts were classified into 
different types Hey, Persons and Phillips (1948) found that radio bursts 
at 73 MHz had a tendency to be associated with flares occurring on the 
eastern half of the solar disk. Hey and Hughes (1955) gave evidence 
that not only the radio but also the Ha radiations were more intense on 
the eastern than on the western half of the disk. This effect is 
referred to as the East- West Asymmetry. Table 4.24 shows the results 
of all the studies we could find in the literature regarding the 
association of Type III bursts with the position of the flares. The 
asymmetry in the distribution of the burst- associated flares may be only 
an effect due to the asymmetry in the distribution of flares. 

In order to analyze the OGO-V data we plotted the frequency 
distribution of burst-associated flares with respect to their central 
meridian distances. We expressed this coordinate in degrees and grouped 
the events in intervals of 10° starting from 0°. Flares for which any 
coordinate was a multiple of 10° were divided into two, each half being 
assigned to adjacent intervals. 
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We adopted three criteria for the selection of the radio events : 

(A) We excluded all events not associated with a flare within 

25 minutes of the radio burst and the five events (Nos. 43, 53, 54, 55 
and 56) suspected to be associated with flares behind the West limb. 

(See Notes to Table A. 2, Appendix 1.) 

(B) The same as (A) except that we did not exclude the behind- 
the-limb events . 

(C) We excluded those cases in Fig. 4.39 for which It - t | > 

1 3-5 Ha 1 

4.5 minutes or that were behind the limb. This rejected not only events 
No. 25, 47 and 63, but also Nos. 43, 53, 5 1 *, 55 and 5 6 (behind the limb) 
and Nos. 4, 15 , 18, 24, 30, 39, 40 , 4l and 52. See Fig. 4.4l. 

The radio data cover the periods l4 March to 25 April 1968 and 
18 June 1968 to 19 February 1970. During this time there were 15 
radio events observed down to or below 0.35 MHz that we did not analyze 
because they were processed after our stop date. During the above 
mentioned periods there occurred approximately 4600 Ha flares. The 
statistical results of the asymmetry found according to the different 
criteria are shown in Table 4.25. 

Table 4.25 

Ha Flares Occurred on the Western Hemisphere of the Sun 

Burst-Associated Ha-Flare* All Ha Flares 

Criterion A Criterion B Criterion C 

68.1* 65.0 % 48.7* 

# 

The association refers to Type III bursts observed by 0G0-Y. 



UNCONFIRMED FLARE 
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The "burst-associated Hot flares exhibit a predominance on the 
western hemisphere even though the total number of Ha flares shows 
a slight deficit in that hemisphere. We observe that the use of 
different criteria of choice practically does not affect the results. 

The average 9^ and standard deviation a of the distribution of 
the burst-associated flares in longitude of Fig. i+ . 1+1 are given in 
Table it. 26. 


Table it. 26 

Distribution in Longitude of Burs t-Associ ate d Ha-Flares 


Criterion 

9° 

_F 

0° 

B 

30 W 

51 

C 

22 W 

43 


These results are the opposite of what has been found from 
ground-based observations, as can be seen in Table lt.2it. To explain 
this difference it could be argued that the 0G0-V results refer only 
to very strong events and therefore are not directly comparable with 
the others, however, the asymmetry found by Swarup et al. (i960) and 
Riihimaa (1963) refer also to major radio events. 

The effect also seems to be independent of the phase of the 
solar cycle. In fact the ground-based observations reported in 
Table 4.24 were obtained during the cycle N°19 that began as a minimum 
on April 1954 and reached the maximum around December 1958. We 
do not have an explanation for this western hemisphere predominance. 

A north-south asymmetry seems to exist also, but it has been 
less studied. Figure it, 42 and Table 4.27 show our results that indicate 
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a high concentration of burst-associated Ha flares on the northern 
hemisphere. During the same period the predominance of all Ha flares 
was also high there. The concentration of the burst-associated Ha flares 
in two bands of latitude as seen in Fig. 1+.1*2 is not unexpected because 
the same phemonenon is characteristic of sun spots. 

The data of Boorman et al . (1961) indicates a northern 
predominance of burst-associated Ha flares but we do not know the 
distribution of all flares during the period of observation. Riihimaa 

(1963) reports north-south asymmetry without mentioning which hemisphere 
was the preferred one. 

This north-south asymmetry was observed also by Harvey and 
Bell (1968), who found that the percentage of microwave bursts associated 
with the northern hemisphere exceeded substantially what was expected 
from the asymmetric distribution of all flares. The data analyzed by 
them covered the period June 1957 to December i960. 

We searched unsuccess fully for a correlation between the 
distribution in latitude of the burst-associated flares and the position 
of the Earth with respect to the solar equator. 

Table k. 2 J 

Ha Flares Occurred on the Northern Hemisphere of the Sun 

Burs t-Associ ate d Ha-Flare* All Ha Flares 

Criterion A Criterion B Criterion ~C ~ 

l6 ' Q% 78 . J% 7 e.J% 63 . 0 ^ 

IT “ 

The association refers to Type III bursts 


observed by OGO-V. 
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4.5-^ X-Ray Flares . In order to study the association of our 
Type III events with x-ray flare events we have used the data on 
Outstanding Events observed from a spacecraft by the University of 
Iowa (Ul) or by the National Research Laboratory (NRL) when the first 
were not available. These data are published in the Solar-Geophysical 
Data, ESSA (now NQAA). The list of outstanding x-ray events is given 
in Table A. 4, Appendix No. 1. 

The time reference for the Type III bursts is again the onset 
of the 3.5 MHz burst. The reference for the x-ray flares is their 
time of start. The distribution of the differences in time of start 
is shown in Fig. 4.1*3) where we have considered only differences with 
absolute value less than 25 minutes. The results indicate that out 
of 64 radio events 40 (62.5%) were associated with x-ray flares. In 
approximately 31 of these 40 cases (T7-5JC)’ the radio bursts followed 
the x-ray, and the most frequent delay was between 1.5 and 2.5 minutes. 
This was close to the results of Teske et al . (1971) who found that on 
the average the soft x-ray rise began 1.4 minutes before the Type III 
bursts in the metric and decametric range. In 25 out of the 40 cases 
( 62 . 5 %) the radio bursts and the x— ray flares were separated by 4.5 
minutes or less . 

We have studied the positional association of x-ray flares with 
our Type III bursts. For this we have had to use the association of 
Hct-flare with the Type III bursts as an intermediate step. 


-2lk- 
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In order to have two degrees of reliability in the temporal 
association we have adopted the criteria (A) and (C) used in the 
analysis of Ha flares. Section U - 5 - 1 • Figures 4.44 and U . 45 show the 
results according to criterion (A) and (C), respectively . We 
observed that some of the Ha flares did not have an associated x-ray 
flare and, conversely, that two x-ray flares are not accompanied 
by Ha flares. In these last two cases the active regions of origin 
were behind the west limb. It is interesting that the radiation in 
the optical range was not visible when the flare was behind the- limb 
but the x-ray radiation was detectable. The explanation is that the 
source of x-ray must be high enough over the photosphere. Teske et_ 
a l • (l97l) have studied the occurrence of Type III bursts accompanied 
by soft x-radiation in the absence of Ha flares. They considered it 
as a phenomenon resulting from a physical process and not from the 
geometry of behind-the-limb events. The real situation may probably 
involve both. 

The presence of an asymmetry in Figs. 4.44 and A. 1+5 is obvious. 
The results of the analysis are shown in Table 4.28 and 4.29. 

Table 4.28 

Type III Bursts Associated with X-Ray Flares 
Occurred on the Western Hemisphere of the Sun 
Criterion A : Criterion C 


70.0 % 


13 . 2 % 



Her (56 CASES, FIG. 4.39) 
X-RAY (40 CASES, FIG. 4.43) 


2.16 



CRITERION A. 


Ha (47 CASES) 
X-RAY (25 CASES) 


•217 
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Table 4.29 

Type III Bursts Associated with X-Ray Flares 
Occurred on the Northern Hemisphere of the Sun 

Criterion A Criterion ■ C 

Bh.0% 85.0% 

Comparing these results with those in Table 4.25 and 4.27 we 
observe that the preference for the wesrt and for the n 0rth seems to be 
more marked for the x-ray flares than for the Ha flares , however we 
cannot malce a definite statement because we do not know the distribution 
in position of all the x-ray flares. As in the case of Ha flares, a 
change in criterion does not appreciably change the statistical results 

for the x-ray flares. The more rigorous conditions tend to accentuate 
the asymmetry slightly. 

— ' 5 ' ^ Particles . These particles are generated 

during or after Ha flares and fall under the general designation of 
solar cosmic rays. The main components are protons and electrons with 
energies of approximately > 1 MeV and > 1*0 KeV, respectively. These 
particles have been classified in the literature according to their 
energy E as cosmic rays or relativistic particles (E > me 2 ) and as 

subcode or nonrelativistic (E < me 2 ). The latter have been In turn 
subdivided into "fast" and "slow." 

The association of solar bursts in the metric range with energetic 
solar particles seems to have been first recognized in connection with 
radio bursts of Types II and IV (Maxwell et al . , 1964). With the advent 
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of artificial satellites in the last decade it has been possible to 
detect these particles directly. This is of special interest to the 
understanding of the mechanism of generation of Type III bursts because^ 
as we have seen in Section U.3.1+, the nature of the exciter particles 
is not definitely known. It is generally accepted however, that they 
are either protons or electrons. The ranges of their energy correspond 
to the ranges of velocity determined for the exciters particles. A 
velocity range between 0.05 and. O .87 c corresponds to protons with 
energies between 1 MeV and 1 GeV, and to electrons between 1 KeV and 
toO KeV. 

The observations of the solar cosmic rays have shown that they 
are guided by the interplanetary magnetic field. It seems logical 
to assume then that the field also guides the Type III exciter 
particles, therefore they would be expected to travel along a spiral 
path. Observations of Type III bursts at very low frequencies (£50 KHz) 
are especially important because the sources may be close to 
1 Astronomical Unit and the times of arrival of the bursts can be 
compared directly with those of the energetic particles detected at 
the Earth. If the exciter particles are identified in a particular 
instance then their velocities would be known. Ideally we could then 
work a fitting process similar to the one described in Section It. 2. 5 to 
determine velocities except that now it would be inverted: instead of 

having the electron density model fixed and the velocity as a variable 
parameter we could have the predictions of times of arrival for fix 
velocities and variable electron density models. In practice we should 
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be warned against delays in time of arrival due to propagation effects 
that can make the identification uncertain, McCracken (1969). 

The fact that Type III bursts can be continuously followed as 
soon as they leave the Sun, makes them valuable to settle questions 
regarding the origin of the two species of particles and the difference 
in the characteristics of their propagation, as discussed by Lin and 
Anderson (19 6j). As a consequence they should be valuable also in the 
study of flare mechanisms or post-flare particle acceleration mechanism. 

In this section we will briefly study the correlation of particles 
events with Type III bursts observed by 0G0-V. 

SolfeC Protons . The phenomenon of interest here is the arrival 
at the Earth of clouds of energetic solar protons. We will refer to 
them as "solar proton enhancements" to differentiate from the so-called 
"Solar Proton Events." This latter name has been given in the 
literature to extremely complex solar events involving the occurrence 
of a wide range of phenomena, solar and terrestrial. During the period 
analyzed here the Solar-Geophysical Data, ESSA, (now NOAA) , gives 
provisional data for four Solar Proton Events: 31 October to 3 November 
1968; 18-20 November 1968; 25 February to 1 March 1969; 30 March to 
2 April 1969. We shall not study these events in detail. 

In the study of solar cosmic ray protons we have used the data 
of the experiment that John Hopkins University and the Goddard Space 
Flight Center had on Explorer 3I* and have on Explorer 35. The 
detector has three windows with energies greater than 60 , 30 and 10 MeV. 
The information is also published in the Solar-Geophysical Data. 
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The data are given in the form of hourly averages and thus are not 
accurate enough to study differences in times of arrival. Our criterion 
for association is necessarily a loose one. We rounded the time of 
start of the radio event to the closest hour, and then selected those 
radio events that preceded the particle event by at most 3 hours, or 
that occur while the particle event was in progress. It takes these 
enhancements hours or days to decay and this makes it more difficult 
to establish an association. The solar proton cosmic ray events 
selected in this way correspond to what the solar particle investi- 
gators call "pronqpt" event, Lin and Anderson (1967). The length of 
the spiral path between' the Sun and the Earth orbit that we have 
adopted (see Section 4.2.2) is approximately 1.3 A.U. This distance 
is traveled in about 22.0 and 12.5 minutes by particles with velocities 
of 0.05c and 0.89c, respectively. During the time between events No. 44 
and 59 » there were no solar proton observations. Out of the 48 radio 
events that occurred during proton-monitoring time 27 (56$) were 
temporally associated with proton enhancements: 13 (27$) preceded them 
and 14 (29$) occurred while the particle event was in progress. The 
positional distribution of the radio events possibly associated with 
proton enhancements is shown in Fig. 4.46. From this figure we deduce 
that 29.5 out of 43 (65.6$) Ha flare- associated radio events occurred 
on the western hemisphere while 29 out of the 27 (74.0$) radio events 
associated with solar protons occurred on that hemisphere. With respect 
to the north-south distribution the situation is as follows: 32 out 

of 43 (74.5$) Ha flare-associated radio events occurred on the North 
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while 18 out of the 27 {&k. 2 %) radio events associated with solar 
proton enhancements occurred on that hemisphere. We do not know 
what was the distribution of all solar proton enhancements. 

Solar Electrons . Anderson and Winckler (1962) detected an 
important x-ray event, and though they did not measure high-energy 
electrons, they implied that electrons associated with the x-ray flare 
were also responsible for the Type III radio bursts. Direct 
measurements of the solar cosmic ray electrons with E > 1*0 KeV 
supported the idea but in connection with solar radio bursts in general 
(van Allen and Krimigis , 1965; Anderson and Lin, 1966). 

The fact that flares producing solar electrons are usually 
associated with Type III bursts was established by Lin and Anderson 
(1967). However the influence of the association of solar electron 
events with Type I bursts, on the distribution on the solar disk, 
asserted by these authors, has been questioned by McCracken and Rao 
(1970). Lin (1968) proposed an identification of solar electrons as 
the exciter particles of Type III bursts. He found that out of 2 k 
electron events with observational radio coverage, 20 (83#) were 
associated with Type III bursts. The radio bursts were ass ume d to be 
associated with the flare if they occurred within -2 minutes of the 
period between the start and maximum of the flare. He noticed that 
the associated Type III bursts tended to be exceptionally intense 
bursts. In a subsequent series of papers Lin found that the emission 
of ~40 KeV electrons is a distinct phenomenon from the emission of 
energetic protons and relativistic electrons. He proposed that low-energy 
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electrons (E < me ) are correlated with Type III bursts while protons 
are correlated with Type II and IV radio bursts (Lin, 1970a,b). 

There are no systematic publications of data on solar cosmic 
ray electrons. The data that follow were obtained from a letter from 
R. P. Lin to F , T. Haddock (March 31, 1971) • Lin examined the 
observations of electron events made from Explorer 35 and IMP-5 during 
the time covered by the 6k radio events analyzed in this work. In 
spite of the fact that he worked with preliminary radio data there are 
very few differences when using the final data. The final results are shown 
in Fig. 1+.1+7. His criteria were as follows: 

(a) A Type III burst was considered associated with an Ha flare 
if the burst occurred within ±3 minutes of the flare start to maximum. 

(b) An electron event was considered to be associated with 

a Type III burst if the electron onset at 1 A.U. was within 1+0 minutes 
of the burst. 

(c) The background radiation was considered too high and the 
comparison invalid if a flux of 10 electrons cm -2 sec" 1 sr" 1 could not 
be observed above the background. 

Out of 55 Ha flare-associated radio events, 2 6 (1+7 .3#) were 
associated with electron events. Out of these 55, 38.5# occurred in 
the western hemisphere; 100# of the ratio bursts associated with 
electron events occurred on that hemisphere. The north-south 
distribution was such that 1+5 out of the 55 (8l.8#) Ha flare-associated 
radio events occurred in the north while 2l+ out of the 2 6 (92.3#) radio 
bursts associated with electron events occurred on that hemisphere. 


WITH ELECTRON EVENT 
WITHOUT ELECTRON EVENT 
HIGH BACKGROUND 



(ELECTRON DATA FROM LIN, 1971 ) 
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It is interesting to observe that all of the electron-associated 
bursts originated on the western hemisphere. This general effect 
could have been expected because of the guidance of the particles by 
the interplanetary magnetic field. In Fig. It. 1*7 we may interpret that 
the electrons produced in flares to the East of approximately 7° West 
of the solar central meridian are guided away from the Earth by the 
spiral field that does not hinder the radio emission. This limiting 
angle depends evidently on the tightness of the spiral and on the 
angular width of the cone of emission of the particles. 

Of special interest is the association of four cases with origin 
located behind the limb. We saw in Section U.5.1 that the two of 
them that were just over the limb (Nos. 1*3 and 53) were associated with 
x-ray flares. The one (No. 1*3) that occurred during time of solar 
proton monitoring was indeed associated with a proton enhancement. 

5 . ^ Ground-Based Radio Observations . In order to find what 
type of radio emission accompanied our Type III bursts , we have used 
the information on spectral radio observations published in the Solar- 
Geophysical Data, ESS A; in the Geophysics and Space Data Bulletin, or 
in data sheets provided by Weissenau Observatory. The burst types as 
defined in the Solar-Geophysical Data are: 

I: Storm bursts. 

II. Slow- drift bursts. 

III. Fast-drift bursts. 

IV. Prolonged continuum. 

V. Brief continuum (normally following Type III bursts). 
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CONT: Continuum in close association with Type III burst 

storms, often with reverse drift bursts and often, 
but not always, associated with noise storms on metric 
wavelength (used by Boulder and Sagamore Hill). 

UNCLF: Unclassified activity. 

The sysbol U stands for a U-shaped burst of Type III. 

The Table A. 5 Appendix 1 shows the individual association of 60 
out of the 61* (93. 8%) bursts observed by OGO-V that occurred during 
ground-based patrol. Out of those 60 events 57 (95-0$) were accompanied 
by activity observed from ground. Table 1+.30 summarizes the results. 


Table it . 30 

Association of 57 Bursts Observed by OGO-V 
with Ground-Based Radio Observations 


Ground-Based Observations 
Type I 
Type II 
Type III 
Type IV 
Type V 
CONTINUUM 
UNCLASSIFIED 


OGO-V Observations (%) 
15.8 
38.6 
98.3 
21*. 6 
36.9 
21.1 
15.8 


The fact that 95$ of the 60 events with ground-based coverage were 
found to have high frequency counterparts is not surprising because 
of the selection criterion. 
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We believe that the events No. 55 and 56 were not accompanied 
by activity observed from ground because they originated at about 120°, 
or further, west of the central meridian. The events No. 1+3, 53, and 5I+ 
are also believed to originate behind the limb but closer to it than 
the events No. 55 and 56. This may explain the fact that they were 
accompanied by high-frequency bursts observed from ground. From 
Table 1+.30 we see that the most ■ frequent type of the associated bursts 
is the type III, as could be expected. 



CHAPTER V. DISCUSSION OF THE RESULTS 

The purpose of this research was to analyze the data on solar 
radio bursts obtained by the OGO-V satellite. Since the wavelengths 
corresponding to our three lowest frequencies of observations exceed 
one kilometer, we have designated the bursts detected in those channels 
as kilometer-waves . 

The data search covers approximately 9200 hours between March 1968 
and February 1970. This period includes the maximum of solar cycle 
No. 20. 

Our data have shown that the Type III bursts reach frequencies as 
low as 0.05 MHz. The high level of radio interference limited the 
number of events observed at that low frequency to one. The detection 
of this event has far reaching implications because it allowed for the 
first time the probing of'the interplanetary space between the Sun and 
the Earth's orbit in an almost continuous way. 

We concentrated on the study of 6k Type III solar radio events 
reaching frequencies equal or lower than 0.35 MHz. This selection 
criterion appears to have led to the choice of the most intense radio 
events. We measured times of start, times of decay and amplitudes of 
these 6k events. We also investigated various solar phenomenon 
associated with them. 

The analysis of the frequency drift-rate led us to a simple 
empirical expression relating frequency drift-rate and frequency. In 
connection with this relationship we introduced the concept of 
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frequency drift-rate index. Combining our results with those of 
other observers we showed that the expression is valid between 550 
and 0.075 MHz. Our formula is a generalization of that proposed 
earlier by Wild. 

We developed a new technique to study the arrival times of 
the bursts . This consists in the prediction of arrival times at the 
Earth and the comparison with the observations. The predictions are 
based on several fundamental assumptions of which we will mention 
the most important. First, we have assumed that the bursts are 
emitted at the fundamental and/or second harmonic of the local plasma 
frequency. We considered that at present this is the best assumption. 
Second, we assumed that the exciter particles travel with constant 
velocity along an Archimedes spiral situated on the plane of the 
ecliptic. The constancy of the velocity seems to us the most tenable 
hypothesis at the present time; we have not found evidence strong 
enough to reject it. The adopted shape of the trajectory seems 
reasonable in view of the results of studies of solar high-energy 
particles and the solar wind. For simplicity we assumed that the 
trajectory lies m the ecliptic. We have found no evidence in favor 
or against this. The lack of good time resolution in the lowest frequency 
channels prevented us from detecting out-of-plane effects. Third, 
to have a self-consistent model we adopted an electron density distri- 
bution derived from the frequency drift-rate expression mentioned above 
and from some of the assumptions now under discussion. Also, the 
electron density distribution was chosen to fit the observations near 
the Sun and near the Earth. In this assumption we depart from the 
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current belief that the Type III bursts occur in dense solar streamers 
and we assume that they occur in the solar wind itself, at least for 
these low-frequency bursts. These two propositions may not be 
incompatible because our results refer to distances from the Sun 
larger than about 20 solar radii and the others refer to distances 
in general less than that. Fourth, we disregarded the effects of 
the coronal plasma on the propagation of the radio waves. This may 
be the most debatable of our assumptions. We adopted it because of its 
simplicity and because the changes in arrival times stemming from it 
would be undetectable in our data due to the lack of good time resolution. 
Finally, we did not consider the effects of magnetic fields except 
for the guiding of the exciter particles along a spiral trajectory but 
not spiraling about it. Again, this assumption was made to simplify 
the analysis . 

The use of the technique for determining the exciter particles 
velocity revealed an unexpected preponderance of the second harmonic 
emission below about 1 MHz. We found that as the frequency of 
observation decreases the times of arrival of the bursts fits better 
the arrival time of the second harmonic than that of the fundamental 
emission. At 0.05 MHz about 100% of the radiation seems to be emitted 
as second harmonic. This phenomenon practically doubles the distances 
from the Sun at which the observed radiation was supposed to originate 
and it places the source of the 0.05 MHz burst near. the Earth's orbit. 

The recognition of this phenomenon led us to derive velocities of the 
same order as those measured from ground-based observations by direct 
methods , without requiring models several times denser than the solar 
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wind. The distribution of the velocities has an average of 0.32c and 
a standard deviation of 0.09c. There is some suggestion in the 
distribution for two preferred ranges of velocities centered approxi- 
mately at 0.l8c and 0.36c. 

We found a simple three-parameter formula to express the 
electron density of the solar corona as a function of distance from 
the Sun between approximately a few solar radii and the Earth's orbit. 

The parameters can be determined from the coronal electron densities 
near the Sun and the Earth, and from the measured frequency drift-rate 
index of Type III bursts. Incidentally, the formula can be used also 
with good approximation to fit optical determinations of electron 
densities in individual streamers. 

The OGO-V Radio Astronomy experiment was not designed to measure 
absolute flux density values. However, by making certain approximations 
we attempted a study of the burst spectra. The basic assumptions were 
that the antenna-spacecraft structure behaves like a short center-fed 
dipole and that the bursts' radiation is unpolarized. 

We found that, in general, the derived flux densities increase 
with decreasing frequency. The peak-power ' spectra usually exhibit one or 
two peaks. We suspect that the double peaks may be associated with 
the presence of second harmonic emission at lower frequencies and with 
directional effects at the emission or during propagation. We computed 
energy densities and found a function that fits approximately the 
average energy density spectra observed. This function is a 
generalization of that proposed previously by Wild. The measured 
flux densities and energies exceed by several orders of magnitude the 
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values obtained from ground-based instruments at higher frequencies. 

The main peak of the few events studied appears to occur on the average 
at about 0.3 MHz. 

We measured the times of decay of selected bursts. The 
measurements were made in idealized curves drawn through their complex 
time profiles , without an attempt to separate fundamental from second 
harmonic emission. We defined the time of decay at a certain frequency 
as the time between the peak of the burst and one e-th of it, 
regardless of whether the decay was exponential or not. Combining our 
observations with others obtained at frequencies up to about 200 MHz 
we established an empirical relationship between decay time and 
frequency. This is practically the same that Wild found earlier in 
the limited range between 80 and 120 MHz. At 100 MHz it takes a burst 
about one second to decay while at 0.05 MHz it takes it about 1000 
seconds. Our results confirmed the inadequacy of the proton-electron 
collisional damping mechanism to explain the decay time of Type III 
bursts at these low frequencies. Using radio observations in the 
full frequency range available and our electron density models we 
established an empirical expression according to which an effective 
collision frequency must increase with distance to produce the observed 
burst durations. This is approximately the same relationship that 
theoreticians have required to make the theory fit the temperature 
observations of solar wind protons and electrons near the Earth. 

We studied the association of the kilometric-wave ' bursts 
with Ha flares. X-ray flares, solar proton enhancements, solar electron 
events, and ground-based observations. The two most interesting 
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results are . first , all of the hursts associated with solar cosmic 
ray electron events originated on the western solar hemisphere and 
second, the Type III solar bursts observed from OGO-V were associated 
with the flares occurring predominantly in the western solar hemisphere 
This last result is opposite to the asymmetry found from ground-based 
observations at higher frequencies. 

Our data indicate that the lowest frequency reached by a burst 
depends on the position of the exciter particles origin with respect 
to the Sun's central meridian. The further westward from the east 
limb that a flare occurs the lower is the lowest frequency observed. 

We have interpreted this either as a directional effect at the emission 
or as refractive effect during the propagation of the radio waves. 

Our results on the prevalence of second— harmonic emission, on 
the times of decay, and on the flux densities may have important 
implications for theories of the emission mechanism of Type III bursts 
and of the propagation of the radio waves through the corona. 

The consistency of our results , within the accuracy of the 
measurements , lends support to some of the assumptions made for the 
analysis, notably, the validity of the local plasma hypothesis, the 
contancy of the exciter particles velocity, and the spiral shape of 
their trajectory. 



CHAPTER VI. RECOMMENDATION FOR FUTURE ANALYSIS 

The analysis of the kilometric-wave bursts observed by 
OGO-V suggests several topics for future investigation. Some are 
refinements of the research presented in this work, others are new 
areas of study. 

Except for problems of theoretical character most of the 
rest would require new instrumentation. We estimate that a stepping 
radiometer like the one used in this research has much to offer yet, 
therefore we will not discuss other types. We will assume that the 
instrument would be aboard a spacecraft stabilized in the same way 
as OGO-V. 

The emphasis in a new instrument should be laid upon the 
following topics that will be discussed subsequently: 

(a) Higher time resolution, 

(b) A less noisy spacecraft, 

(c) Higher system sensitivity, 

(d) Selection of observing frequencies, 

(e) Bandwidth requirements, 

(f) Internal functions. 

Conditions (a) and (b ) would allow a better identification of 
the components of an event at the different frequencies especially 
with regard to the nature of the emission: fundamental at a certain 
frequency or second harmonic at that frequency. A higher time 
resolution is needed mainly in the high frequency end, near 3.5 MHz, 
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where the time scale of the burst is still short. The accuracy in 
the measurements of the times of arrival could help in establishing 
whether the exciter particle trajectories are radial or spiral. 

Conditions (b ) and (c) would permit the detection of a large 
number of weak events down to the lowest frequencies. As we have 
shown, the spectrum drops rapidly below 0.1 MHz. It would be of 
interest to know to what extent, if any, the selection criterion used 
for OGO-V data affected our conclusions. Of special importance would 
be the possibility of working with less complex burst events, if the 
weaker bursts are indeed simpler. In our analysis complexity was a 
major source of confusion in the identification of components and modes 
of emission. Of particular interest would be the measurement of decay 
time of the radiation in the fundament al mode at a certain frequency and 
in the second harmonic mode at the same frequency. 

Conditions (a) through (d) would allow a more accurate 
determination of the range of frequencies in which the predominance 
of the second harmonic sets in. It would be of interest to investigate 
if there is a correlation between these frequencies and the position 
of the associated flares. This might result in information regarding 
the directivity pattern at the emission or regarding refractive effects 
on the propagation of the radio waves. 

The selection of frequencies, (d), may depend on the main objective 
of the experiment. For example, if we want to study the second harmonic 
prevalence we may want to have several channels more or less closely 
spaced m the upper half of the range, from 3-5 MHz to 0.35 MHz. 

Otherwise we may want to have them concentrated in other ranges. 
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Condition (e) refers to the possible need of a narrower 
receiver bandwidth at the lower frequencies, below about 0.1 MHz, if 
the; time scale of the events at those frequencies is less than the 
time it takes the exciters particles to traverse the range of plasma 
frequencies determined by the extremes of the receiver bandpass. 
Measurement of the rise time of the bursts may give information about 
the ambient solar plasma and/or the spatial distribution of the exciter 
packet. 

We could search for the possibility that the frequency drift 
rate depends on the intensity of the burst, as suggested by Rabben 
(1960b). We could test whether the lowest frequency observed in a 
burst is a function of the position of the associated Ha flare, as 
suggested by our results , and whether the suggestion of bimodality 
in the distribution of exciter particles velocities is real or not. 

The most important function included in (f) is the inflight 
calibration. The OGO-V radiometer stability was better than expected, 
therefore a similar instrument would not need such frequent calibration 
which would reduce the probability of coincidence with a burst. 

An improved radiometer would probably make it necessary to 
incorporate in the analysis refinements related to the effect of the 
solar plasma on the propagation of the radio burst radiation. This 
would improve upon the accuracy of the deduced velocities and would 
test the validity of the ideas proposed by other workers, for example, 
the "echo hypothesis." 
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So far we have assumed that the lowest frequency of observation 
is 0.05 MHz. A lower frequency would be better, especially if 
conditions (b) and (c) are fulfilled. In this range the bandwidth 
could become a problem because it should be narrow enough not to smear 
the time profiles . 

Of special importance would be the study of those lype III 
events that reach 0.05 MHz or lower, and are associated with energetic 
solar particle events observed near the Earth by space instruments. 

A precise knowledge of the mode of emission, of the particle velocity 
and trajectory would allow the calculation of times of arrival of 
the particles at the Earth and the comparison with the direct detection 
of solar energetic particles. This would answer the question about 
the nature of the particles. 

In dealing with individual events it would also be of interest 
to undertake a joint study of selected iype III bursts with a ground- 
based radio observatory that has the capability of measuring directly 
the velocities. This would give information about the coronal electron 
density distribution. For similar reasons it would be of interest to 
study individual events observed by 0G0-V and other spacecraft that 
may provide supplementary information, for example, IMP-6. 

The analysis in this study points out the need of theoretical 
work, especially in two areas. The first is an extension of the theory 
of lype III mechanisms down to very low frequencies. The theories 
should account for the second harmonic phenomenon and for the high 
intensities observed at very low frequencies. Of special importance 
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for the theories would be experimental measurement of the size of 
the Type III sources, the scattering of the radiation, the absolute 
or relative intensities of the radiation at the different frequencies 
and modes, the energy distribution of the exciter particles, etc. 

The second area is the study of the mechanisms that shorten the decay 
time of the bursts. We have given quantitative information on an 
effective collision frequency. This should be combined with other 
observations that give information concerning the size of the electron 
density inhomogeneities or the state of turbulence of the solar wind. 
For example, studies based on the occultation of radio sources by the 
solar oorona, on the observation of pulsars, and on radar observations 
As far as we know radio observations and deep space probes provide the 
only means available today to study the corona between the Earth and 
a few solar radii from the Sun. 



APPENDIX NO. 1 


Table A.l 


List of Events Observed by OGO-V Below 0.6 MHz 


Event 

No. 

Date 


Time ^ ^ 
(U.T • ) 

Lowest 

Frequency 

(MHz) 

(2 ) 

Spacecraft v ’ 
Distance 
(10 3 km) 

r 

9 ° 

Note 

1 

28 Mar. 

1968 

0321.0 

0.1 

54 

LT\ 

OJ 

54 


2 

l4 Aug. 

1968 

1323.0 

0-35 

l4o 

173 

31 


3 

6 Sept. 

1968 

0916.6 

0.35 

149 

143 

39 


4 

13 Sept. 

1968 

1302.2 

0.35 

109 

120 

50 

(3) 

5 

26 Sept. 

1968 

1925.1 

0.2 

131 

115 

45 


6 

27 Sept. 

1968 

1033.8 

0.35 

149 

126 

35 

(3) 

7 

29 Sept. 

1968 

0929.6 

0.1 

129 

112 

46 


8 

29 Sept. 

1968 

1621.4 

0.35 

145 

118 

4i 


9 

4 Oct. 

1968 

0000.7 

0.1 

47 

52 

60 


10 

4 Oct. 

1968 

2122.0 

0.35 

145 

113 

4i 


11 

6 Oct. 

11968 

1723.8 

0.2 

72 

77 

57 


12 

11 Oct- 

1968 

0306.3 

0.1 

113 

123 

22 


13 

17 Oct. 

1968 

08l6.4 

0.2 

106 

. 85 

51 

(3) 

l4 

23 Oct. 

1968 

2001.4 

0.1 

136 

105 

28 


15 

27 Oct. 

1968 >1312.8 
<1313.6 

0-35 

77 

58 

56 

(M 

16 

29 Oct. 

1968 

0953.9 

0.2 

103 

107 

19 


17 

1 Nov. 

1968 

1659-7 

0.35 

69 

46 

58 


18 

2 Nov. 

1968 

0949.6 

0.35 

l4i 

82 

42 


19 

4 Nov. 

1968 

0518.0 

0.35 

51 

23 

60 


20 

13 Nov. 

1968 

1246.4 

0.1 

142 

83 

30 


21 

18 Nov. 

1968 

1027.7 

0.03 

148 

75 

55 
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Table A . 1 (cont . ) 


( 2 ) 

/, \ Lowest Spacecraft ' 
Event Time' ' Frequency Distance 


No. 

Date 

(U.T.) 

(MHz) 

(10 3 km) 

r 

9° 

Note 

22 

23 Nov. 1968 

>0531.5 

<0540.5 

0.1 

142 

62 

42 

(4) 

23 

25 Nov. 1968 

0354.6 

0.35 

76 

28 

57 

(3) 

24 

9 Dec. 1968 

0622.7 

0.2 

148 

52 

35 


25 

20 Dec. 1968 

1546.3 

0-35 

~113 

60 

10 


26 

24 Dec. 1968 

2227.7 

0.1 

147 

37 

33 


27 

3 Jan. 1969 

0316.4 

0.35 

68 

342 

60 


28 

6 Jan. 1969 

0907.3 

0.2 

137 

15 

43 


29 

17 Jan. 1969 

1704.7 

0.35 

133 

18 

26 


30 

24 Jan. 1969 

0717.0 

0.1 

116 

348 

49 


31 

8 Feb. 1969 

04i5.3 

0.35 

50 

15 

-5 


32 

9 Feb. 1969 

1431.9 

0.35 

147 

347 

56 


33 

9 Feb. 1969 

1724.7 

0.2 

l46 

349 

33 


34 

10 Feb. 1969 

0343.8 

0.35 

130 

355 

24 


35 

23 Feb. 1969 

o44o.9 

0.1 

128 

342 

23 


36 

24 Feb. 1969 

1427.2 

0.2 

123 

319 

48 


37 

24 Feb. 1969 

2308.1 

0.2 

142 

327 

4o 

(3) 

38 

25 Feb. 1969 

0909.7 

0.1 

145 

334 

32 


39 

27 Feb. 1969 

1403.8 

0.35 

143 

325 

4o 


4o 

27 Feb. 1969 

l4i2.0 

0.35 

i43 

325 

4o 


4i 

12 Mar. 1969 

1752.1 

0.35 

145 

314 

37 


42 

13 Mar. 1969 

2117.7 

0.35 

72 

335 

1 


43 

30 Mar. 1969 

0249.6 

0.1 

88 

267 

57 

(cont. ) 
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Table A.l (cont. ) 


Event 

No. 

Date 

Time 
(U.T. ) 

Lowest 

Frequency 

(MHz) 

(2) 

Spacecraft v ' 
Distance 
( 10 3 km) 

r 

9 ° 

Note 

44 

3 May 1969 

1942 . 5 

0.2 

145 

264 

36 


45 

4 May 1969 

0719.8 

0.1 

133 

271 

25 


46 

6 May 1969 

0639.5 

0.1 

143 

250 

39 


47 

11 May 1969 

1820.3 

0.2 

l 44 

258 

33 


48 

25 May 1969 

o 64 i .5 

0.35 

122 

253 

20 


49 

28 May 1969 

2310.2 

0.1 

91 

209 

58 


50 

29 May 1969 

0020.8 

0.2 

99 

214 

55 


51 

29 May 1969 

1451.5 

0.2 

l 4 i 

235 

4o 


52 

29 May 1969 

1938.4 

0.2 

144 

238 

36 


53 

31 May I969 

2101.5 

0.2 

124 

224 

48 


54 

1 June 1969 

1700.1 

0.1 

l 4 i 

24 0 

30 

( 3 ) 

(M 

55 

2 June 1969 

0854.2 

0.35 

91 

251 

8 


56 

2 June 1969 

0914.5 

0.1 

89 

252 

7 


57 

11 June 1969 

1619.8 

0.1 

142 

223 

39 


58 

l8 June 1969 

2109.5 

0.1 

106 

198 

54 


59 

l 8 Ju&e 1969 

2338.5 

0.1 

116 

203 

51 


60 

5 July 1969 

1225.0 

0.2 

l4o 

207 

29 


61 

25 Oct. 1969 

1534.9 

0.1 

126 

102 

20 


62 

7 Nov. 1969 

0346.0 

0.1 

i4o 

82 

33 

( 3 ) 

(M 

63 

25 Nov. 1969 

1303.4 

0.2 

137 

66 

30 

64 

19 Feb. 1970 

0801.2 

0.2 

137 

338 

29 

(cont . ) 
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Table A . 1 (coni.) 


Notes : 

(1) Time of start of the event at 3-5 MHz. 

(2) Approximate values taken from the 0R/3A. The distance is geocentric. 
0 is the angle from Earth-Sun line to the projection of earth- 
spacecraft line onto ecliptic plane. 9 is the angle between Earth- 
spacecraft line and ecliptic plane. In Fig. 2 . 9 , 0 = $ SEU, 9 = ^ 
TEU. 

(3) Burst probably observed also in next channel below. 

(4) Event No. 15 - Calibration at the beginning of the burst. 

Event No. 22 - Data gap at the beginning of the burst. 

Events No. 5^ and 62 - Time corresponds to the start of the 

event at 0.9 MHz. The event was not observed at the two higher 
frequencies . 



APPENDIX NO. 1 


Table A. 2 

Association of OGO-V Events with Hen Flares, Heliographic Latitude B 

o 


Event 

No. 

of the Solar Disk 

Position of the 
Associated Flare 

SID 

Center 4 ** 

McMath 

Plage 

No. 

at the Time 

Flare Start 
Time 

(U.T.H+) 

of the Events 

Range of 
Flare 
Importance 

B (°) 

0 v ' Note 

1 

S 13 

W 47 

9273 

0321 

IN, IB 

-6.72 

2 

N 12 

W 79 

9567 

1322 

IB 

+6.62 

3 

N 12 

w 51 

9630 

0917 

-N, 2B 

+7-25 

4 

S 16 

E 90 

9668 

1309 

-F 

+ 7.22 ( 1 ), ( 2 ) 

5 

S 21 

E 55 

9692 

1925 

-F, IN 

+6.91 

6 

N 10 

E 12 

9687 

1033 

-B 

+6.87 

7 

N 13 

w 13 

9687 

0930 

IN, 2B 

+6.78 

8 

N 17 

w 52 

9678 

1619 

IN, 2 B 

+6.78 

9 

S 16 

W 38 

9692 

2357 

2N, 3B 

+ 6.54 (3) 

10 

N 20 

E 4 o 

9705 

2120 

2B 

+6.54 

ll 

N l4 

w 19 

9700 

1722 

IN, IB 

+6.43 

12 

N 18 

w 46 

9705 

0304 

IN 

+6.12 

13 

N 21 

w 01 

9726 

08l6 

-N, -B 

+ 5*68 

14 

N 20 

W 90 

9726 

2001 

-N, -B 

+ 5-18 

15 

S 17 

E 17 

9740 

1300 

2N 

+4.81 

16 

N 14 

w 88 

9733 

0953 

-N, -B 

+4.62 (l) 

17 

S 12 

1 

w 61 

9740 

1656 

-N 

+4.32 (2) 

18 

s l 4 

w 66 

9740 

0940 

2B 

+4.22 

10 

s 15 

w 00 

074o 

031 5 

IB, 2B 

+ 4.01 

20 

N 18 

W 22 

0760 

12') 6 

-N, IB 

+ '). 02 


(cont. ) 
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Table A. 2 (cont.) 


Event 
No . 

Position of the 
Associated Flare 

5(°) -\(°) 

McMath 

Plage 

No. 

Flare Start 
Time 

. (U.T.)(t) 

Range of 
Flare 
Importance 

B (°) 
0 v ' 

Note 

21 

N 21 

w 90 

9760 

1026 

IB 

+2.43 


22 

S 12 

E 27 

9780 

0538 

-N 

+1.82 

(2) 

23 

S 11 

E 01 

9780 

0335 

-N 

+1.57 


24 

N 11 

W 46 

9798 

0613 

3F, IN 

-0.21 


25 

s 19 

N 22 

W 12 

w 100 

9829 

9821 

-- 

-- 

-1.6l 

(M 

26 

N 19 

E 28 

9842 

2228 

IB, 2B 

-2.11 


27 

N 24 

E 68 

9855 

0320 

-N 

-3.30 


28 

N 20 

E 28 

9855 

0907 

IB 

-3.64 


29 

N 17 

E 46 

9873 

1704 

-B, IB 

-4.79 


50 

N 20 

w 08 

9879 

0706 

3N, 2B 

-5.42 

(5) 

31 

N 12 

E 06 

9911 

0415 

-N 

-6.50 


32 

N 18 

E 02 

9918 

1431 

-N, -B 

-6.55 


33 

N 18 

E 00 

9918 

1723 

-B 

-6.55 


3 i+ 

N 18 

w 06 

9918 

0343 

-N, -B 

-6.6l 


35 

N 12 

w 09 

9946 

o44i 

-N, IB 

-7.11 


36 

N 11 

w 26 

9946 

1432 

-N, IB 

-7.14 


37 

N 12 

W 31 

9946 

2310 

2B 

-7.14 


38 

N 14. 

W 38 

9946 

0909 

2N, 2B 

-7.16 


39 

N 13 

w 64 

9946 

1355 

IN, 2B 

-7.19 

(6) 

4o 

N 13 

w 64 

9946 

1355 . 

IN, 2B 

-7.19 

(6) 

4l 

N 12 

w 80 

9966 

1739 

2B, 3B 

-7.21 


42 

S 14 

E 45 

9988 

2117 

-N, -B 

-7.20 


43 

N 09 
N 22 

E 50 

w no 

10114 

9994 

-- 

-- 

-6.63 

(M 


(cont. ) 


-2i 16- 


Table A- 2 (cont. ) 


Event 

No. 

Position of the 
Associated Flare 

sn x(°) 

McMath 

Plage 

No. 

Flare Start 
Time 

(u.T.Ht) 

Range of 
Flare 

Importance 

B (°) 
0 ' 

Note 

44 

N 06 

w 59 

10057 

1945 

IN 

-3.9^ 


45 

N 10 

w 61 

10057 

0720 

-F, IN 

-3.83 


46 

N 10 

w 90 

10057 

0637 

--F 

-3.63 

(2), 

47 

-- 

-- 

— 

-- 

-- 

-3.08 

(M 

48 

N 11 

E 21 

10109 

o64o 

-F, IF 

-1.47 


49 

N 11 

w 63 

10109 

2310 

IN 

-1.11 


50 

N 11 

w 64 

10109 

0020 

IB 

-0.99 


51 

N 11 

W 73 

10109 

1452 

IB 

-0.99 


52 

N 10 

w 76 

10109 

1938 

IB 

-0.99 


53 

N 12 

w 110 

10109 

-- 

-- 

-0.75 

(M 

5^ 

N 12 

W 120 

10109 

-- 

-- 

-0.63 

(M 

55 

N 12 

w 130 

10109 

-- 

-- 

-0.51 

(*0 

56 

N 12 

w 130 

10109 

-- 

-- 

-0.51 

(M 

57 

N 10 

W 20 

10134 

1620 

2N, 3B 

+0.57 


58 

N 07 

E 09 

10148 

2109 

-N, -B 

+l.4i 


59 

N 08 

E 05 

10148 

2337 

-N, -B 

+1.4i 


60 

s 13 

E 16 

10181 

1226 

IN, 2B 

+3.34 


61 

N 10 

w 18 

10381 

1535 

-N 

+5.02 


62 

N 13 

E 11 

io4o6 

0324 

IF, 2N 

+3.72 


63 

N 11 

w 47 

10452 

1225 

IN 

+1.60 

(M 

64 

s 13 

w 31 

10579 

0802 

£5 

1 

-6.99 
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Table A .2 (cent.) 


From Solar-Geophysical Data, Environmental Science Services Administration 
(ESSA), now National Oceanic and Atmospheric Administration (NOAA). 


** 

From the American Ephemeris and Nautical Almanac, 
(t) 

Time closest to the start of the 3*5 MHz burst. 
Notes : 


(1) Unconfirmed. 

(2) Only one observation reported. 

(3) Time corresponds to previous day. 

(4) Because of the selection criterion all the radio events studied were 
extremely strong and except in very few cases they were always related 
to intense and complex optical flares. Usually there were several 
flares from the same McMath plage, starting within a few minutes 

of our radio events. The times given here indicate how close the 
group of flares was to the radio event. Unless we had no choice 
we did not consider flares that were already in progress for times 
longer than about 10 minutes since our radio events started. The 
Importance is given as a range when there were discrepancies in its 
estimation. The minus (-) sign indicates a subflare. The position 
corresponds to the average of the group. 

Event No. 25 -- Unconfirmed flare -F at S 19 W 12. At 1501 it was 
in progress. The McMa'th plage region No. 9821 was very active from 
December 16 until December 19 when it disappeared behind the limb. 

The longitude 100°W is extrapolated. 

Event No. 43 -- Unconfirmed flare -B N 09 E 50 starting before 02 49 . 
Ground base radio observations give evidence that the flare occurred 
at the McMath plage No. 9994 behind the limb (V. L. Badillo, J. E. 
Salcedo, 1969 )* The longitude 110°W is extrapolated. 

Event No. 47 -- No flares between 1711 and 2056 hrs. No flare patrol 
between 1900 and 1905- All confirmed flares on the day of this event 
and on the previous day came from the eastern hemisphere the Sun, 
almost all of them from McMath. No. 10088. The estimated average 
position is E 50. 

Events No. 53, 54, 55, 56 -- No flare occurred near the time of these 
events except 3 unconfirmed weak subflares on June 2, starting between 
0845 and 0850 and ending before 0915, at N 12 E 79* Evidence given 


(cont . ) 
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Table A. 2 (cont.) 

in the text suggests that these four events originated at the 
McMath plage No. 10109 when it was behind the limb and that they 
belong a sequence started with event No. 49. This has been 
supported later on by solar electrons observations by R. P. Lin 
(1971). He found correlations with these four events. Longitudes 
are extrapolated. 

Event No. 63 -- This confirmed flare continued after 1246. 

(5) Earliest time of observation. At this time the flare was in progress. 

( 6 ) Same flare. 



APPKNIVI X NO. I. 


Table A.') 

Velocities Derived from the Model of Fig. 4.16 


Event 

No. 

Unambiguous 


Ambiguous 


P 



V P L 

1 

0.36 




2 


0.17 

0.33 

1.94 

3 

0.43 




]+(!) 


-- 

-- ' 


5 


0.20 

0.40 

2.00 

6 


o.l4 

Q. 35 

2.50 

7 

0.36 




8 

0.18 




9 


0.19 

0.36 

1.90 

10 


0.20 

0.42 

2.10 

11 

0.19 

•• 



12 


0.21 

0 . 4 i 

1-95 

13 

0.39 




14 ( 7 ) 

0.46 




15 


0.20 

0.34 

1.70 

16 

0.32 




17 


0.13 

0.28 

2.15 

18 


0.17 

0.30 

I.76 

19 ( 2 ) 


0.29 

0 . 4 o 

1.38 

20 


0.17 

O.29 

1.71 

21 ( 3 ) 


0.18 

0.36 

2.00 

22 


0.18 

0.34 

1.89 

23 

0.19 




24 


0.24 

0.36 

1.50 

25 ( 2 , 4 ) 


0.17 

0.30 

1.76 

26 


0.16 

0.30 

1.88 

27 


0.18 

0.36 

2.00 





(cont. ) 
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Table A. 3 (cont. ) 


Event 

No. 

Unambiguous 


Ambiguous 


P 

P L 

9. 

6 /e 
p h /f l 

28 

0.20 




29 


0.20 

0.42 

2.10 

30 


0.15 

0.24 

1.60 

31 


0.20 

0.32 

1.60 

32 

0.49 




33 

0.43 




34 


0.25 

o.4i 

1.64 

35 

0.32 




36 

0.18 




37 

0.30 




38 

0.26 




39 

0.36 




4o 

0.46 




4l 

0.4i 




42 


0.10 

0.26 

1.63 

43(4) 

0.32 




44 

0.26 




45 

0.36 




46(2) 


0.18 

0.36 

2.00 

^7(5) 


0.13 

0.22 

1.69 

48(2) 

0.33 




49 

0.29 




50 


0.16 

0.32 

2.00 

51 


0.18 

0.25 

1.39 

52(2) 

0.37 




53(6) 


0.18 

0.37 

2.05 

54(6) 

0.34 




55(6) 


0.15 

0.30 

2.00 

56(6) 

0.4o 




57 


0.18 

0.30 

1.67 

38 

0.37 




50 

0.32 





(eont. ) 
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Table A -3 (cont. ) 



Unambiguous 


Ambiguous 


Event 

No. 

6 

P L 


6 /p 
M H /P L 

60 

0.28 




6l 


0.18 

0.34 . 

1.89 

62 

0.36 




63 


0.16 

0.32 

2.00 

64 

" 0.19 




Notes : 






(l) Useless to determine velocity. 


(2) Somewhat large spread in the velocities of the components. 

(3) November l8, 1968 event. 

(4) Two possible associated flares. For event No. 25 we accepted 9 = 100° 

for event No. 43 we accepted 9^ = 110°. See Table A. 2. 

(5) The velocity was determined assuming an average 9 = -50. 

See Table A. 2. 

(6) The associated flares were assumed to be behind the west limb. 

See Table A. 2. 

(7) Almost at the end of this work we found that due to mistake in the 
date this . event was misidentified. The position of the wrongly 
associated Ha flare was S 12 E 60 while the correct one was N 20 W 90. 
This error leads to a velocity of 0.22c. The wrong information was 
used in Figs. 4.23, 4.24 and A. 2 from which the computational results 
contained in Tables 4.15., A. 6 and A. 7 were obtained. We estimated that 
the use of the correct values would not change the result appreciably 
therefore we did not modify them. 



appendix no. 1 


Table A. 4 


Outstanding 

X-Ray Events ^ Associated 

with the Radio 

Events Observed 


by 

OGO-V 



Radio 
Event No. 


X-Ray 

Event 



Start 

Max. 

End 

Observer^ 


U.T. 

U.T. 

U.T. 

1 

03^9 


0522 

N.R.L. 

2 

1324 

1352 


U.I. 

3 

0917 

0921 


IT 

6 

1033 

1039 


tt 

7 

0925 

0942 


It 

8 

1616 

1624 


t! 

9 

0010 


0011 

N.R.L. 


2324 ( 3 ) 

0011 


U.I. 

10 

2132 


2133 

N.R.L. 


-- 

2155 


U.I. 

n 

1721 

1736 


It 

15 

1307 

1341 


ft 

l 6 

0944 

0959 


II 

18 

0947 

1005 


It 

19 

0520 


0540 

N.R.L. 


0518 ( 4 ) 

0529 


U.I. 

21 

1026 

1057 


ft 

22 

0538 

0543 


ft 

24 

0613 


0711 

N.R.L. 


0610 

0620 


U.I. 

26 

2227 


2240 

N.R.L. 

28 

0907 


0918 

N.R.L. 


0906 

0913 


U.I. 

29 

1704 


i7.ll 

N.R.L. 


170 

IfOft 


ll.l . 





(cont . ) 
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Table A .4 (cont . ) 


Radio 



X-Ray Event 


Event 

Start 

Max 

End 

(?) 

Observer 

No. 

U.T. 

U.T. 

U.T. 

30 

0709 


1031 

N.R.L. 


0658 

0740 


U.I. 

31 

04i3 

o4i8 


It 

32 

1432 


1438(6) 

N.R.L. 


1430 

1^35 


U.I. 

33 

1724 


1736 

N.R.L. 


1723 

1726 


U.I. 

3k 

0343 

0346 


It 

35 

o44o(5) 


0508 

N.R.L. 


o442 

0451 


U.I. 

36 

1^33 


1445 

N.R.L. 

37 

2330 


0112 

N.R.L. 


2306 

2322 


U.I. 

38 

0856 

0918 


tt 

39 

l4oo(5) 


1735 

N.R.L. 


1356 

1412 


U.I. 

4o 

i4oo(5) 


1735(6) 

N.R.L. 


1356 

1412 


U.I. 

4i 

1739(5) 


1954 

N.R.L. 


1737 

17^5 


U.I. 

^3 

0257(5) 


0904 

N.R.L. 


0246 

0302 


U.I. 

44 

1959(5) 


2016 

N.R.L. 


1947 

2006 


U.I. 

45 

0730 


0744 

N.R.L. 

46 

. 0638 


0716 

II 


0636 

0648 


U.I. 

50 

0019 


0031(6) 

N.R.L. 


0019 

0027 


U.I. 

51 

1449 


1518(6) 

N.R.L. 


1449 

1500 


U.I. 


(cont. ) 


-254- 


Table A. 4 (cont. ) 


Radio 


X-Ray 

Event 


Event 

Start 

Max. 

End 


No. 

U.T. 

U.T. 

U.T. 

(2) 

Observer v ' 

52 

1939 


2013 

N.R.L. 


1938 

1946 


U.I. 

53 

2057 


2300 ( 6 ) 

N.R.L. 


2103 

2126 


U.I. 

57 

1617 


1834 

N.R.L. 


l 6 l 6 

1630 


U.I. 

60 

1229 


1244(6) 

N.R.L. 


1225 

1238 


U.I. 


Notes : 


(1) Published in Solar-Geophysical Data, E.S.S.A. 

(2) N.R.L. stands for Naval Research Laboratories and U.I. stands 
for University of Iowa. 

(3) Flare started in the day previous to the day of the maximum. 

(4) Start time given corresponds to rapid rise. Onset occurred at 

0458 . 

(5) The flare was in progress when observation started. 

( 6 ) The flare was in progress when observation ended. 
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Table A . 5 


Association of Ground-Based Radio Observations with the Radio Events 

i 




Observed by 

OGO-V 




Event 








No. 

I 

II 

III 

IV 

V 

CONT. ' 

UNCLF 

1 


sj 

sj 


sj 


sj 

2 



sj 


sj 

sj 


5 



sl 





4 

No activity 

recorded. 






5 



sl 





6 



sj 


sl 



7 



sj 

sl 

sj 



8 


sj 

sl 



sj 


9 


sj 

sj 

I 


sj 


10 



sl 



sl 

sj 

11 


sj 

sj 





12 



sj 





13 





sj 



l 4 

sl 







15 



sj 

sj 




16 



sj 


sj 



17 


sj 

sl 

sl 



sj 

l8 


si 

sj 

sj 



sj 

19 


sl 

sj- U 


sj 

sl 


20 



sl 





21 


sj 

V 

sl 




22 



sj- U 


sj 



23 

si ' 


sj 


sj 


sj 

24 



sl -u 


sj 



25 



sj 





26 


sl 

sj 

sj 

sj 


sj 








(cont . ' 
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Table A . 5 (cont. ) 


Event 

No. I II III IV V 


27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


No spectral observations available. 

sj 

si x/ sj si 


No spectral observations available. 


si 


x / x/-U 

>/ V 

X J si sj 

sj 


sj 

sj 


V J 

si 

si 

si 

si 

sj 

sj 


si 

sj 

J 

si 

si 

si 

si 

sj 

sj 

sj 


x/ 

si 

si 

si 

si 


sj 

si 

si 

sj 

sj 

sj 


No spectral observations available. 


si 

sj 

sj 

sj 

si sj - u 


sj 

sj 


sj 

sj 

sj 

si 


No activity recorded. 

No activity recorded. 

sj si - u si si 


sj 


si 


CONT. 


sj 


sj 


sj 

si 


si 

si 

si 


(cont 


UNCLF 


sj 


sj 


sj 
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Table A. 5 

(cent . ) 





Event 

No. 

I II 

III 

IV V 

CONT. 

UNCLF 

59 

.sj 

i 

c 

V ' 



6o 


n/ 




6i 

s/ 

s/ 




62 

No spectral 

observations 

available. 



63 

V 

s/ 




64 


s/ 





0 


APPENDIX NO. 2 


Derivation of the Electron Density Distribution Function^ Form 1 

The assumption made in the following derivation were discussed 

in Section 4.2.2 and they will be briefly reviewed here. With reference 

to Fig. 4.6, S represents the Sun and E the Earth. Exciter particles 

are ejected from a point on the Sun located at an angular distance © 

P 

from the central meridian. The particles travel along a spiral path 
with constant velocity v and produce a burst at frequency f at a point 
P. These radio waves travel along a straight trajectory to the Earth 
with a velocity equal to the velocity of light in vacuum. At a later 
time a burst at frequency f - Af is produced at point P , and it 

l 

propagates to the Earth in the same conditions as the burst at P. The 
burst from P^ is observed at the Earth with a time delay respect to 
the burst coming from P. 

We want to find a differential equation of the form 
®^+Q(r) = 0 . 

We start by using the chain law and write: 
dN = dN df dt 

dr df dt dr ’ (4-1) 

W(' wi 11 proceed to find these derivatives. 

For convenience we reproduce here the following equations : 
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df 

dt 


-io a f a 


) 


f 


3 ( > N 


,/ 2 


) 



(3 = — = constant 

c 


(4.2) 

(^. 7 ) 

(4.8) 

(4.10) 

(4.11) 


Let: 


t' = time at which the particles reach point P, 
t = time at which the radio waves emitted from point P reach 
the Earth, 

t = time at which the radio waves emitted from point P reach 
l . 1 

the Earth and 

c = velocity of light in vacuum. 


P E 

l 

v c 




Then: 


t = t ' + 


PE 


PP 


The delay in arrival at the Earth is then: 
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Let us call 


PP = ds , 

1 

PE = x , 

¥~e - pe = dx . 
1 


Then 


dt 


ds [ dx 
v c 


(A. 2) 


and 


dt _ ds dx 

dr vdr cdr (A* 3) 


From the cosine law: 


x 2 = r 2 + i 2 - 2ir cos 0 . 


Using Eq. (A. 8) we can write: 


x 2 = r 2 + l 2 - 2lr cos (9 - e j ) 


(A.b) 


so (dx/dr) can be obtained. 

Also in polar coordinates 


ds 2 = (r dQ) 2 + dr 2 


and from Eq. (4.8); 


r 2 e 2 


. / 2 


(A. 5) 


ds 


+ 1 


dr 
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Next, from Eq. (4.7) we have: 


df 

dN 



Substituting Eq. (4.7) into (4.2) we have: 


(4.6) 


df 

dt 


-10 a (9 j N l/2 ) a 


(A. 7) 


Therefore substituting Eqs. (A. 3), (A. 6), (A. 7) and (4.11) into (A.l) 
we have : 


M = . 2 !Q a (a if' 1 N (a+i)/2 ( dx 1 ds 'N 

dr c V 9 V V dr + p dr ) 


(A. 8) 


We can next perform a definite integration: 


r N(r) N - ( 0 * 0/2 = _ 2 10 a (9j) a-i 

N ( r ) 


x(r) 


x(r ) 
o 


dx + i J 


-s(r) 

s (r ) 


ds 


(A- 9) 


The quantity x is given by Eq. (A.4). From Eq. (A»5)> except 
for an integration constant: 



(A. 10) 
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The limit r^ depends on where we want to apply the boundary condition 

on the electron density N(r ) = N . 

o o 

Let us also define 


and 


P = 



(A. 11) 


D » 


10 


-a 


(a - l) (9j) 


a-l 


2/(a-i) 


(A. 12) 


Then we can rewrite the result of the integration of Eq. (A. 9) as 
follows : 


N -1 / p _ n" 1 /? 
o 




(A- 13) 


The meaning of s(r) - s(r Q ) is the length of the 

trajectory of the particles between points at distances r and r while 

o 

x(r) - x(r Q ) is the difference in paths along which the radio waves 
propagate from those points to the Earth. This equation is perfectly 
general in two dimensions in the sense that it is valid for any trajectory 
and for any form of propagation. For the case of the radial trajectory 
the limit e = 0 should be taken in Eq. (A. 5) before the integration rather 
than at Eq. (A. 10). 

Finally we can rearrange Eq. (A. 13) and write: 


N = 


D 


w 


(A. 1*0 


N~ ) + p ( v s(r) " s(r o ) ) + " x ( r D ) 


where p is defined by Eq. (A. 11), D by Eq. (A. 12), p by Eq. (A. 11), x by 
Eq . (A . 4 ) and s by Eq. (A. 5). 



APPENDIX NO. 3 


Derivation of the Frequency Drift Rate from the 
Electron Density Distribution, Form 2 

We start from Eq. (A.l) that we write now in a different form: 


dN 

df _ dr 
dt dN dt_ 
df dr 

We will find these derivaties as functions of f. 
From Eq. (4.12) 


(A. 15) 


dN 

to = " P 


(r - b) 


P+1 


(A. l6 ) 


Equating N in Eqs . (4.7) and (4.12) we have 


(r - b) P (9j) 2 


(A. 17 ) 


After solving for (r - b) and substituting in Eq. (A.l6) we can write: 


dN 

dr 


- P 


“I (p+l)/p 


9 2 j 2 A 


(A. 18) 


From Eq. (4.7) 


dN 

df 


2f 

9 2 j 2 


(A. 19) 
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Using Eq. (4.11) we write Eq. (A. 3) as follows: 


dt 

dr 


1 

c 


1. ds_ dx 

8 dr dr 


(A. 20) 


where ds/dr is given by Eq. (A. 5) and dx/dr can be derived from Eq. (A. It), 
Let us define two functions as follows: 


and 


H [r(f ) ] = 

l dr 


G [r(f ) ] = 

l dr 


(A. 21) 


and 


Then 


H [r(f ) ] = 

l 


e 2 r 2 


+ 1 


1/2 


G [r(f ) ] = 

l 


r - e rsin 

Tg - 

L F £_ 

- £ cos 

1 1 

CD 

1 

tJ 

r 2 + £ 2 - 

2£r cos £ 

*r - Tj _ 

1/2 


(A. 22) 


From Eq. (A. 17) we could express r in terms of f and then 
substitute in Eq. (A. 22). However this would make the equation too 
complicated therefore we will only indicate this dependence by defining 


two new functions : 


G(f» e, 9 p ) = G^rCf)] 


H(f, e, 9 p ) = H ^ [r(f ) ] 


and 


(A. 23) 



-265- 


Now we can write Eq. (A. 20) as follows: 


dt 

dr 


1 

c 


6 H(f, e, Q p ) + G(f, e, Q f ) 


(A.2l») 


Replacing Eqs. (A. 18). (A. 19) and (A.2U) into Eq. (A.l 5 ) we hare: 
f (p+2)/p 


df pc 

dt p 


■ 

1/p 

- 

A 9 2 j 2 


F H < f > e > S F ) + G(f, e, 9 f ) 


(A. 25) 


If we express M in (cef 3 ), r in B 0 . and f in Md take 
logarithms of the modulus in both sides of Eq. ( A . 25 ) we can write: 


1 df 
1 dt 

I = log 


2R 




-E£ L 


2R (81 10" 6 A) l/P 


( p ) log f - lo S [H(f, E, 0 p ) + 6G(f, e, 9 )] . (A. 


26) 


This is of the form: 


log 


df 


dt 


- a + a log f + A 


(A. 27) 


Except the term A, the equation is the same as Eq. (h.l). 

relationship between p and „ is the same as in Eq. (A.ll). A 

represents the departure of Eq (a Pi) -p™™ * 

q {A ' 2T) from a straight line in a log-lo^ 

plot and we could cal i D . 

old call it a correction term. It is interesting to 

note that the constant "a" depends on a •*. 

depends on the density model through A and p, 

on the mode of emission j and on the velocity 6. Qhe constant a depends 
only on the density model through p. Abe correction term A depends on 
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all the parameters: velocity, geometry of trajectory (e and 9„) model 

and mode. 

The function G 1 in Eq. (A. 22) has a singularity when the 
numerator and denominator approach zero simultaneously. This happens 
whenever the trajectory of the exciter particles passes through the 
Earth and when r = £. 

A series expansion of the numerator and denominator about r = £ 
shows that the function G^ approaches a limit whose sign depends on 
whether r approaches £ from larger or smaller values. This phenomenon 
introduces a jump discontinuity in the predicted drift rates whenever 
particles go through the Earth. 
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APPENDIX NO. 4 

Statistical Tests and Correlations Related to the Exciter Particles Velocity 
(a) y 2 Test 

We wanted to know the likelihood that the observed distributions 
of velocities (Figs, 4.23 and 4.24) came from normal distributions. We 
divided the normal curve into six categories in such a way that to each 
one there would correspond a theoretical frequency equal to or larger 
than five (Dixon and Massey, 1951 » p. 190). The degrees of freedom 
(D.F. ) turn out to be three. Figure A.l shows the percentage of areas 
in each category. The term n is the total number of cases and P is 
the probability that y 2 >_ X 2 Q -[ 3S < The estimates of the parameters were 
taken from Table 4.15. The results are shown in Table A. 6. (See Note 
No. 7 of Table A. 3.) 

Table A. 6 

Results of the X 2 Test with 3 D.F. 


Distribution 

n 

x 2 

A obs 

P 

Unambiguous 

32 

1.743 

0.62 

Ambiguous-High 

31 

1.645 

0.64 

Ambiguous-Low 

31 

2.925 

0.40 


For example for the Unambiguous distribution the test says 
that in 62 "of the cases out of 100 we should expect a fit as poor 
as this or poorer" (Mode, 1961, p. 208). According to Fisher (1958, 
p. 80), "if P is between 0.1 and 0.9 there is certainly no reason to 


- 267 - 



- 268 - 



-269- 


suspect the hypothesis tested," therefore we can accept that the three 
distributions may be samples from normal populations. The P for the 
Ambiguous -High is higher than that of the Ambiguous-Low , however this 
is not enough justification to accept it in favor of the Ambiguous -Low. 


(b) Detailed Study of the Unambiguous Distribution 

In order to study more the Unambiguous distribution we calculated 

its moments m up to the fourth (i = 4). 

The amount of skewness and flattening are given respectively by 

the parameters £ and £ where (Aitken, 1952): 

1 2 


and 


5 

i 



(sign m ) 

3 



The values obtained are: 


m = 0.320 , 

1 

m = 0 . T3U , 

2 

m = -6.0 x 10 5 

3 

and 

m = 12.0 x 10 5 

4 

Therefore: 

l = -0.91 

1 

and 

c = 

2 


2.22 . 
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K is negative because m <0. From Fig. i<. .23 it is obvious that the 

1 3 

distribution has negative skewness. The normal distribution has £ = 3 

2 

therefore our distribution is platykurtic, that is, flattened at the 
center as compared to the normal distribution. 


(c) F-Test 

Next, following the extreme alternatives of velocity chosen, we 
compared the 3 H and $ L distributions with the 3 (Unambiguous) distri- 
bution. For this we ran an F-test to test the hypothesis that the two 
distributions came from two normal parent distributions with the same 

variance (but not necessarily the same mean). Calling a 2 and a 2 

H 

the variances of the Unambiguous and Ambiguous-High distributions 
respectively, the F-test gave: 


F 



2.64 


and 


Prob (0.39 < F < 2.63) = 0.01 . 


That is, the differences in a 2 and a 2 are large enough to be considered 

n 

a rare event, it happens in less than 1% of the cases. We may conclude 
that the 3 and 3^ distributions are unlikely to be samples from the same 
normal population under the stipulated conditions. This does not 
contradict the conclusions from the y 2 test that says that each of 
distributions is likely to come from different but normal populations. 

Because the 3 L distribution departs from the 3 distribution even 
more than the 3 H we decided not to test it. 
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It is probably possible to select from the ambiguous velocities 
some high and some low in such a way as to mimic the Unambiguous 
distribution, considered as true, however we do not feel confident about 
this method of resolving the ambiguity. 

(d) Frequency of Occurrence Versus Position of Associated Flare 

In another attempt to find an indication that the Ambiguous and 
Unambiguous events could belong to different populations we studied 
the distribution of the associated flare position. The histograms of 
Fig. A. 2 illustrate the distributions. (See Note No. 7 of Table A. 3.) 

We found that 70.4$ of the 32 unambiguous cases and 64. 5 % of the ambiguous 
cases occur on the western hemisphere. The statistical parameters of 
the two distributions are shown in Table A. 7 where n is the number of 
cases. 

Table A. 7 


Parameters of the Distributions in 

F 


Distribution 

n 

9° 

y F 

0 

Unambiguous 

32 

33-8 

44.8 

Ambiguous 

31 

26.2 

55.7 


We used the F-test to test the hypothesis that the variances 
are equal. From Table A. 7 we obtained F = 1.51 and from statistical 
tables (Wine, 1964), we found: 

Prob (F > 1.6l) = 0.10 
Prob (F < 0.62) = 0.10 
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l 

I 



(a) UNAMBIGUOUS CASES (32) 


i 

i 



(b) AMBIGUOUS CASES (31) 

FIG. A. 2 LONGITUDE DISTRIBUTION OF THE H3 FLARES ASSOCIATED WITH 

TYPE III BURSTS HAVING AMBIGUOUS AND UNAMBIGUOUS VELOCITIES. 


so that at the 20$ confidence level there is no reason to doubt the 
hypothesis. By general standards this level is considered high. 

Next, we applied the Student's t-test and obtained t = 0.58. 

Then from the tables of Johnson and Jackson (1959) arid for 6l D.F. 
we find: 

Prob ( 1 1 J > 0.58) = 0.56 . 

That is, we can accept the hypothesis that the distribution in position 
of both the ambiguous and unambiguous bursts came from the same 
normal population. 

(e ) Search ' for Correlations 

In pursuit of the same goal we tried to find correlations that 
would characterize the unambiguous cases. For example in order to look 
for a possible effect of the position of the associated flare on the 
velocity we plotted B versus 0 for the unambiguous cases . • The plot is 

r 

shown in Fig. A. 3. The data seems to suggest a depletion of high 

velocities on the eastern half of the disk and a depletion of low 

velocities on the western half. Figure A.U, shows a plot of Bu 

B x versus 0 . If the suggested correlation existed it would permit the 
L r 

resolution of some cases for 0 larger than about 60 ° , however, the 
effect is not strong enough as to provide a discrimination. Besides, 
if it were, it could mean that some of our assumptions are not correct 
because the velocity of the exciter particles is unlikely to depend on 
the angular position of the Earth with respect to the associated flare. 
Figure U.27 was also studied in the search for possible correlations. 
Again, the result was not significant. 
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